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Abstract 


This report documents the derivation of minimum-dimension sets of discrete- 
coordinate and hybrid-coordinate equations of motion of a system consisting of an 
arbitrary number of hinge-connected rigid bodies assembled in tree topology. 
These equations are useful for the simulation of dynamical systems that can be 
idealized as tree-like arrangements of substructures, with each substructure con- 
sisting of either a rigid body or a collection of elastically interconnected rigid 
bodies restricted to small relative rotations at each connection. Thus, some of the 
substructures represent elastic bodies subjected to small strains or local deforma- 
tions, but possibly large gross deformations; in the hybrid formulation, distributed 
coordinates, referred to herein as large-deformation modal coordinates, are used 
for the deformations of these substructures. The equations are in a form suitable 
for incorporation into one or more computer programs to be used as multipurpose 
tools in the simulation of spacecraft and other complex electromechanical systems. 



Large- Deformation Modal Coordinates 
for Nonrigid Vehicle Dynamics 


I. Introduction 

A. Background and Motivation 

Although the influence of vehicle nonrigidity on spacecraft attitude stability and 
control has been noted in singular cases since the very beginning of the space age 
(Refs. 1, 2), it is only in recent years that this influence has been widely recognized 
as a systems problem of paramount importance to the success of many missions 
considered crucial to the space program. 

The determination of the influence of spacecraft nonrigidity on mission perform- 
ance is not a single problem but a family of problems. It is therefore appropriate 
that the engineer have at his disposal an arsenal of analytical and computational 
weapons, and not just a single approach or a single computer program. This report 
focuses accordingly on the exposition of a particular method for simulating a cer- 
tain class of spacecraft subject to large changes in attitude and configuration. The 
formulation provided here seems to be more comprehensive than any previously 
published for application to spacecraft amenable to mathematical modeling as a 
collection of point-connected rigid bodies assembled in a topological tree (with 
contiguous bodies sharing a common point and no closed rings of bodies); but not 
every spacecraft should be modeled in this way. In order to establish the position 
of this contribution in the spectrum of methods currently available for the simula- 
tion of nonrigid spacecraft, the literature in this field must be reviewed and the 
need for extensions of existing procedures established. 
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1. Energy sink and reaction force methods. Attitude anomalies observed after 
the launch of Explorer I (Fig. 1) were explained (Ref. 1) by means of an approxi- 
mate dynamic analysis in which the vehicle nonrigidity was ignored in generating 
a preliminary representation of rotational motion, and then the oscillatory deforma- 
tions of flexible appendages on the vehicle (the wire turnstile antennas in Fig. 1) 
were estimated by prescribing the appendage base motion to be that obtained from 
the rigid vehicle approximation. Estimates of the energy dissipation rate were then 
obtained from the predictions of relative motions of the appendage components, 
and finally, the kinetic energy of the original rigid-body model of the vehicle was 
reduced over time in correspondence with the estimated rate of dissipation of 
energy in the nonrigid portions of the vehicle. In this fashion, it was determined 
in Ref. 1 that nonrigidity and consequent energy dissipation would cause Explorer I 
to change asymptotically from initial spin about its longitudinal axis to eventual 
spin about a transverse axis, thus explaining the observed behavior of this satellite. 


Since the nonrigidity of a spacecraft is acknowledged by the analyst who follows 
the practices of Ref. I only in the process of predicting the energy dissipation rate, 
the vehicle mathematical model which he adopts is often characterized by the 
necessarily imprecise term quasirigid, which implies that over a “short” time inter- 
val, the overall vehicle motion is essentially that of a rigid body. The procedure 
described here has been widely employed in the dynamic analysis of spinning and 
dual-spin spacecraft (Refs. 3-5), and it is known in its various guises as the energy- 
sink approach. In a variation of this procedure applicable also to systems of quasi- 
rigid bodies, the energy loss calculation is replaced by an estimate of the reaction 
forces on the spacecraft due to the relative motion of components (Refs. 6, 7). 


Although the energy-sink and reaction-force methods have proven extremely 
useful, particularly in the interpretation of observed attitude anomalies (Refs. 1, 5) 
and in preliminary design (Refs. 3, 4, 6, 7), they are inherently approximate, and 
their predictions generally require more formal corroboration. In some industrial 
applications, an energy-sink or reaction-force analysis might suffice to establish, in 
a given case, the noncritical nature of vehicle nonrigidity, but in this study, a deeper 
investigation is presumed necessary. At issue is the question of how to proceed. 


2. Discrete-coordinate methods. The practice of modeling a spacecraft as a 
collection of discrete rigid bodies (or particles and rigid bodies) is limited only by 
the skill with which the analyst can devise his model, the size of the computer 
(and computer budget) available to integrate numerically the resulting ordinary 
differential equations of motion in the discrete coordinates of the system, and the 
accuracy and computational stability of his integration program. Since there need 
be no limitations on the relative motions of the constituent parts of the spacecraft, 
this discrete-coordinate approach is conceptually ideal for spacecraft undergoing 
large changes in orientation and configuration, such as the dual-spin vehicle shown 
in Fig. 2 in the process of deploying solar panels while maintaining antenna point- 
ing control. This vehicle would seem to require at least eleven bodies in the mathe- 
matical model, and more bodies would be required if it became necessary to 
accommodate the flexibility of the individual solar panels or the antenna. Addi- 
tional problems might be presented by fluids or flexible components in the rotor. 
Practical limitations in this approach arise because the costs of computer simula- 
tion increase rapidly with the number of bodies in the system, and problems of 
computational inaccuracy and instability increase correspondingly. 
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The objective of formulating equations of motion in a form most conveniently 
obtained and most readily solved has preoccupied dynamicists for centuries. The 
accomplishments of Euler, Lagrange, and Hamilton in this area are well known; 
Lagrange found, for example, that if he restricted the dynamical system to one 
involving only holonomic constraints, and if all forces other than those imposed to 
satisfy constraints were conservative, then he could write scalar equations of 
motion in the form 


d /8L\ _ d L 
dt\dqi) dqi 


i = 1, • • • , n 


where the Lagrangian L is the difference in kinetic and potential energies, and 
q u ■ ■ ■ , q n are the generalized coordinates of the system. Similarly (but less spec- 
tacularly), modem dynamicists have discovered that if they impose other sets of 
restrictive conditions on the dynamical system, they too can obtain generic equa- 
tions of motion in a form particularly convenient for practical application; in a 
modern context, the equations are often judged by their suitability for digital 
computer programs for numerical integration. 

The specific dynamical system for which equations of motion have most success- 
fully been formulated generically consists of a collection of rigid bodies which are 
point-connected (i.e., interconnected in such a way that each body shares at least 
one point in common with some other body), with no closed loops or rings of point- 
connected bodies in the system. This system has become known as the set of point- 
connected rigid bodies in a topological tree. 


The eleven-body model suggested for the dual-spin spacecraft in Fig. 2 provides 
a natural example of such a system of point-connected rigid bodies in a topological 
tree. This is a special case because contiguous bodies in every case share a common 
line and not just a common point ; such systems are called hinge-connected. If it is 
deemed essential that the dynamical model accommodate the flexibility of each 
of the six individual segments of the solar panels in Fig. 2, then this can be accom- 
plished within the framework of the point-connected topological tree model in a 
variety of ways, as illustrated by Figs. 3a, b, c. In each model, a series of line hinges 
has been introduced, with suitable restraints on rotation (spring torques and damp- 
ing torques). In model (a), only the bending flexibility of the panel is accommo- 
dated. Model (b) also permits the panel to twist about its central axis. Model (c) 
permits twisting and bending about two orthogonal axes. With the introduction 
of another family of hinges normal to the page, one could accommodate in-plane 
deformations of the panels. 

Should it become necessary to connect rigid bodies with linear springs in addi- 
tion to the basic point connections, this can be accomplished within the framework 
of the point-connected topological tree system equations, even if the linear springs 
provide the closure of a chain of bodies into a ring; the linear spring forces must 
then be included with the external forces in the system equations. 


Note that in Figs. 2 and 3, the bodies are hinge-connected rather than merely 
point-connected. It is always possible to substitute a hinge-connected system for 
a point-connected system simply by introducing auxiliary bodies, such as the small 
blocks shown in Figs. 3b and 3c; such bodies can be idealized as massless and 
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Fig. 3. Rigid-body model of flexible panel segments from Fig. 2 


dimensionless, and their introduction does not alter the minimum dimension of 
the problem because it does not change the number of degrees of freedom in 
the system. 

The disadvantages of employing discrete-coordinate equations of motion of a 
system of point-connected (or hinge-connected) rigid bodies in a topological tree 
are obvious: (1) No degree of sophistication in mathematical modeling can provide 
proper representation of actual local deformations, as required for the calculation 
of stresses; and (2) the number of rigid bodies in the system soon grows too large 
for computer simulation. (If model 3c is used for each of the six panels in Fig. 2, 
the vehicle model involves 95 hinge-connected bodies!) 
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The advantages of employing discrete-coordinate equations of motion of a sys- 
tem of point-connected (or hinge-connected) rigid bodies in a topological tree can 
also be clearly stated: (1) A minimum-dimension set of equations of motion can 
easily be formulated generically for an n-body system and programmed once and 
for all for reasonably efficient digital computer numerical integration; and (2) rela- 
tive rotations between connected bodies can be unrestricted in magnitude, per- 
mitting the simulation of time-varying configurations or large deformations. 

The contributions of the present report in the area of discrete-coordinate formu- 
lations are based upon a very substantial foundation provided by others, whose 
work is briefly reviewed in what follows. 

The publication in 1963 by Fletcher, Rongved, and Yu (Ref. 8) of the equations 
of motion of two point-connected bodies representing a gravity-stabilized satellite 
stimulated the parallel derivations by Hooker and Margulies (Ref. 9) and by 
Roberson and Wittenburg (Ref. 10) of the corresponding equations for an arbitrary 
number of point-connected rigid bodies in a topological tree. Reference 9 provides 
vector-dyadic equations, and Ref. 10 contains equivalent but independently derived 
matrix equations. Both sets of equations involve all torques transmitted between 
bodies at connection points, and this can be a substantial liability if relative motions 
are constrained (as by a line hinge). As a consequence of their retention of unknown 
kinematical constraint torques and the equations necessary to provide solutions for 
these unknowns, neither Ref. 9 nor Ref. 10 offers a set of equations of motion of the 
minimum dimension required by the number of kinematical degrees of freedom 
in the system. 

The equations provided in Refs. 9 and 10 fulfilled in a timely manner a keenly 
felt industrial need, and several aerospace organizations developed computer pro- 
grams based on these equations or others of similar character derived indepen- 
dently. Procedures were proposed and implemented by Velman (Ref. 11) and by 
Fleischer (Ref. 12) to enable the computer to integrate the equations without actu- 
ally solving for the constraint torques. Russell (Ref. 13) devised an approach to 
deriving the equations without ever introducing the constraint torques; he formu- 
lated vector equations of motion for subsets of rigid bodies selected so that only 
one joint connected each subset to a body outside of the subset, and then dot- 
multiplied each subset equation by unit vectors dictated by the degrees of freedom 
at that one joint of the subset linked to an external body. Russell provided a pro- 
cedure for systematically deriving equations rather than a generic system of equa- 
tions suitable for complete pre-programming. Hooker (Ref. 14) then indicated in 
a brief note a procedure for operating on the Hooker-Margulies equations (Ref. 9) 
in such a way as to obtain a minimum-dimension set of equations. Hooker’s sug- 
gestion required coordinate transformations, summation of selected subsets of 
equations corresponding to those used by Russell, and dot-multiplying by unit 
vectors established by the details of the joint geometry. In his published note 
(Ref. 14), these operations are outlined but not presented in explicit detail. 

Computer programs based on generic formulation of equations of motion of 
discrete-coordinate systems have been widely used in the aerospace community 
since 1965, and in a few cases (Refs. 11, 15, 16), public documentation is available. 
Although these programs have served a useful and even necessary function, they 
have not become a panacea for all problems of simulation of nonrigid spacecraft, 
because of the disadvantages noted previously. The deficiencies of the rigid-body 
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models in representing local strains and stresses are rarely of concern to the atti- 
tude control specialists who have used this approach, but the prohibitive cost of 
computer simulations involving many bodies has effectively restricted the use of 
these programs. It soon became evident that more efficient approaches would be 
required. 


3. Vehicle normal-coordinate methods. The problem of accurate representation 
of the vibrations of a complex elastic structure with a modest number of coordi- 
nates has traditionally been resolved among structural dynamicists by the use of 
distributed or modal coordinates. When such coordinates are chosen so that each 
one represents a normal mode of vibration, in which the structure can vibrate 
without the excitation of other modal coordinates, then the modal coordinate 
equations are uncoupled, and one can with impunity solve the scalar equations of 
motion one at a time. Estimates of the motion can then be based on the super- 
position of as many (or as few) of the modal coordinate solutions as may be desired. 
This approach, called the vehicle normal-coordinate method, is rarely appropriate 
for vehicles of the complexity of modern spacecraft, however, and in recent years, 
other alternatives have been developed which combine some of the computational 
efficiencies of the vehicle normal-coordinate method with some of the advantages 
of generality offered by the discrete-coordinate method by using combinations 
of discrete and distributed coordinates. Any approach which combines discrete 
coordinates describing the translations and rotations of some bodies or reference 
frames of the system with distributed or modal coordinates describing the small 
relative motions of other parts of the system is called a hybrid-coordinate method. 


4. Hybrid- coordinate methods. Hybrid-coordinate methods are many and varied, 
being represented, for example, by Refs. 17-21. The last of these references pro- 
vides an extensive discussion of each of the three approaches (discrete coordinate, 
vehicle normal coordinate, and hybrid coordinate) in the stage to which these 
methods had evolved by 1970. The application to complex spacecraft of the hybrid 
coordinate method described in Ref. 21 is documented in Refs. 22-24, and an exten- 
sion of Ref. 21 providing for the representation of a flexible appendage with a 
distributed-mass finite element model is given in Ref. 25. Quantitative compari- 
sons of simulation results and computer time for discrete-coordinate and hybrid- 
coordinate simulations of a Viking spacecraft may be found in Ref. 16, and similar 
comparisons of discrete-coordinate and vehicle normal-coordinate simulations are 
available for the Radio Astronomy Explorer satellite in Ref. 15. 


Because distributed coordinates are used to describe the deformations of flexible 
substructures in the hybrid coordinate approach, the internal mathematical model 
of the flexible structure is obscured; it might be an elastic continuum (Ref. 18), or 
a collection of particles interconnected by massless elastic elements (Refs. 22-24), 
or a collection of interconnected elastic elements possessing mass (Ref. 25). The 
question might even be ignored in formulating the equations (as in Ref. 20), but 
when the time comes for applications, some decision must be made. In the formula- 
tions which have been published thus far, no one has suggested the adoption for 
an elastic substructure of a model consisting of a point-connected set of rigid 
bodies in a topological tree, but, as suggested by the flexible solar panel segment 
models in Fig. 3, this is a realistic option. The advantages associated with this 
choice (to be amplified later) stem from the possibility of beginning with a descrip- 
tion in terms of discrete coordinates describing relative rotations of contiguous 
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bodies, and therefore representing local deformations. This distinction permits the 
use of modal coordinates to represent large gross deformations. 


The several variations of the hybrid-coordinate approach which have emerged 
thus far share certain limitations. (1) The critical task of modal analysis of flexible 
appendages has in all publications to date been based on restricted boundary con- 
ditions; in the most general formulation (Ref. 25), the flexible substructure is 
assumed to be attached to one rigid body with arbitrary motion, and in other 
formulations in which the flexible substructure is connected to more than one 
additional substructure, the system is assumed to be nominally inertially stationary. 

(2) Flexible appendages are restricted to small displacements from a reference 
state, implying small deformations both in a local sense and in an overall sense. 

(3) Except in restricted cases (e.g., Ref. 20), hybrid-coordinate formulations have 
not been incorporated into comprehensive generic computer programs (although 
such programs are under development at JPL and elsewhere). 

5. The need for new procedures. When one considers the limitations of the 
discrete-coordinate, vehicle normal-coordinate, and hybrid-coordinate procedures 
in their present stages of development, it becomes clear that none of the methods 
described in the references is suitable for certain classes of vehicles. As an example, 
note that a flexible substructure experiencing small strains and large overall defor- 
mations can at present be accommodated only by a discrete-coordinate formula- 
tion, which must be based on a crude model to avoid prohibitive expense (Ref. 15). 
As an illustration of a different class of problem beyond the scope of present 
methods, note that none of the published formulations will suffice in the accurate 
simulation of the dual-spin vehicle shown in Fig. 2 while its solar panels are being 
deployed, unless the panel segments are presumed rigid. The discrete-coordinate 
approach is ideal for a gross approximation of the deployment operation (using 
rigid panel segments), and the hybrid coordinate approach in Ref. 21 is ideal for 
the simulation of flexible panels after deployment. But if antenna pointing con- 
trol requirements during deployment demand the accommodation of panel seg- 
ment flexibility in the mathematical model, then the discrete-coordinate approach 
becomes computationally cumbersome and perhaps prohibitively expensive, while 
the existing hybrid-coordinate approaches are inapplicable because of theoretical 
restrictions. New procedures are therefore required. 


B. Scope of Study 

The objectives of this report include the documentation in scalar detail not 
available elsewhere of a minimum-dimension set of discrete-coordinate equations 
of motion of a hinge-connected set of rigid bodies in a topological tree, and the 
introduction of a new procedure for modifying these equations by coordinate 
transformation and truncation, so as to obtain a hybrid-coordinate formulation of 
the system equations of motion, employing large deformation modal coordinates 
in combination with discrete coordinates. 


The discrete-coordinate equations presented here (see Section II) are the results 
of the application of a variation of Hooker’s procedure (Ref. 14) in a systematic 
fashion that results in generic matrix (or scalar) equations, formulated in the spirit 
of the Roberson-Wittenburg equations (Ref. 10) but with constraint torques 
removed. 
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The hybrid-coordinate equations presented here (see Section III) have been 
obtained by (1) isolating discrete-coordinate equations for subsets of rigid bodies 
identified collectively as a flexible appendage or flexible substructure, (2) lineariz- 
ing in the variables which describe the relative angular motions of contiguous 
bodies within the flexible substructure, (3) transforming these substructure vari- 
ables into distributed coordinates for the substructure, and (4) truncating these 
distributed coordinates. As noted in detail in Section III, the linearization employed 
here remains a useful engineering approximation as long as contiguous bodies 
within the flexible substructure experience “small” relative rotations, even if the 
rotations between the first and last bodies in a chain become “large” (say, 90 deg). 
This is not the case for traditional “small deflection theory” as applied to beams, 
etc., because this theory assumes small relative displacements and rotations of all 
differential elements. The modal coordinates in Section III are useful even when 
the actual elastic substructure experiences large elastic deflections from a reference 
state as long as the elongations and shear deformations of differential elements 
remain small. 


The useful results in this report consist of three versions of the equations of 
motion of an arbitrary vehicle modeled as a collection of point-connected rigid 
bodies in a topological tree: (1) unrestricted discrete-coordinate equations (see 
Eq. 1), (2) partially linearized discrete-coordinate equations (see Eq. 30), and 
(3) hybrid-coordinate equations (see Eq. 118). These equations can all be incor- 
porated as subroutines in a single computer program, or they can become mem- 
bers of a family of computer programs, with use of the hybrid-coordinate program 
or subroutine requiring additional input from a separate program designed for 
eigenvalue-eigenvector analysis of linear subsets of differential equations. 


II. Discrete-Coordinate Equations of Motion 

A. Mathematical Model 

Figure 4 illustrates a set of eleven hinge-connected rigid bodies, labeled accord- 
ing to a convention designed to facilitate the processing of the equations of motion. 
This example, which portrays a system of a degree of complexity that would mini- 
mally suffice for the spacecraft shown in Fig. 2, will be useful in illustrating and 
interpreting the general mathematical model and the labeling conventions adopted 
here for its description. (In reading the lengthy list of conventions and symbol 
definitions, frequent reference should be made to Fig. 4.) Only after entertaining 
the complete cast of symbols used in what follows will consideration be given to 
the much smaller list of symbols to which values must be assigned by the user of 
the projected computer program. 


The mathematical model consists of a set of n + 1 rigid bodies interconnected 
by n hinges*; the indicated numbers of bodies and hinges imply a tree topology. 
Any interbody connection forces in addition to those at the n hinges must be 
treated as forces external to the entire system. (If there were a linear spring con- 
necting the mid-point of body in Fig. 4 to some point on body / 0 , for example, 
it would be necessary to replace the spring with a pair of equal and opposite forces 


*The word “hinge” as used here implies a connection which maintains a line common to both 
bodies; such a connection is sometimes called a line hinge. 
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Fig. 4. Eleven rigid bodies interconnected by ten hinges 


applied to the spring attachment points on 4 e and 4 0 ; the magnitude and direction 
of these forces would have to be expressed in terms of the unknown kinematical 
variables of the system.) 

1. Definitions of fundamental symbols 

Def. 1. Let n be the number of hinges interconnecting a set of n + 1 bodies. 

Def. 2. Define the integer set £8 = {0, 1, • • • , n). 

Def. 3. Define the integer set & = {!,••• >»}• 

Def. 4. Let 4 0 be a label assigned to one body chosen arbitrarily as a reference 
body, and let 4 U ■ ■ - , 4„ be labels assigned to the remaining bodies in 
such a way that if is located between 4 0 and then 0 < / < k. 

Def. 5. Define the kth neighbor set £8* for k e £8 such that r e £8* if 4 r is attached 
to 4 k . 

Def. 6. Define the dextral, orthogonal sets of unit vectors b^, b\, b* so as to be 
imbedded in 4>. for k e £8, and such that in some arbitrarily selected nomi- 
nal configuration of the system b* = b> for a — 1, 2 , 3 and k, j € £8. 
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Def. 7. Define {b 4 } as the column array of unit vectors* 



for ke £B. 

Def. 8. Define {i} as a column array of inertially fixed, dextral, orthogonal unit 
vectors i x , i 2 , i 3> so that 



Def. 9. Let C be the direction cosine matrix defined by 

{b 0 } = C {i} 

A 

Def. 10. Let to 0 = {b°} T <o° be the inertial angular velocity vector of so that w° 
is the corresponding 3 by 1 matrix in vector basis {b 0 }. 

Def. 11. Let c k be the mass center of 4 k for fee £3 (see Fig. 4). 

Def. 12. Let /i k be a point on the hinge axis common to 6 k and for j < k and 
keP (see Fig. 4). 

Def. 13. For j e £B kl let p ki = {b k } T p k> be the position vector with respect to c k of 
the labeled point on the hinge axis (either fij or /i k ) which is common to 
4 k and 4j, so that p 4j is the corresponding 3 by 1 matrix in vector basis 
{b 4 }. (See Fig. 4, and note the special cases p 21 = p 23 = p 32 = p 34 = p 43 = 0.) 

A 

Def. 14. Let g 4 = {b k } T g k be a unit vector parallel to the hinge axis through y lk , so 
that g k is the corresponding 3 by 1 matrix in vector basis {b 4 }. 

Def. 15. For keff 1 , let y k be the angle of a g 4 -rotation of 4 k with respect to the body 
attached at /t k , that is, a rotation during which a right-handed screw fixed 
in 4 k with its axis parallel to g 4 advances in the direction of g k . Assign the 
value zero to y k for ke£P when the system is in its arbitrarily chosen nomi- 
nal configuration, for which b£ = b£, a = 1, 2, 3; k, je£B. 

A 

Def. 16. Let X = {i} T X be the position vector of the system mass center CM with 
respect to an inertially fixed point $, so that X is the corresponding 3 by 1 
matriv in an inertial vector basis {i}. 

Def. 17. Let m k be the mass of 4 k for keff}. 

Def. 18. Let I* = {b 4 } 7 !* {b 4 } be the inertia dyadic of 4 k for c k , so that Z 4 is the 
corresponding inertia matrix in vector basis {b 4 }. 


*Note that braces { } are used in this report both to identify column arrays of vectors and to 
enclose sets of integers; the distinction is apparent as soon as the objects within the braces are 
identified as vectors (boldface) or scalars. 
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A 

Def. 19. Let F* = {b k } T F k be the resultant vector of all forces applied to except 
for those due to interbody forces transmitted at hinge connections, so 
that F k is the corresponding 3 by 1 matrix in vector basis {b*}. (Thus, any 
interbody forces due to spring connections, etc., which are separate from 
the n hinge connections would contribute to F*, as would all forces from 
sources external to the system.) 

Def. 20. Let T k = {b*'} r T s be the resultant moment vector with respect to c k of 
all forces applied to 4 except for those due to interbody forces trans- 
mitted at hinge connections, so that T k is the corresponding 3 by 1 matrix 
in vector basis {b*}. 

Def. 21. Let r k be the scalar magnitude of the torque component applied to in 
the direction of g k by the body connected at /t k . 

Def. 22. Let the Kronecker delta symbol 8 a p be defined by 

S af} = 1 - \ (a - p) 2 [5 - (a - /?) 2 ] <*,/? = 1 , 2,3 

so that 

( lif« = j8 

= i 

l 0 otherwise 

Def. 23. Let the Levi-Civita density symbol € 0 p» be defined by 


{a — P) (/? — 6) (8 — a) a, fl, 6 1, 2, 3 

so that 

! + 1 for a, p, 6 a cyclic permutation of 1, 2, 3 

0 for a ~ p or p = 0 or 6 = a 

— 1 for a, p, 6 a cyclic permutation of 1, 3, 2 

Def. 24. Let the tilde symbol (~) signify in application to a 3 by 1 matrix V with 
elements V„ (9 = 1,2,3) transformation to a skew-symmetric 3 by 3 
matrix V, whose elements are given by* 

Thus, the matrix V has the expanded form 



-V, Vf 

o -V, 

V: o. 


2. Definitions of derived symbols 

Def. 25. Let F = (i} T F = 2 F*, so that F is the 3 by 1 matrix corresponding to 

I Mffi 

F in vector basis (i). 


12 


*When lowercase Greek subscripts appear twice in a single term, summation over the values 
1, 2, 3 is implied, so that, for example, ««.sV 0 = e 01( jV, + t nSt V, + t n ^V,. 
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Def. 26. For r e (£3 — k)* let iV kr denote the index of the body attached to and 

A A 

on the path leading to 6, and let N ** = k and N fc = N k0 (to simplify nota- 
tion). The term network elements will be applied to the scalars typified 
by N kr , which are (n + l) 2 in number. 

Def. 27. Define the scalar e sk such that for keffi and seff, 

1 if /i s lies between / 0 and / k 

0 otherwise 

The term path elements will be applied to the scalars E 8k , which are 
n (n + 1) in number. 

Def. 28. Let 3i k j be the branch set of integers r such that r e £8« if j — N, cr . Thus, 
£Bkj consists of the indices of those bodies attached to on a branch 
which begins with 4j. Since there are n hinges, there are 2n branch sets 
appearing in pairs, such that £3/* = ffi — ffikj- 

Def. 29. Let C ri be the direction cosine matrix defined by 

{b r } = C r * {V} 

for r, /' e£3. In the nominal state, y* = 0 for keff, and all unit vectors of 
corresponding subscript are aligned, so that C ri then is U, the 3 by 3 unit 
matrix. 

A 

Def. 30. Define Q7l — 2 so that Offt is the total system mass. 

it® 

Def. 31 . Define 0 7l k j = 2 m r , so that Q Tlkj is the total mass on the branch of the 
system attached to 6 k and commencing with 6 f . 

Def. 32. For re {£3 — k}, define h kr = p k!itr , and define Li k — 0. Thus, L fcr is the 
vector from c* to the hinge point attached to d k and on the path leading 
from 6 k toward S T . (In Fig. 4, for example, L 56 = L 57 = L 58 = p 56 and 
L 3r = 0 for re £8.) 

Def. 33. Define (b k } T D kk = D kk = — 2 k. k ’mj/Qryi for keffi. 

m 

Def. 34. Let b k be a point called the barycenter of 6 k and fixed relative to such 
that D kk is the position vector of c k with respect to b k . 

Def. 35. Define [h k ) T D ki — D*’ = + L, k ’ for k, feffi. Thus, is the vector 

from bit to the hinge point fixed in S k and on the path leading to Sj. 

Def. 36. Define the dyadic 

V* = I* + m k (D kk • D**U - D kk D kk ) + 2 QW*, (D* 7 • D*>U - D k ’D k ’) 

*‘For notational brevity, the set SB — {k} is designated S3 — h. 
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where U is the unit dyadic, and define the corresponding constant 3 by 3 
matrix in vector basis {b k } as 

= {b*} - - {b^} 7, 

= I k + m k (D kkT D kk U - D kk D kkT ) + 2 Qflkj (D k > T D k >U - D k ’D k ’ T ) 

it83k 

Def. 37. For / e {£8 — k), define the dyadic 

<j>fcj 4 —Cyyi (D j7c ‘D^U - D^D*’) 
and define a corresponding variable 3 by 3 matrix as 

$*>' = {b J } *0*' • {b*} r 

= -0M [ {b , ‘} • (D> kT {b>} • {b*=} r D k >) {b*} r 
- {bO • {b >} T D ik D k ’ T {b*} • {b*} 11 ] 

= ~7m [{b’} • (D ikT C ik D ki ) {Wy - D’ k D k ’ T ] 

= -cm (C ik D ikT C’ k D ki - D ,k D kjT ) 

Def. 38. Define the 3 by 3 matrix 

a 00 = 2 2 C°»4W = 2 C ok & k C kB 

kiS3 m ki® 

-cm 2 2 (C ok D ikT C’ k D ki - C°’D’ k D k ’ r ) C k0 

ki£D /effi-k 

Def. 39. For ke£P, define the 3 by 1 matrices 

flo* = 2 2 e kr C°WC rk g* = 2 skrC 0r <s> rr C r Y 

re® je® riff 

- cm 2 2 (C 0r D’ rT C’ r D ri - C°’D ir D r ’ T ) C rt g* 

riff jlSi-r 

and let a k0 = <&- 

Def. 40. For /, k e <£?, define the scalars 

a hi = f22 e kr e is C k WC'igt 

n£P st£P 

= g kT 2 e kr e ir C kr <f> rr C ri g 

riff 

-Q7l^ T 2 2 <*r£j S (C kr D srT C sr D rs - C ks D sr D r,T ) C r > g j 

riff 3 iff-r 


14 
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Def. 41. For ke(B, define the 3 by 1 matrices* 

A k =T k + 2 D ki C k ’F> - & k 2 *rkyr (C ka Z>C 0r + 2 e, r y.C“ g’C") g r 

jt(B riff siff 

- [C'W* + 2 B ik yiC k >WC*] & k [C fc V + 2 WiCVgj] 

jiff iiff 

+ Qfl 2 b ki C ki [C>°tfC oi + 2 e ri yrC ir g r C ri ] ! D ik 

ji£R-k riff 

- 2 C k ”t> ki 2 trjyr [C k WC 0r + 2 e sr y a C k °g‘C° r ] g r 

jt£P re£P st£? 

where the superscript 2 indicates the square of the matrix in square 
brackets. This expression for A* can be written in a form which reveals 
more explicitly its dependence on time-variable quantities by substituting 
in A k the expansion of $ k ’ recorded in Def. 37. The resulting expression is 

A h = T k + 2 D ki C ki F> - <t>™ 2 £rky r ( C ko Z°C 0r + 2 e sr y s C ks g , C sr ) g r 

ilffi riff si ff 

- [C ko tfC ok + 2 e ikyiC k ’g’C’ k ] [C So (b° + 2 *ik 7 jC w g’] 

jiff jiff 

+ (ftl 2 {D k ’C k > [C’^C 01 + 2 e r jy T C ir g r C ri ] 2 D’ k 

iiffj - k riff 

+ (UD ikT C’ k D ki - C k) D’ k D kiT ) 2 tnyr [C^C 0 ' + 2 s,r7»C*‘g»C 8 '] g r ) 

rtd> $e£P 

3. Augmented bodies and barycenters. The barycenter b k of body 6 k is defined 
in Def. 34 in mathematical terms which admit an interesting interpretation as a 
physical property of an abstraction called the augmented body. If one imagines 
the body 6 k to be augmented by the addition of a particle of mass Q7l ki at the con- 
nection point of 6 k with 6, (for ft(Bk), then the result is called the kth augmented 
body, here referred to as 6 ' k . Since 0 7l k j is the total mass of the branch of bodies 
attached to 6k commencing with 6j, the mass of each augmented body is 077, the 
total system mass. 

Definition 33 provides 


-OWD** = 2 mjU’ 

je£B 

which reveals the fact that — D** is the position vector of the mass center of 6 k with 
respect to the mass center c k of b k . Thus, the barycenter b k is the mass center of the 
kth augmented body 6 k . With this interpretation comes the relationship 

2 m r D* r = m^D** + 2 m r D* tr = 0 

r«SJ re® - It 

Since D' tfe is the position vector from the barycenter b k to the mass center c k of 6k, 
and (from Def. 35) D k > is the position vector from b k to a particle of mass Q 7l k j 
located at the connection point of 6 k leading to body 6j, for / e (Bk, one can interpret 
4> kk from Def. 36 as the inertia dyadic of the augmented body 6 k with respect to the 


*It may be helpful to note that identifying superscripts are never attached to the right of a scalar 
in this report; a superscript on a scalar denotes its exponent (so that a b kj means the scalar a kj 
raised to the fifth power). In contrast, nonscalar quantities are not exponentiated unless noted 
explicitly (as in Def. 41), so that a superscript to the right of a matrix, vector, or dyadic is then 
an identifier and not an exponent. 


JPL TECHNICAL REPORT 32-1565 


15 




barycenter b k . Such physical interpretation is not possible for the dyadic defined 
in Def . 37, since this dyadic involves vector D jfc fixed in and vector D s ’ fixed in Sk- 

in illustration of the concepts just described. Fig. 5 shows the augmented body 
So associated with body S Q in the system of Fig. 4. 

B. Matrix Dynamical Equations 

In terms of the symbols defined in the preceding section, a minimum-dimension 
set of n + 3 scalar equations of rotational motion of a set of n + 1 rigid bodies inter- 
connected by n line hinges can be written as: 


000 

001 

002 

0On 


'<1°' 


' 2 C ok A k 

kegS 

0io 

011 

012 

0in 


Vi 


g lT 2 «i*C‘*A* + ti 

ke& 

020 

021 

022 

02 n 


72 

— 

g 2T 2 ertC*A* + t 2 

k *0 

0»o 

0» 1 

0»2 

0nn 


Y» 


g nT 2 enkC nk A k + Tj\ 
keg 


This result is proven in Appendix A. 

Note that the coefficient matrix on the left side of Eq. (1) is, from Defs. 38-40, 
a real, nonsingular,* symmetric matrix which varies with time as the kinematical 


’Singularity would imply that at least two rows (or columns) of the matrix are dependent, in 
which event some linear combination of the indicated scalar equations would provide a matrix 
equation with a zero row (or column) in the coefficient matrix of the column matrix of the most 
highly differentiated variables. Thus, singularity would indicate an incomplete set of equations 
unless the variables are constrained. Since the variables in Eq. (1) are independent, the coeffi- 
cient matrix on the left must be nonsingular. 
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variables y x , ■ • • , y„ change in value. If these angles and their time derivatives 
remain in the immediate neighborhood of their nominal zero state, and o>° remains 
in the immediate neighborhood of its nominal constant value, then upon linear- 
ization in the variational coordinates, the indicated coefficient matrix becomes a 
constant, symmetric, nonsingular matrix, which can be inverted in advance of 
numerical integration of Eq. (1). 

If for preliminary analysis one wishes to replace the variable y > (/' e($) by a 
known function of time, then one can remove y, from the angular acceleration 
matrix on the left side of Eq. (1), removing also the corresponding row and column 
of its coefficient matrix, and then rewrite on the right side those terms just deleted 
by the removal of the column on the left. The result is a set of n + 2 scalar differen- 
tial equations in which time appears explicitly. 

The equations of mass center translation of the total system can from Newton’s 
second law and Defs. 3, 16, 25, and 30 be written as the vector equation 

{iy f = (yn {i} T x 

or the matrix equivalent 

F = (VhX (2) 

Equations (1) and (2) constitute a complete set of dynamical equations, but they 
are not fully descriptive of the system motion until they are augmented by control 
equations specifying r x , • • , r„ and certain external forces and torques, and aug- 
mented also by kinematical equations as provided in the next section. 

C. Matrix Kinematical Equations 

The kinematical variables adopted in the previous sections are as follows: y k for 
ke& (Def. 15); C ri for r,jeQ 3 (Def. 29); C (Def. 9); a. 0 = {b 0 } -co° (Def. 10); and 

A 

X = (i) *X (Def. 16). Although the equations of motion have been expressed in 
terms of these quantities, the latter are not all independent. Relationships among 
kinematical variables developed in this section must therefore either be considered 
in conjunction with the dynamical equations or be substituted into them to remove 
redundant variables whenever a solution is sought. 

The direction cosine matrix C (Def. 9) which establishes the inertial attitude of 
4 is related to the inertial angular velocity matrix o>° of <f 0 (Def. 10) by 

C = -S°C (3) 

A • 

The relationship <3° = CC T is used here to define the angular velocity to 0 , as on p. 96 
of Ref. 26 (where the symbol C is the transpose of that used here). 

The direction cosine matrix C ri (Def. 29) relates sets of orthogonal unit vectors 
fixed in and and hence depends upon those angles y a for which /i a lies 
between 6 r and />, and also upon the corresponding unit vectors g° defining the 
intervening hinge axes. For the special case in which 6 r and are contiguous, it 
is always possible to identify them in one sequence or another as k and N k (as 
introduced in Def. 26), and then to express C kN ‘ and C Nkk in terms of the single 
angle y k and the single matrix g\ as follows. 
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The orientation of with respect to is established by a simple right-handed 
rotation through angle yk about a hinge axis with positive direction of rotation 
established by g k = {b k } T g k . As will be demonstrated, the direction cosine matrices 
then become 

C kNt = U cos y k — 'g 1 sin y k + g k g kT (1 - cos y k ) (4) 

and 

C Ktk = U cos y k + g k sin y k + g k g kT (1 — cos y k ) = (C kNl ) T (5) 

The proof of Eqs. (4) and (5), which can be found (in different terminology) in 
Ref. 26, is reproduced here for convenience. Figure 6 portrays two bodies and 
in some relative orientation for which y k ■=£■ 0. The unit vectors shown on the 
sketch are fixed in the two bodies in such a way that at all times X oc 2 = g* = 
Pi X p 2 , and when y k = 0, then {b*} = {b v *} (as in Def. 7), and also a 1 = p, and 
oc 2 = p 2 . Thus, for some constant values of p 8 , q s , and r s for s = 1, 2, 3, the equations 

bg‘ = p a ai + q 8 <x 2 + r s g k 

and 

b* = p a Pi + q s P 2 + r«g* 

must apply. Substituting the relationships 

Pi = cos y*oq + sin y k a 2 
p 2 = -sin yfcoq + cosy fc a 2 
into the preceding equations, one finds that 

b, = (p s cos y k - q s sin y k ) oq + (p, sin y k + q s cos y k ) a 2 + r s g k 
or, since otj X g k = — a 2 and s 2 Xg* = a 1 , 

b* = b^‘ cos y k — b'J* X g k sin y k + b*‘ • g k g k (1 — cos 9) (6) 



Fig. 6. Simple rotation of 4 k relative to ^ 
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From the definition of C r >, the element is given by 


c k s * u ' = K-K' ( 7 ) 

permitting, in combination with Eq. (6), the determination of the elements 

C™* = cos y k + (g *) 2 (1 - cos y k) (8) 

C k J‘ = g\ sin yk + g\g\ (1 - cos y k ) (9) 

C\f ‘ = -g* sin yk + g\g\ (1 - cos yk) (10) 

C“‘ = -g$ sin y fc + g^giS (1 - cos y«c) (11) 

C\* k = COS y t + (g|) 2 (1 - COS yk) (12) 

C*"* = g^ sin yk + g*g* (1 - COS yk) (13) 

c™‘ = g\ sin yk + g\g\ (1 - cos ys) (14) 

C\1 k = -g* sin yk + g^gt (1 - cos yk) (15) 

C™' = cos yk + (g5) 2 (1 - cos yk) (16) 


Although the validity of Eq. (4) can be established by expansion and comparison 
with Eqs. (8)-(16), this step may for some readers be more apparent with the 
intervening representation of Eqs. (8)~(16) in indicial or Cartesian tensor notation. 
With the definitions of the Levi-Civita symbol e a tc (Def. 23) and the Kronecker 
delta symbol (Def. 22), one can record the elements of C kNk generically as 

= 8„u COS y k + €,«,.g* sin y fc + g*g* (1 - cos y k) (u, v = 1, 2, 3) (17) 

and then recognize the validity of Eq. (4). Finally, one might prefer to construct 
directly from Eqs. (6) and (7) the direction cosine dyadic 


C k!il = U cos yk — U X g* sin yk + U • g*g* (1 — cos yk) (18) 

such that 

= (19) 

Any of these approaches can be used to verify Eq. (4), and perhaps to provide 
useful information or perspective as well. Equation (5) is, of course, available 
directly by transposition of Eq. (4), using 


C Nkk = (C w f 


(20) 


Equations (4) and (5) provide the direction cosine matrices between contiguous 
bodies only, and the equations of motion are more conveniently expressed in terms 
of the more general matrices C r, '(r, /e£ 3) provided in Def. 29. A direction cosine 
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matrix relating the orientations of two bodies connected by a chain of intermediate 
bodies can be written as the ordered product of the direction cosine matrices relat- 
ing contiguous pairs of bodies in the chain. This product can be written sym- 
bolically as 

C r > = n^C M,i (21) 

where the symbol 11%" is understood to imply the following algorithm: 

(1) Define p = N r j and construct C rp from Eq. (4) if r > p and from Eq. (5) 
if r < p. 

(2) Define q = N pj and construct C vq from Eq. (4) if p > q and from Eq. (5) 
if p < q. 

(3) Proceed until an integer u emerges such that j = N uj , finally constructing 
C ui from Eq. (4) if u > / and from Eq. (5) if u < 

(4) Multiply the matrices obtained in the sequence 


C r > = C rp C m ■ ■ • C ui 


( 22 ) 


D. System Specification 

Once the appropriate equations from the three preceding sections are embodied 
in a computer program, the user need not deal directly with all of the symbols and 
concepts introduced in the formulation of these equations. He can provide the 
computer with the very limited body of input required to specify his system, and 
then concentrate on interpreting the numerical integration output. Although the 
specific problems of programming the digital computer to process the input are 
deferred to a later report, the definition of the required input is to be established 
immediately. 


It is assumed that the engineer in his wisdom and experience has devised for 
his physical system a mathematical model consisting of a system of n + 1 bodies 
interconnected by n hinges (such as that illustrated in Fig. 4). He has labeled the 
bodies as indicated by Def. 4 (Section IIA), and fixed in his mind a nominal con- 
figuration (for which y fc = 0 for k e & and b* = b’„ for j,ke£B and a = 1, 2, 3). He 

must then provide the computer with the following information: 

Computer input (required) 

(1) The integer n (see Def. 1). 

(2) The n network elements N h (k e £?) (see Def. 26). 

(3) For kz(B and je& 3 fc , the 3 by 1 matrices p ki (see Def. 13). 

(4) The n 3 by 1 matrices g* for ke & (see Def. 14). 

(5) The n + 1 masses m* for (see Def. 17). 

(6) The n + 1 3 by 3 inertia matrices I k for ke£B (see Def. 18). 
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(7) For k e £B, the 2 (n + 1) 3 by 1 matrices F k and T k , either as explicit functions 
of system variables (possibly zero) or in the form of differential equations 
characterizing control laws (see Defs. 19 and 20). 

(8) For keFP, the n scalar functions r k , either as explicit functions of system 
variables (possibly zero) or as differential equations representing control laws 
(see Def. 21). 

(9) Initial values for the 2ra scalars y* and for keFP (see Def. 15), the nine ele- 
ments of the direction cosine matrix C (see Def. 9) and the three elements 
of the angular velocity matrix <o° (see Def. 10). 

(10) Initial values of X and X, to be included if F = QflX (Eq. 2) is included in 
the simulation program. More often, Eq. (2) will be solvable for vehicle mass 
center motion X(t) without reference to the rotational equations (Eq. 1), 
and X(t) will be input to the rotational dynamics program as an explicit 
function of time. Very often, Eq. (1) will not involve X(t) or its time deriva- 
tives, and then, of course, no input regarding X is required. 

Computer input (optional) 

(1) A set of attitude variables (attitude angles, Euler parameters, or other kine- 
matical quantities) to characterize the inertial attitude of / 0 ; if this option is 
exercised, an expression must be provided in terms of these variables for the 
direction cosine matrix C introduced in Def. 9, and initial conditions on 
these variables and their derivatives (consistent with constraints) might be 
adopted rather than initial values of C and o>°. Equation (3) can then be 
replaced by a kinematical equation appropriate for the variables selected. 

(2) Any kinematically prescribed variables, such as y,- (f) for some j e Q>. As noted 
in Section IIB, it is possible to modify Eq. (1) so as to accommodate the sub- 
stitution of an explicit function of time for a variable, such as y, ( t ), which 
in Eq. (1) is treated as an unknown. If this option is elected, the computer 
program input must include the specification of variables to be prescribed 
and the appropriate time functions. 

No additional input is required. The computer can be programmed to construct 
numerically all intermediate concepts (such as the neighbor sets £3* and the remain- 
ing network elements N hr for k,re£P, the path elements e sk for ke£B,seFP, the 
barycentric position matrices D kh and D kj for k,je£B, the augmented body inertia 
matrices $**, and the inertia-like matrices $ k ’ for k, / e£8). The computer can then 
evaluate numerically the functions displayed in Eq. (1) and, with initial values 
prescribed, accomplish the numerical integrations required for Eq. (1) and any 
necessary kinematical equations (such as Eq. 3). This process is illustrated in gen- 
eral terms by example in the next section. Explicit scalar equations of motion are 
obtained from Eq. (1) for a simple three-body system in Appendix B. 

E. Sample Problem Formulation 

Consider once again the eleven-body system portrayed in Fig. 4. Imagine that 
equations from the preceding sections have been programmed for digital computa- 
tions. In this section, an outline will be provided for the specific input required by 
this system, and the internal operations of the computer will be traced conceptually 
to the point of integrating equations of motion and printing out results. 
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According to the input list in the preceding section, the following system param- 
eters would suffice for a simulation of the vehicle model illustrated in Fig. 4: 


Required input 

(1) n = 10 

(2) Ni = 0 N 6 = 5 

N, = l N 7 = 6 

N 3 = 2 N s = 7 

N 4 = 3 N 9 = 8 

N 5 = 0 = 6 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


( 10 ) 


Numerical values for the 3 by 1 matrices p 01 , p 05 , p 10 , p 12 , p 2i , p 23 , p 32 , p 3 \ 
p* 3 , p 50 , p 56 , P 65 , P 67 , p 6 ’ 10 , p 76 , p 78 , p 87 , p 89 , p 98 , p 19 ’ 6 . In this case, p 21 = p 23 = 
p32 — p34 = ^43 = a nd other values are nonzero, but are not recorded here. 


Numerical values for the ten 3 by 1 matrices g k for k = 1, 2, • • • ,10. Here 
we adopt g 3 = g 5 = g 6 = g 7 = g 8 = [0 —1 0] T , g 10 = [ — 1 0 0] r , g 9 = 
[ - V272 VST* 0] T , and g 1 = g 2 = g 4 = [0 0 l] r . 

Numerical values for the eleven scalars m k for k = 0, 1, 2, • • ■ , 10, not fully 
recorded here. In this example, the masses m 2 and m :i are taken as zero in 
order to illustrate the possibility of simulating 3 deg of rotational freedom 
between <f 4 and 


Numerical values for the eleven 3 by 3 matrices P for k = 0, 1, 2, • • • , 10, 
not fully recorded here. In this example, I 2 and I 3 are taken as zero, and all 
matrices except P are assumed diagonal. 

Eleven values for each of the 3 by 1 matrices F k and T k , for k = 0, 1, • • • , 10. 
Although these are often complex functions representing environmental 
interactions and controller influences, in this example, the system is assumed 
to be free of all external forces, so that F k = T* = 0 for all ke$>. 

Ten scalars n, • • • , n 0 , which for this example are given by r« == —k a y a for 
a = 5, 6, 7, 8, 9, 10 and rp = 0 for /3 = 2, 3, 4. The scalar t, is unspecified here 
because will be prescribed as an explicit function of time. 

Initial values for y„ • • • , y 10 , y u ■ ■ ■ , y 10 , C and o>°, which for this example 
are given by y 2 — ys — y* ~ 0, y 5 = y$ = y? = y& — y$ = yio = tt /8 rad, y 4 = 
10tt rad/s, y 2 = 73 — y 5 ^ ys = yT = ys = y 9 = yio = 0, C = U, and o>° = 0. 
The function y x (t) is prescribed, as noted below. 

Initial values of X and X; in this example, these are irrelevant because they 
do not appear in Eq, (1), and it is assumed that mass center trajectory data 
can be obtained more efficiently from a special-purpose program for the inte- 
gration of Eq. (2). 


Optional input 

(1) Inertial attitude variables for here assumed to be the set of body three-axis 
angles 0 U 9 2 , 0 3 so defined that b 9 = i a when 0i = 62 = 0 S for a — 1, 2, 3, and 
in general, the attitude of {b°} can be obtained by rotating {b°} from a state 
of alignment with {i} first through an angle 0 i about ii = b®, then through 
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an angle 0 2 about the newly displaced b®, and finally, through an angle 9 3 
about b®. The resulting direction cosine matrix is given by (Ref. 26, p. 61, 
after transposition) 



C2C3 

c 2 s 3 

s 2 

c = 


SiS 2 S 3 “I - C3C1 

CiS 2 S 3 €3$! 


^Ci<S 2 C3 “1"' S 3 S 1 

S1C2 

CiC 2 


(23) 


where c a = cos 9 a and s a = sin Q a . 

Although one could simply input Eq. (23) into a computer program de- 
signed to integrate the kinematical equation C — — <5 °C (Eq. 3), this option 
has little merit when contrasted with the alternative of replacing Eq. (3) in 
the program with 

6 = Pa,® (24) 

where 0 = [0, 0 2 0 3 ] r and (from Ref. 26, p. 147, with P = M~ lT ) 


P 


A J_ 

C2 


c 2 

0 

0 


$ 1$2 


CiC 2 




Ci«9 2 


Ct 


Equation (23) is then no longer required input in this example since it does 
not appear in Eq. (1). 


Note that P is undefined for c 2 = 0, that is, for 0 2 = n/2, 3t t/2, etc. This 
singularity is characteristic of all three-variable systems for attitude descrip- 
tion. The choice of the particular set of attitude angles adopted for this 
example is based on the conviction that, for the given parameters and initial 
conditions, the variable d 2 will remain away from those values for which P 
is singular. 

(2) Kinematically prescribed variables, in this case being y, = sin St rad, with 
S given a numerical value not recorded here. 


This completes the input required by the computer. Since, ultimately, the differ- 
ential equations, (1) and (24), must be integrated, the computer must generate 
numerical values for all parameters appearing in these equations. Before this 
process is outlined, the equations to be integrated will be rewritten as a single 
system of first-order differential equations, in a form most suitable for computer 
operations. 

Define the 24 by 1 matrix x = [0, 0 2 63 a,® a,® a,® y 2 y 2 • • • yio yioF and recast the 
required equations as the state equation 


Vx = W (25) 

Note that the kinematically prescribed variables y 3 and y t have been excluded 
from the state variable x. The second-order equations in Eq. (1) have been cast 
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as first-order equations in Eq. (25) by the expedience of including y 2 — x s , y 3 — 

x 10 , * • • , yio = x 2i in the state variable x. This step necessitates the inclusion in 
Eq. (25) of identities such as x 8 = x 7 (or y 2 = y 2 ). For the particular arrangement of 
variables indicated for x, the matrices Vi and W in Eq. (25) are, for equivalence 
with Eqs. (1) and (24), given by 


Vi- 


and 


- 

1 1 

0 

1 

0 I 

0 

1 0 

u 

1 0 1 

0 

! 

o 1 

0 

1 0 


1 1 

0 

1 

0 1 

0 

1 0 


1 1 

0 

1 


0 

1 

0 

I doo I 

0 

1 

#02 

0 

1 #03 


J I 

0 

1 


0 

I 

000 

| 000 I 

I 


0 

0 

0 

000 

' &20 

0 


#22 

0 

0 23 

000 

1 000 | 

0 


0 

1 

0 

000 

1 #30 1 

0 


a 32 

0 

#33 









1 j 

1 oooj 




• 


000 

0 


0 

0 

0 

Looo 

1 # 10 , ol 

0 


010,2 

0 

010,3 







Po>° 


W = 


0 
0 
0_ 

0 
0 

.±L 


2 - Oom 


72 

g 2T 2 £2kC 2k A k + t 2 — 0 21 y! 
ktff 

73 

g 3T 2 E 3k C 3k A k + t 3 — 0 31 yi 
ta 9 


0 

0 

0 


0„, io 


0 0 

0 0 2 |, 10 

0 0 

0 0 3 , 10 


0 0i 


y 2 

72 

73 
73 


yio 

.7io. 


(26) 


(27) 


7i° 

L g lor 2 e 1#f *C w -*A‘ + T 10 -a 1 o tl y 1 J 

Uff 

Of course, the rows of V (and of W) could be shifted around by adopting a different 
ordering of the elements of x. The indicated choice is governed by considerations 
that become evident in Section III. 


The only task remaining in preparation for integration is the identification of 
procedures for evaluating the elements of V and W from the input. Functions such 
as P, t 2 , t 10 , g 2 , g 10 are given explicitly as input, and need no further consideration. 
Direction cosine matrices appearing explicitly in W are obtained from Eq. (21) or 
the algorithm following that equation, and to apply that algorithm, the network 
elements N fcr (k, re 13) are required. These quantities are provided as input for 
r = 0, and are generated by the computer from knowledge of the branch sets 
£3*i (ke£B, and/e£8k), which in turn are obtained from the neighbor sets £3* (ke£3), 
which in turn are available from the input quantities Nic(=N k0 ). 


24 


JPL TECHNICAL REPORT 32-1565 



Specifically, the neighbor sets £3* may be obtained from the input network ele- 
ments N k by the following interpretations of Def. 5, written here both verbally and 
in terms of set theory symbolism: £8 0 is the set of all integers / such that N,- = 0 
(or £B 0 = {/|Nj = 0}), and for ke&, is N k plus the set of all integers f such that 
Nj = k (or — ( N k + {j\Nj = k}} for ked>). 

By scanning the input values of N x through N 10 for this sample problem, the 
computer can readily find the following neighbor sets: 


mo = {i,5] 
m 1 = {o,2} 
m 2 = { 1,3} 

m 3 = {2,4) 
m, = { 3 } 
£8 5 = {0,6} 


m& — {5, 7, 10} 

m, = { 6 , 8 } 

m& = { 7 , 9 } 

m* = {8} 
m w = {6} 


The validity of these results can be established by inspection of Fig. 4, with Def. 5 
in mind. 

As a second major step, the elements of the 2 n branch sets 03 ki can be obtained 
by the following general algorithm: 

(1) Set k = 0, and from the previous calculation identify the elements of Oik as 
a u a 2 , • • • , a„ k , where n k is the number of elements in Q3 k - 

(2) Set <7 = 1, and record the calculated elements in 03a n \ then identify the ele- 
ments of {£B a<l ~ k) as b u b 2 , ■ ■ ■ , b„ a , where n aq is the number of elements 
in { m aq -k). 

(3) Set r = 1, and obtain the neighbor set Sib r from the calculations; then identify 
the elements of {<£8& r — a q } as c lt c 2 , • • • , c nb , where n& r is the number of 
elements in {(B ir — a q }. 

(4) Proceed as in (2) and (3), finding £B Cs for s = 1 and identifying the elements 

of (£3c„ — b r } as d x , d 2y ■ • ■ , d n „ , where n c is the number of elements in 
(£Sc s — b r }, proceeding in this manner, identifying elements e u ' ' ' , Cn df 
and f x , ■ ■ • , etc., until an empty set is encountered. (This sequence of 

operations cannot involve more than n steps.) 

(5) Upon obtaining an empty set, say the set ( £B lv — e u j for v = 1, change the 
most recently generated index (here v) from 1 to 2, and proceed with the 
indicated sequence of operations until an empty set is again obtained. Then 
set the most recent index (here v ) to 3 and proceed again to an empty set. 
Continue this process until v = n eu has been considered. 

(6) Change u from 1 to 2, and repeat the process from the point at which u was 
previously given value 1. Continue until u = n rf( has been considered. 

(7) Change t from 1 to 2 and repeat, continuing in this manner until the index r 
introduced in step (3) exhausts the n a? elements of the set {£Ba q — k), for 
q = k = 1. 
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(8) Populate the set 83ka q with the elements of the following sets: 

{ a «} 

{£Ba q ~ k) 

2" a o {$&, - a q ) 

r = l 

{$c„ -b r ) 

r = i s = 1 

• - • 2 n *« {£ 8 /, - «„} 

r=i « = i t? = i 

At this point, the branch set $ oa] has been obtained. 

(9) Return to step (2), change q from 1 to 2, and repeat the process to find $ 0 a 2 ’ 
repeating again until $ 00no is obtained. 

(10) Return to step (1), change k from 0 to 1, and repeat the entire process until 
the final branch set £B na is obtained. 

n n 

Although the indicated algorithm for finding the 2n branch sets may appear to 
be lengthy and arduous, it is exactly this process that a man accomplishes in a 
moment’s time when he examines a sketch such as Fig. 4 and records the branch 
sets “by inspection.” If you will close your eyes to the sketch in Fig. 4 and consider 
only the neighbor sets generated for that system from item (2) of the computer 
input, you will be able to obtain the following branch sets from the algorithm 
provided: 


$oi = {1,2, 3, 4} 

£3o 5 = {5,6,7, 8,9, 10} 
fiJi.o = {0,5,6,7,8,9,10} 

$> 2 = {2,3,4} 

$ 2 i = {1,0, 5, 6, 7, 8, 9, 10} 

$ 2 3 = {3,4} 

® 32 = {2, 1,0,5,6,7,8,9,10} 

$34 = {4} 

<$43 = {3,2, 1,0,5, 6, 7, 8,9, 10} 
<$50 = { 0 , 1 , 2 , 3 , 4 } 

$ 56 = {6,7,8,9,10} 

$65 ~ {5,0, 1,2, 3, 4} 

$ 67 = {7,8,9} 

$6,10 = {10} 

$ 76 = {6,5,0, l,2,3,4}i 
$78 = {8,9} 
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= {7, 6, 5,0, 1,2, 3, 4} 

£8sd = {9} 

£Soh = {8, 7, 6, 5, 0, 1, 2, 3, 4, 10} 

£B w .e = {6,5,0, 1,2, 3, 4, 7, 8, 9} 

These results are easily confirmed by examination of Fig. 4. By having the com- 
puter generate this information internally, we minimize user input and user error. 


The next step is construction of the network elements N kr for re£P and k efB ; 

A 

these elements were input for r — 0, since N k = N k0 . This is a simple matter since, 

from the definition, N kk = k and, for r=}=k, N kr — j for refB kj . Hence, for a given 
k and r, one merely searches 8b ks for s = 0, 1, • • • until one finds a value s = ; for 
which 0$ ks includes r, and then assigns this value j of s to N kr . In application to 
Fig. 4, the result is as shown in Table 1, which includes the input N k0 . 


Table 1. Network elements N kr for Fig. 4, k, r«£B 


\. r 

Ik \ 

0 

(input) 

j 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

0 

1 1 

1 

1 

1 

5 

5 

5 

5 

5 

5 

1 

0 

1 

2 

2 

2 

0 

0 

0 

0 

0 

0 

2 

1 

1 

2 

3 

3 

1 

1 

1 

1 

1 

1 

3 

2 

2 

2 

3 

4 

2 

2 

2 

2 

2 

2 

4 

3 

1 3 

3 

3 

4 

3 

3 

3 

3 

3 

3 

5 

0 

1 0 

0 

0 

0 

5 

6 

6 

6 

6 

6 

6 

5 

1 5 

5 

5 

5 

5 

6 

7 

7 

7 

10 

7 

6 

j 6 

6 

6 

6 

6 

6 

7 

8 

8 

6 

8 

7 

1 7 

7 

7 

7 

7 

7 

7 

8 

9 

7 

9 

8 

8 

8 

8 

8 

8 

8 

8 

8 

9 

8 

10 

6 

I 6 

6 

6 

6 

6 

6 

6 

6 

6 

6 


With the network elements available, the computer can, by means of the algo- 
rithm culminating in Eq. (22), construct any direction cosine matrix appearing in 
Eqs. (25)-(27) as a product of matrices whose numerical values are available in 
terms of input parameters g k and kinematical variables y k (keQi) from Eq. (4), 
Eq. (5), or Eqs. (8)-(16). For example, Eq. (22) provides 


C 27 = C 21 C ,0 C° 5 C 56 C G7 


05,10 — 05C06, it) 

040 — 043032021010 

and each of the matrices on the right sides is available from Eq. (4) or Eq. (5) (or its 
elements can be written directly from Eqs. 8-16). For example, from Eq. (4), 

C 21 = U cos y 2 — g 2 sin y z + g 2 g 2T (1 — cos y 2 ) 
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and from Eq. (5), 


C 05 = U cos y 5 + g 5 sin y s + g 5 g 5T (1 — cos y 5 ) 

In addition to the direction cosines, Eqs. (25)-(27) involve many other func- 
tions to be generated internally by the computer. Among these are the path ele- 
ments Est for s e£P and ke£B. These are readily determined from the branch sets 
$>r s (re£B,S£&) and the network elements A 7 ., (for se &) as follows: If k e£3, v?s , 
then £ S jfc = 1; otherwise, e.^ = 0. In application to the system of Fig. 4, this algorithm 
provides the results shown in Table 2. 


Table 2. Path elements e sk for Fig. 4, se£P, k e£B 


Is 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

2 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

5 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

6 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

7 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

8 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


All that remains for the computation of the terms in Eqs. (25)-(27) is the evalua- 
tion of the matrices A k (k e£3) shown in Eq. (27), and the evaluation of the symbols 
appearing in the upper left partition of V in Eq. (26). All of these symbols represent 
functions for which explicit expressions are provided in Defs. 38-41 (Section IIA). 
These definitions, in turn, involve symbols requiring reference to the preceding 
Defs. 33 and 35-37, and these require additional reference to Defs. 30-32; but 
finally, all quantities are available in terms of those provided in the input. Thus, 
the coefficient matrices V and W in Eq. (25) can be evaluated numerically for time 
zero, and the numerical integration process can begin. 

The discussion of this example is continued in Section HIE, where partially 
linearized equations are recorded and coordinate transformations and truncations 
are considered. 


III. Hybrid-Coordinate Equations of Motion 

A. Rationale 

In many applications, it can be anticipated that some of the kinematical vari- 
ables appearing in Eq. (1) will remain in some sense “small.” It may be, for example, 
that some of the angles y u • • • ,y n represent relative rotations of contiguous 
bodies connected by an analytical hinge, rather than a physical hinge (see Fig. 3, 
for example). Then, the characteristics of the hinge (stiffness, etc.) may be based 
upon an idealization which retains its validity only for small structural strains. 
Nonlinear terms in the angle of relative rotation then become meaningless and 
should be dismissed, with the understanding that solutions indicating large values 
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of such variables are also meaningless. Further argument must be advanced in 
order to justify the additional assumption that nonlinear terms in the time deriva- 
tives of these small angles are also negligible; this step would be justified, for 
example, if it could be anticipated that the angle would experience only low- 
frequency oscillations. 

Linearization of kinematical variables and their time derivatives might alterna- 
tively be justified in an entirely different way, relying upon mathematical theorems 
rather than physically based arguments. If an exact solution of Eq. (1) can be 
found, then this solution can be adopted as a nominal motion. By transforming to 
a new set of kinematical variables which describe the deviation from the nominal 
motion, and linearizing in these variables and their time derivatives, one can obtain 
a set of linear differential equations which, in many cases, quite rigorously establish 
the Liapunov stability properties of the nominal motion. The new variables are 
often called variational coordinates, and the transformed equations are then known 
as the linearized variational equations. 

In many spacecraft applications, an exact solution of Eq. (1) will not be avail- 
able, but still there will exist a desired motion which can informally take the place 
of the nominal solution in the preceding paragraph. The resulting linear equations 
may have considerable engineering significance, even though they are no longer 
rigorously indicative of motion stability properties. 

Linearization based on mathematical arguments is at best formally indicative of 
local stability properties of solutions; thus, all variational coordinates and their 
time derivatives are taken to be arbitrarily small. Since the sum of a finite number 
of arbitrarily small quantities is still arbitrarily small, a formal interpretation of 
linearized variational equations representing the deviation from an exact solution 
of Eq. (1) does not permit large relative motions of even physically separated 
bodies of the system. 

In engineering applications, however, linearization of a variable is generally con- 
sidered to be an acceptable practice as long as the linear term in the variable is 
“substantially larger” than additive terms of higher degree. In this sense, lineariza- 
tion is a process whose range of validity is somewhat ill-defined. For example, if 
we anticipate that y 2 will have a solution approximating y 2 = 0.2 cos 2 1 rad, then it 
is quite reasonable to replace sin y 2 = y 2 — y|/3! + y|/5! — • • • by y 2 (since 
sin 0.2 = 0.199 = 0.2) and to replace cos y 2 = 1 — y 2 /2! + y*/4! — • • by 1 (since 

cos 0.2 = 0.98 s 1); but one might hesitate to replace y 2 — y* = —0.8 cos 2 1 
— 0.16 sin 2 2t by y 2 = —0.8 cos 2t. 

Evidently subjective judgments are involved in the engineering interpretation 
of linearized equations. Within this framework, however, one can interpret linear 
differential equations obtained from Eq. (1) as descriptive of large gross deforma- 
tions of a vehicle experiencing small strains as manifested by small angles of rela- 
tive motion of contiguous pairs of bodies. If, for example, for the eleven-body 
system illustrated in Fig. 4, it could be assumed that angles y 5 , y 6 , y-, ys, yo, and y ]0 
remain “small,” then equations of motion linearized in these variables would still 
permit the dynamical description of “large” relative rotations of physically sepa- 
rated bodies such as <f 9 and / 0 , which might experience rotations on the order of 
magnitude of 1 rad (0.2 rad for each of five hinge angles), within the reasonable 
limits of the equations. If there were fifty bodies in the chain rather than five as in 
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the example cited, large gross deformations could be accommodated even if every 
relative rotation angle remained well below 0.1 rad. 

When any justification can be found for the linearization of any of the kine- 
matical variables appearing in Eq. (1), this step should be taken, since the mathe- 
matical and computational benefits of linearization are quite substantial. The 
primary immediate benefit is the elimination of variable terms from the coeffi- 
cient matrix on the left side of Eq. (1) (see Defs. 38-40). If «° and all angles 
yi, • • ■ , y n and their time derivatives can be assumed to remain arbitrarily small, 
then the indicated coefficient matrix is formally constant, and much computer 
simulation time can be gained by computing its inverse only once, in advance of 
numerical integration. When the angles are expected to be small in an engineering 
sense, less rigorous arguments can be marshalled to gain the same computational 
advantage. 

Potentially more significant advantages can be gained when coordinate trans- 
formations can be found which permit the substitution for the linearized variables 
of new coordinates, some of which may with impunity be deleted entirely from the 
dynamical description. In order for this coordinate truncation to be justifiable, the 
new coordinates must be distributed or modal coordinates which are to some 
degree uncoupled, and some of these coordinates must be demonstrably incon- 
sequential to the dynamic response. The subject of coordinate transformation and 
truncation has been explored extensively (see Refs. 21 and 25, for example), and 
previous results can be applied directly once some or all of the scalar equations 
implied by Eq. (1) have been written in linearized form. 

If all of the kinematical variables in Eq. (1) are linearized, then the transformed 
coordinates are vehicle normal mode coordinates, while if only a subset of the 
kinematical variables in Eq. (1) are linearized, then a hybrid-coordinate formula- 
tion results. 

B. Partial Linearization of Discrete-Coordinate Equations 

The immediate objective is to isolate a group of unknown kinematical variables 
appearing in the rotational equations of motion (Eq. 1), and to linearize these equa- 
tions in those variables and their time derivatives. 


In the special case in which all angles of relative motion between contiguous 
bodies remain small, it becomes advantageous to introduce a chain of three imagi- 
nary, massless bodies connected at one end to one of the rigid bodies of the system. 
This rigid body is then labeled and the chain of imaginary bodies is labeled se- 
quentially d 2 , and <f„. All other rigid bodies of the system are labeled according 
to previously established conventions, which apply now to both real and imaginary 
bodies. Unit vectors g 3 , g 2 , and g 1 are fixed along hinge axes in the chain of imagi- 
nary bodies in such a way that in the nominal configuration they are mutually 
orthogonal and g 1 X g 2 = g 3 . The hinge torques x 1 , x 2 , and x 3 are all set equal to 
zero. With these interpretations, Eq. (I) continues to apply without change, with 
the understanding that the “n + 1 rigid bodies” to which it applies now include 
the three imaginary bodies. In order to minimize computations, it is convenient to 
specify these imaginary bodies as having their mass centers all coincident with the 
mass center c 3 of d 3 , .and to let the hinge axes parallel to g 3 , g 2 , and g 1 all pass 
through c 3 . Figure 7 provides an illustration of a system for which imaginary bodies 
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have been introduced. The new (imaginary) reference body 4 0 now serves as a 
floating reference frame which accommodates the gross motion or mean motion 
of the collection of bodies in a sense to be established by the coordinate trans- 
formations in the sections to follow. The coordinate truncation which follows the 
transformation eliminates the coordinate redundancy introduced with the creation 
of three imaginary bodies. 

In the more general case, not all of the angles of relative rotation can be assumed 
(with their time derivatives) to remain small. Then, / 0 is selected as some rigid 
body of the system which is attached to some other body by a hinge where large- 
angle rotations are not precluded. 


In every case, it will be possible to identify some number of angles (say, v) for 
which linearization is in some sense justified, and to collect the indices of these 
angles in the set (A- If all angles of the system defy linearization, then v = 0 and 
(A is empty. There is then no recourse but to integrate Eq. (1) directly (numeri- 
cally, of course). If all of the angles of the system permit linearization (so that 
“the system” has been augmented to include imaginary bodies), then v = n and 
(A = {1, 2, ,n) = 0‘s More generally, we must expect 0 £ v ^ n. 
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In any case, Eqs. (1) can be rewritten in the form 


and 


a w o>° + 2 a ok : Y'fc = 2 C ok A k 

ktff keiB 


(28a) 


a»o«° + 2 — g lT 2 encC lk A k + Ti (ie&) (28b) 

It should be noted that, for a wide class of systems, it can be recognized in advance 
that eik = 0 for i e ffl and keff — ff\. This is the case whenever none of the hinge 
points whose indices comprise lies on a path from to any hinge point whose 
index is not in ff\. In physical terms, this means that the small-angle rotations are 
confined to terminal appendages, each of which is attached to only one body whose 
index lies outside of £ 7\ . In programming, it may be convenient to have the option 
to declare an appendage terminal and reduce the sum in Eq. (28) from the range 
kef? to the range k e ff[. 

If in the example depicted as Fig. 4, the angles y 5 , y 6 , y 7 , ys, y», and y 10 remain 
small enough to justify linearization, then v = 6 and (since the flexible substructure 
is a terminal appendage), whenever i = 5, • • • , 10, the summations over the set 
ke£P in Eq. (28) could be replaced by sums over £$. 

If, on the other hand, it became necessary to permit nonlinear terms in y 10 (still 
keeping y 5 , • • • , y 0 small), then one would have v = 5, and nonzero summations 
over the set Q — Q\. (The path elements e c .io and e 5 .i 0 would be nonzero, since /»,. 
and /«, lie on the path between / and 

The next step is the linearization of Eq. (28) in the v variables y k and their time 
derivatives for keffl. To this end, we expand each symbol in Eq. (28) which may 
involve such variables into three parts, the first being free of these variables (indi- 
cated by overbar), the second being linear in the variables (indicated by overcaret), 
and the third containing terms above the first degree in the variables (indicated by 
three dots). In particular, we write (for any i, k) 


C ik = c ik + C ik + • • • 

(29a) 

Ta — T« + To + ■ ■ ■ 

(29b) 

- , A 

dik — Oi k + Oi k + ' ’ ' 

(29c) 

— f\ 


A k = A k + A k + ■ ■ ■ 

(29d) 


and then determine explicit expressions for the new barred and careted symbols 
from the definitions of a ik (Defs. 38-40), A k (Def. 41), and the expansion for C ih 
(Eq. 21). In terms of the symbols introduced in Eq. (29), the linearized form of 
Eq. (28) becomes 

(a 0 o + a 00 ) <1° + 2 ( fl o k + Ook) y* + 2 &ok'h = 2 [(C ok + C ok ) A k + C ok A k ] 

— !A kefft kt(B 

(30a) 
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(i *&) 


(30b) 


(a i0 + &io) + 2 ( a ik + a-ik) yk + 2 ^iky* — 
keea-cn kem 

g iT 2 *ik [C*A* + ( C ik + C ik ) A k ] + u + 

kcff 

If the flexible substructure is one or more terminal appendages, then in Eq. (30b), 
whenever i e <71, the sum over the range ke£P can be reduced to the sum over k e <7\, 
since e ik = 0 for k e & — <71. 


It is convenient both for immediate computations and for later coordinate trans- 
formations to have Eqs. (30) written in matrix form. To this end, we can define 
column matrices y 2 , y*, • ■ • , y v+1 (with N an odd number) to consist of angles 
with sequentially numbered indices within the appendage set <7i, so that all angles 
y j for j e <7l appear in some matrix y k for some even integer k. Similarly, we define 
y 1 , y 3 , y 5 , ' ' ' , y" (for N odd) to accommodate those angles of the system not 
assumed small, so that all angles y,- for j e — (A appear in some matrix y k for k 
odd. As a general convention in what follows, an even index on a matrix estab- 
lishes an identification with a flexible substructure having small angles of rela- 
tive rotation, and an odd number is identified with a set of unrestricted angles. 
Equations (30a, b) then take the matrix form shown in Eq. (30c). 


_+ &°° | a 01 + a 01 | a 02 I a 03 + &° 


a 10 

4- S 10 

1 a 11 

4- a 11 1 a 12 | 

a 13 4- a 13 | • 

• 1 

a iv + 

T "a l ’ N+l 

a 20 

4- & 20 

! a 21 

+ a 21 1 a 22 , 

a 23 + cL 23 J ■ 

a 33 4- cL 33 | • 

1 

i 

a 2 * + a 2 ” 

| a 2 ’" +1 

a 30 

+ fc° 

| a 31 

+ a 31 1 a 32 j 

1 h 

: i : I 

a 3X + S 3V 




1 

1 

J 

• 1 ‘ 

1 

. i 


1 • 

1 ; 


a oti + & ON | a 0 ' w+l 


^jF' + a ’ 10 I a N1 + b N1 1 a K2 | a * 3 + & N3 ' • • • ' a NX + A™ a y - N + 1 

qW+i,o a N+1, ° I W* 1 ' 1 + a w+1 ’ 1 ^ "a‘ v +i.2 I' + + ^ w+1 - 3 ' '• ■ ” I a N * 1,N 4- Qw+i.x ■gw+i.w+i 


•*1 

y 


Ly 


y_ 

:W + i 


r 2 1 (c° k 

+ C° k ) A k + C ok A k ] 

- 2 a-oky'k 

ke£B 


keg 

R l (o>°, y 1 

~ • • -y ~ r 

>7 > 7 » 

,r+\t) 

R 2 (o>°, y 1 

. . . tf+l *1 . . . 

> Y > I > 


R 3 (o>°, y 1 

• • • ,y* +1 ,yS • • % 





m R N * 1 (< a °, 

y\ ■ ■ ■ , y w+ \ y 1 , • 



(30c) 


The many previously undefined symbols appearing in Eq. (30c) can be con- 
structed by comparison with Eq. (30b), with the set keQ defining given or pre- 
scribed variables. The purpose of writing the coefficient matrix on the left side in 
partitioned form is to facilitate the repeated inversion that is required for numerical 
integration. Any terminal appendages in the system can always be labeled so that 
the appendage bodies bear the highest indices in the system; this puts the corre- 
sponding angles in the matrix y x+1 . (If only internal flexible subsystems are present, 
there is nothing in y w+1 , and if all angles are small then y' v is empty.) Inspection of 
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the matrix ajm.w+i then reveals that it consists solely of scalar elements a ik for 
i, k e £$jv+i (the terminal appendage set), and it can be shown that these quantities 
are all constants. The right sides of Eq. (30c) will depend explicitly on t when 
either torques or kinematical variables are prescribed functions of time. 


Both for computational reasons and for purposes of coordinate transforma- 
tion, it is necessary to replace Eq. (30c) by a first-order matrix differential equa- 
tion. This can be accomplished in a variety of ways, and the optimum arrangement 
for computations is often inappropriate for the transformations to be discussed 
in sections following. Guided primarily by the latter consideration, we define 
(for j, k = 1, • • • , N + 1) 


r' 


A 



A 0 „ = a 00 + a 01 


K, = [0i a°i + d»i] (/ odd); A, 0 = AJ, 

A oj = [0j ~a oi ] (/ even); A j0 = A T oj 

A ^ = ["Y|^-+H (/ ’ t0dd) 

A ik = — ] (/ odd, k even); A kj = A T jk 

a r 8jkU j oi 
Aff = L""ol^ T j (f ’ fceven) 


and rewrite Eq. (30c) as 


> 

o 

o 

tb- 

O 

> 

c 

tc 

■ ' ' 1 Ao.W+l 

[_ 


• 0 

(O 


2 [(C<* + C° k ) A* + C<*&] - 2 Mfe 

kiS _ 

Aio | A ai | A 12 

1 1 

' ' ' i A liW+ i 

I 


f 1 


y 1 

R 1 

A 20 [ A 21 J A 22 

. . . 1 a 

^2, IV + 1 

1 _ _ 


f 2 

= 

• 2 

y 

R 2 

1 1 — 

Ajf + 1,0 j Atf + 1 .! | A.V + 1,2 l 1 Ay + 1, iV + 1 


pW + 1 


~ R N+1 


(30d) 


Whether we choose to work with Eq. (30d), or its equivalent, Eq. (30c), or the 
more explicit equations) in i_Eqs. (30a, b), any progress requires explicit expressions 
for the symbols C a , C ik , A k , A k , a ik , and aa for all permissible i and k. (We also 
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need T; and r*, but these can be written only for specific systems.) The required 
expressions can be obtained by combining Eq. (29) with the definitions of and 
A k in Section II A and the expression for C' k in Section IIC. 

Specifically, from Eq. (21) we have an expression which provides C ri (for r, f e£8) 
in terms of direction cosine matrices relating contiguous bodies, and they, in turn, 
are always available in expanded form from Eq. (4) or Eq. (5). The basic ingre- 
dients of C ri and & ri are therefore available from Eqs. (4) and (5) as 


C k - w * = c *- Nk 

for k e & — (ft 

(31a) 

C*'»* = u 

lor ke (ft 

(31b) 

A 

C k - N * = 0 

for k e & — (ft 

(31c) 

C k ’»' =-y k g k 

f or k e (ft 

(31d) 


When C r ’ in the symbolic form of Eq. (21) is expanded according to the algo- 
rithm culminating in Eq. (22), this latter expression then becomes 

_ a _ a _ A 

C ri = C rp C M • • • C ui gg ( C rp + C rp ) (C p? + C m ) ■ ■ ■ ( C ui 4- C ui ) 

£* C rp C M • • • C ui + C rp C M ■ ■ ■ C ui 2 (e kj - e kr ) g*y fc (32) 

ke (A 


so that, with indices generated according to the algorithm following Eq. (21), 


C r > = C rp C M • • • C u > (33a) 

b'i = ■ ■ ■ C“> 2 («.i - e,r)g'y, (33b) 

se!A 

In substituting Eqs. (33) into Eq. (30), simplifications can often be realized by 
recognizing Eqs. (31). If, for example, the flexible substructure consists of a single 
terminal appendage, then, in Eq. (30), the direction cosine terms for r,fe (ft sim- 
plify to 

C" - u C'l = 2 («., - «,r) g*y, (34) 

t(SA 


Equation (30) then becomes 

(flo„ + a on ) 6° + 2 a ok y k + 2 (aok + a ok ) y k = 2 [{C ok + C ok ) A k + C ok X k ] 

ktSA ke&-Cft k® 

(35a) 


(a i0 + a i0 ) <5° + 2 a ifcV + 2 (dot + k) y* — 

kcra kt®- (A 

g iT 2%{^ + [f/+2 Usk - «.,) g'yj] A*} +u + Pi (ic&) (35b) 

ke&l se&l 


Whether we retain the general result m Eq. (30) or only the special equations in 
Eqs. (35), we require expressions for A k , A k , &i k and ~a ik . These symbols can be 
expressed in terms of more basic quantities by expanding Defs. 38-41 A in the man- 
ner indicated by Eqs. (29c) and (29d). In terms of the symbols C rj and C r ’ available 
from Eqs. (33), these definitions lead to the following, in which we have substituted 
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_ A — A 

T* = T k 4- T* 4- ■ ■ • and F k = F k + F k + - in the manner established by 

Eq. (29): 

a 00 = 2 C°W° - 07? 2 2 (C ok D ikT C ik D ki - C°’D> k D k ’ T ) C k0 (36a) 

Jte® 

For 

flofc = 2 E kr C 0r 4> rr C rk g k - cm 2 2 E kr (C 0r D’ rT C ir D Ti - C 0 ’D’ r D T ’ T ) C Tk g? 

rtf? rtf? it8S-r 

(36b) 

«*o = «ok (36c) 

For i,ke£P, 

a ik = g iT 2 eiretr&WC^ 
rtg> 

- VYlg iT 2 2 e ir e*, ( C ir D‘ rT C‘ r D r ‘ - C u D sr D rsT ) C rk g? (36d) 

se^?— r 

For fce£8, 

A k ~ T k + 2 D ki C ki F> - 2 frkyr (C^C 0 *' 4- 2 csrj«C ks g*C* r ) g T 

j*85 riff-W st&-Q ! 

- [C*V>C° fc + 2 $** [C*V 4- 2 e ik yjC ki gi] 

4Qi! 2 {D*'C‘<[CW4 2 er>yrC ,r g r C rj ] 2 D’ k 

1*35 — k re^-fi* 

4- (UD ikT C ik D k> - C ki D’ k D k ’ T ) 2 e ri y, [C^C^ 4- 2 e.rf.C^'g’C"] g'} 

re<£?— £3 j 

(36e) 

So,, = 2 (<W*C*° 4- C ak & k C k0 ) - 07? 2 2 [(C°*D'* T C jk D* # - C°’D’ k D k ’ T ) C k0 

k*35 k*35 i*(3-k 

4- (C ok D’ kT C ik D ki + C ok D> kT C’ k D ki - C°’D> k D kiT ) C*°] (37a) 


For ke&, 

t k = 2 e kr ( C 0r 4> rr C rk + C°'4>"C'*) g* 

rtff 

- 07? 2 2 % [(C 0r D’ rT C’ r D r > - C°’D’ r D r ’ T ) C rk g k 

r*S3 jeffi-r 

4_ (C 0r D> rT C iT D r > 4- C 0r D’ rT C ir D r > - C°’D’ r D riT ) C rk g*] (37b) 

(37c) 

For i,ke{P, 

&ik = g iT 2 e ir e kr (6 ir $ rr C rk 4- C ir ®"£'*) g k 

r*ff 

- 07?g ir 2 2 *ir<* 8 [(C‘ r D ,rr C sr D r ' - C u D ki D ikT )t rk ^ 

Tiff ttP-r, 

4- (& ir D arT C er D r ’ 4- C %r D ,rT t ,r D ra — fc u D ki D ikT ) C rk g k ] (37d) 
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For ke£B, 


fa = T k + 2 D ki ( C k ’F> + C k ’F>) - 2 ewcvv [£*V>0 + C ko 7S°t or 

it® ttff-m 

+ 2 e, r y.C k 'g‘C’ T + 2 e.ry» {£ k ’g‘C‘ r + C* 8 g‘C 8r )] g r 

se^ st£P—(A 

- * kk 2 e rk y, (C‘VC»' + 2 e, r y.C ks g‘C‘ r ) g r 

rt{A stS'-ffi 

- [& ko Z°C ok + C ko 3°C ok + 2 e ik yj (C k >g>O k + C k 'g>C> k ) 

it&-!A 

+ 2 «irf#C«g'C'»] *» [C*<V> + 2 e ik y,C k >g>] 

id ? l it&-!A 

- [C k W k + 2 e ik y,C k ’g^C’ k ] 4>“ [C^° + 2 e^fa’g' + 2 e^C^g’] 

jt9-(A it£P-!A \t(h 

+ <yn 2 j D k >C k > [CWC°> + 2 ErjyrC^g r C^] 2 D )fc 

idB-k * U.9-CA 

+ {[C'W oi + C'VC 0 ’ + 2 £riy r (C’ r g r C r > + O r g r C ri ) 

rtS’-fA 

+ 2 «r^rC^rC^] [C’WC 0 > + 2 ^YrC ir g r C r ’] 

refi3 rtff-CA 

+ [CW> + 2 s r/ y r C> r g r C Ti ] [C'W' + C'V'C 0 ' 

+ 2 e ri yr (& r i r C'> + C"g'C") + 2 6nyrC'f r C^]} D ,k j 

rt&~(A rdA 

+ cm 2 j (UD ikT C ik D ki - C k ’D’ k D k ’ T ) 2 «,/tr [C^VC 08 

fdB-k rtfP-rA 

+ 2 Cry.C^C 88 ] g 8 + ( UD’ kT C’ k D k ’ - C k ’D ik D k ’ T ) { 2 e ri y r [C^C 08 

xff-CA rtffi 

+ 2 £srysC kg g*C >r ] g r + 2 £r>yr[C' [(> 3 <, C‘ ,r + C*VC 08 
+ 2 £»ry, (C* 8 g 8 C 88 + C k 'g‘& r ) + 2 e ar y s C ks g g C’ r ] g 8 } j (37e) 

xff-(A >tffi ’ 

(The superscript 2 in Eq. 36e indicates the square of the preceding bracketed 
matrix.) 


Upon substituting Eqs. (31), (33), (36), and (37) into Eq. (30), one can obtain the 
most general form of the partially linearized set of differential equations stemming 
from Eq. (1). These equations (to be assembled in detail by computer program) 
are then in a form whose dependence upon the small angles y«, • • • , y (l and their 
derivatives is fully explicit; by examining the structure of an appropriate subset of 
these equations, we can determine those circumstances for which transformation 
to distributed coordinates and subsequent coordinate transformation is feasible. 
It should be noted that in the special case of Eq. (30) recorded as Eq. (35), corre- 
spondingly simplified versions of Eq. (36) and (37) can be used. For example, for 
the special case in which the flexible substructure consists of a single terminal 
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appendage, one can simplify Eqs. (36) and (37) by noting that if <*_is the_append- 
age hinge closest to <f 0 , then, for ke&{ and ie& — e ik = e ia and C' k = C*°. Since 
A k and A k appear in Eq. (35b) only for k e (A, this simplification is computationally 
quite significant. 


Before proceeding with the search for coordinate transformations for the small 
variables and their first time-derivatives, it should be noted that the partially 
linearized discrete-coordinate equations are themselves of some computational 
value. Depending upon the extent of the partial linearization, Eq. (30) (or for 
terminal appendages, Eq. 35) may be much more amenable to numerical integra- 
tion than the unrestricted counterpart, Eq. (1). If Eqs. (30) are to be integrated 
directly (without coordinate transformation), then the suggested imaginary bodies 
should not be introduced, even if all angles remain small. 


C. Linear, Constant-Coefficient Differential Equations for Coordinate 

Transformations 

Equation (30b) includes the set of v equations 
2 OifcV'fc + a»oo>° + 2 flifcV'k — g' T 2 eik [C 1,[ A k + C lk A k ] — u = 

ke(A kc3i-(A ktff 

a io w° 4- 2 dayic + g' T 2 eikC' k A k +7* {ie(A) (38) 

fcc£B— £3 ke& 

and these will provide the basis for the desired coordinate transformation. The 
hybrid-coordinate procedure entails the transformation of the set of variables defin- 
ing the substructure deformation (here yi and y; for ieffl) into a set of distributed 
coordinates for which coordinate truncation can be accomplished without jeopar- 
dizing the salient features of the dynamic response. Truncation is never a rigorous 
mathematical process, and its justification must be based on rather subjective judg- 
ments of the acceptability of certain engineering approximations. Great caution is 
therefore necessary in adopting a truncated set of coordinates, and all possible 
effort should be made to select the coordinate transformation for which maximum 
truncation of coordinates can be accomplished with minimum risk. In order to per- 
mit the systematic evaluation of the consequences of truncation, it is important 
that coupling among the scalar equations of motion in the transformed coordinates 
be minimized. (If coupling were entirely eliminated, one could solve each of the 
transformed scalar equations in turn, and obtain the true solution by superposi- 
tion; then, irrelevant coordinates could be discarded with no risk of adverse 
consequences.) 

In the quest for a coordinate transformation which minimizes the coupling 
among distributed coordinates, the theory of linear, constant-coefficient, homo- 
geneous, differential equations provides our only guidance. Although in fact, 
Eqs. (38) are not in this category, in application, it is often reasonable to approxi- 
mate slowly varying coefficients or coefficients with small magnitude variations as 
constants, and to ignore the indirect coupling of Eqs. (38) through kinematical 
variables external to the flexible substructure, focusing only on the homogeneous 
counterparts to Eqs. (38). These approximations will therefore be adopted in this 
section for the purpose of finding a coordinate transformation, after which the 
inhomogeneous equations (30) will be transformed and retained for simulation. 
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The homogeneous counterpart to Eq. (38) is 

2 a ik y k + £ i0 »° + 2 a ik y k - 2 e ifc (C ik A k + C‘*A*) - h = 0 (t e 

iter?) Iceffl-r?) heff 

(39) 

All careted symbols in Eq. (39) are linear in the small-angle variables of the 
substructure. From Eq. (36e), A* depends on to 0 , so that the coefficients of the small 
variables in Eq. (39) are not constant unless o>° = 0. This restriction also eliminates 
the second term in Eq. (39). Moreover, the dependence of these coefficients on 
direction cosines which involve all of the angles of rotation external to the append- 
age requires that these angles be constant if Eq. (39) is to provide a coordinate 
transformation based on constant-coefficient differential equation theory. This 
restriction removes the^third term from Eq. (39) and greatly simplifies Eqs. (36e) 
and (37e) for A* and A*. As the basis for coordinate transformation, we adopt 
instead of Eq. (39) the following approximation of the equation: 

2 fi*y* - g iT 2 «« (C ik 2 k + C ik A k ) - $ = 0 (tea) (40) 

kt® 

where a ik is a constant defined by_Oi* in some nominal state (see Eq. 36d), and 
where, from Eqs. (36e) and (37e), A* and A k are given by 

A k = T k + 2 D ki C k >F’ - c ko &c ok $ kk c ko a° + (yyi 2 D ki C k ’ (C’°S°C 0 ') 2 D’ k 

m m-k 

(41a) 

A k = T k + 2 n ki (C k ’F> + C k ’F>) - $“2 erkyrC k<f a>°C 0r g r 

i*£b rc& 

- [C k ^C ok + C*6* + 2 wi&g'C**] <P kk C k W 

- C k ^C ok ^ kk [C k0 <7,° + 2 tjkyfi* 1 g'] + 041 2 j D k ’t k > (0>°5°C O >) 2 D’ k 

idB-k' 

+ D ki C k i {[C'V’C 0 ' + C'°5°^ 0 ’ 4- 2 e r ,irC> r g r C r >} C>°5°C 0 ' 

rtffi 

+ Z: I0 Z°C 0 ’ [^'WC 01 ' + C io 7S°€! oi + 2 er,yr£ ,f g r C r# ]} D lk 

«£?) 

+ (UD’ kT ^’ k D k > - Z: k iD’ k D kiT ) [ 2 er/yJ&*&V tr g r ] } (41b) 

n £> i ’ 

It is understood without the introduction of more notation that each symbol in 
Eqs. (40) and (41) is constant, except for those depending on y, or its time deriva- 
tives, for / € (A- 


Equation (40) is typical of v scalar equations of motion, which can be written as 
a single matrix equation of the form 


I'y 4" Jy + Ky — 0 


( 42 ) 


where, if the v indices in (A range from « to /x, 


y = [y« * ‘ • y>>] 
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and where I is the v by v symmetric matrix whose elements are a,, for i,je <3\ (see 
Eq. 36 for an). As indicated in the footnote following Eq. (1), the variable coeffi- 
cient matrix with elements a ik is nonsingular, so that the constant matrix I is also 
nonsingular. 


The v _by v matrix / comprises the coefficients of y, in terms involving A k , and 
perhaps r i; for i,j and ke0>. More specifically, a typical element of ] may be 
written for t, j e 63 as 


hi - ~g 


iT 


2 

kt& 


= 

C' k — 


9^ 


Syy 9 y i 


(43) 


The elements of the v by v matrix K are similarly available as 

A 


Kit — —g' T 2 Cik 

kt& 


dC' k = = 9& 

I - — A k + C' k — 
L h t °y 


i _ 


9ft, 

hi 


(44) 


Since coordinate transformations which may usefully be applied to Eq. (42) are 
more readily calculated when the matrices J and K have certain properties of sym- 
metry (or skew symmetry), we might profitably examine the structure of these 
matrices. 


The hinge torques r, for ieff\ will in most cases be representative of structural 
connections, in which case they will usually represent either a linear elastic spring, 
for which 

£ = ~k iyi (45) 

or a viscoelastic connection, for which 

Ti = —kiyi — dy yi (46) 

In both cases, the contribution of the hinge torque to K is symmetric, and the con- 
tribution to } is zero for the elastic spring and symmetric for the viscoelastic con- 
nection. Only in the unusual case in which ftj is established by a control law that 
depends on appendage body relative rotations other than y< will the hinge torques 
contribute any but diagonal terms to the matrices / and K. 


In order to assemble the matrix equation (42), we require the partial derivatives 
dC' k /dyj, dX k /dyj, and dX k /dy ;, for t, je (A and From Eq. (33b), we have 


~ = C ip C pg ■ ■ ■ C uk (e ik - eji ) g’ (47) 

hi 


where again the algorithm following Eq. (21) is used to generate the unspecified 
indices. When the flexible substructure is an external appendage and a is the index 
of the hinge point closest to 4 0 , Eq. (47) becomes 


a c ik = = 

=r C ap c vq ■ ■ ■ 

hi 


C* k {e ik - ei i)W 


(48a) 
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When k e <71, both of these expressions reduce to 

dt ilc 


d yi 


— (ejk Bji) g’ 


(48b) 


Equation (41b) provides, for k e Q> and ; e <71, 

o T* k ~~ — p , _ 

— = — + 2 D k, C k> — 4> kk e jkC k0 Z? > C 0 ’ g’ - e ik C k >g’C’ k <£> kk C k(, Z 

°7i °y> 0 y, 


- C ke S°C ok ^*t ik C kl g i + 0>2 2 e,-, [0**0 (C ,^ g'C' <l C s VC 0, 

+ C‘^C os C* i g’C>’) D‘ k + ( UD* kT C ,k D k ‘ - C k, D ,k D k,T )C k0 s°C 0 > g>] 


(49) 


Also available from Eq. (41b) is the expression 

- = — 4 - v m*/ f?** - _ 
d yi 


* a r k ^ ~ /= 3F® ac* s = \ 

-=r-+2D‘' C l, r + F s ) 

9yy ^ \ Oyy dyj / 


rar fco — = 3C ok- i = — — dC ko 

- X>C° k + C k0 7j° — 4 > kk C ko m° - C ko T5°C ok $ kk — — fti° 


L 9y/ 


9yy J 


d Yi 


D k * — — [C a0 5°C 0 ’] 2 D‘* 

°y i 


+ D k, C k ’ 


era c so _ = = _ ec o, i= _ = 



rac so _ = = ac^T) 


4- C’^C 0 ’ 




C so S 0 C 0s 


D rt 


(50) 


For the special case in which there are no external forces and body /„ is nomi- 
nally at rest in inertial space so that T k = F k = a>° — 0, Eqs. (49) and (50) reduce to 


8& __ 3& 

9y>- 3y; 


as might be expected. If, in addition to these simplifications, we have linearly elastic 
spring torques at the joints, so that tj = — ki y* for i e <7i, then the equation to be 
examined for coordinate transformation (Eq. 42) is simply 


Iy + Ky=0 (51) 

where K is the diagonal v by v matrix with nonzero elements Ka = K t , and the 
elements of the symmetric v by v matrix I are given by Eq. (36d) 

A second special case of interest is one for which all hinge lines are parallel, so 
that for all i € TP, we can replace g< by the new unit vector g; in addition, the 
angular velocity vector o° has the same orientation, so that co° = Qg for some scalar 
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spin rate Q. This special case will be called the transverse deformation case for 
brevity; it should be understood that this term is applicable to three-dimensional 
systems such as that illustrated in Fig. 8. Substantial simplification of the preceding 
equations then follows from the vector identities (for i, j € &>) 


63° X g* = 0 

(52a) 

g» • 6)° X (any vector) = 0 

(52b) 

g‘ • g’ X (any vector) = 0 

(52c) 

In matrix terms, these identities can take various forms, depending on the vector 
basis chosen for matrix representation of the various vectors involved. Since we 
have chosen to employ the matrices defined for i e 9 s by 

£o° = {b 0 } 1 ^) 0 g‘ = 

(53a) 

with vector bases related for i,jeQS by 


{V} = C« (b’} 

(53b) 

we might expect to find Eqs. (52) reflected in our matrix equations in the form 

w°C ol g’ = 0 

(54a) 

g iT C i0 S°( ) = 0 

(54b) 

g iT C'>g’( )=o 

(54c) 


Because more than one change in vector basis might intervene, we might find 
instead of Eqs. (54b) and (54c) equivalent expressions requiring substitution of the 
identities (for i, j e 03) 


C *1 = C ik C ki (54d) 

and 

C ij = c ik C k ‘C*i (54e) 


In the transverse deformation case, the final sum in Eq. (49) for dfc/dyj vanishes 
by virtue of Eq. (54a), as does the third term on the right side. Moreover, most of 
the terms remaining in Eq. (49) vanish when it is substituted into the expression 
for Jn in Eq. (43), because of identities (54b)-(54e). All that remains in Eq. (43) is 


= fdT k 
Jii = ~g r 2 e ik C' k — 
ktff L °y> 


2 b k ‘C k ’^r 

«fs °y 


F*" 
Vi . 


0^ 

3 ii 


- g T 2 e ik C ik mi 2 euD k ’ (C k ’gC i0 ngC 0 ' + C^agC^gC’’) D ,k 

(55a) 
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If, in addition, any external forces (F 8 ) or torques (T*) on the bodies are inde- 
pendent of the relative rotation rates (y ,-), then the_first summation over ke£P 
disappears from J if . The interbody hinge torques (t 4 ) will contribute to Ja if 
“damping” is included in the joint, or the hinge rotation is subject to automatic 
control with rate feedback, but only in extraordinary cases will t ; contribute to 
Jn for r =)= /. For simplicity in what follows, we ignore any contribution of hinge 
torques to ] a. Moreover, we imagine that in the nominal steady state, all direction- 
cosine matrices equal the unit matrix, with no loss in generality. Then we have 


Jif = -2(tong T 2 «« 2 e /t &"ggb* (55b) 

keff sf(B~k 

or 

Jif = 2Q^ng T 2 2 e ik e js D k, gD‘ k g (55c) 

klff stf?-k 


From Eq. (55b), we see that ] t , = 0 if all D qr for q,r e(B are parallel. This restricted 
transverse deformation case will be explored further in what follows. Before pur- 
suing this course, we should note that in the more general transverse deformation 
case represented by Eq. (55c), the matrix / is skew symmetric, so that Jj, = — J ;i . 
To establish this result, note that the scalar must equal its transpose. Thus, 
Eq. (55c) provides 


Jn = -2 <Yhcig T 2 2 ti k e is D ,k gD k, g 

ktff stty-k 


Exchanging i and / in this expression, we find 

Jfi = -2tynng T 2 2 EikEisD sk gD k ’g 

kft? stff-k 


After exchanging the dummy indices s and k, we obtain 


Jn = -2Cty?ng J ' 2 2 ej,e ik D k, gD ,k g 

stff ktff-s 


(56) 


Comparison of this expression with Eq. (55c) establishes the claimed relationship, 
Jfi = — Jif, since in the course of each of the two indicated summations, both 
s and k assume all values in the set 


The next question of symmetry in the transverse deformation case concerns the 
matrix K, whose elements are given by Eq. (44). Because the externa] forces and 
torques represent a variety of influences (environmental interactions, control forces, 
etc.), their contribution to K cannot be described in comprehensive general terms, 
and they will be ignored in the discussion of the symmetry of K. The hinge torques 
Tj are also ignored. With these omissions, and the vector-dyadic identities 


g‘ • to 0 X < b kk • w° = g * ft 2 g X Q kk • g = 0 


(57) 
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which apply to the transverse deformation case, Eqs. (44), (41a), and (50) combine 
to provide 




( C*** - — 

Ki,= -g r 2 eiJ-r— 2 n 2 D*'C*'(C'°gC°') 2 D r * 

uff l °y> rtm -k 

= ~ r a c ks — - 

+ C'Kyyi 2 n 2 D lc * -z — (C 80 gC 0s ) 2 

Mffl-jt L °y> 

A A k 

— / r)C>° — — 3C 0S \ = = 

+ C ** ( IS 08 + C‘°g-r— ) C 30 gC 03 

\ oyj oyj ) 


— _ / a A so = — 3C 0s \T 1 

C*°gC°‘ r — gC°* + C°‘g — ll D«* | 


A 

3 Vi 


(58a) 


The structure of the matrix K is more apparent when in the nominal steady-state 
motion all direction-cosine matrices are the unit matrix. Then we have 

A 

( rtC ik ~ 

Kit = -Qm*g T 2 2 D* r ggD rt 

kiff l ®Vi «$-* 


+ 


2 

sefB—k 


~ rd£ k ‘ 

D ‘ L 3vi 


gg + 


A A ~K 

dC s0 _ _ dC os „ _ 3C®° 

gg + g^~ g +g‘ 


3yi 


3yy 


3y; 



From Eq. (33b), it follows that 

C*o + C 0 ’ = 2 (e P o - *p») gy P g + 2 (e v , ~ £ P o) gy P g = 0 

j*£S pe5<) 


so that Kn simplifies to 


/ dr'** ^ dC 1 * 8 'i 

K t , = —QWCl 2 g T 2 2 D fc * ggD* fc + 2 D k ‘ — ggD‘ k V 

l d Vi «®-k OVi ) 

A.. A. 


te9 


f r3C ifc ~ ~ 3C fcs ”1 ___ 1 

= -0t!n 2 g r 2 s 2 -r — 6* 8 + D** -r- — ggD 8 * > 
1*0 Uffl-*'L 8 Vi oy> J 1 


= -Q?m 2 g r 2 e« { 2 [(«* - »#i) gD* 8 + D*’g ( e „ - «,*)] ggD ek ) 

kef? se£R—k 

= -QW1 2 2 2 (e>. - e,*) g T D k, gggD sk (58b) 

kef? se(B—k 

since g T g — 0. Because e, 0 = 0, the range of s can be modified to obtain 


Ki, = -Qnn 2 2 2 £i k ej,g T D k *gggD** + Own 2 2 2 E ik e jk g T D ke gggD sk 

kef? sefP—k kef? s*£B—k 


(58 c) 
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By exchanging i and / in this expression, one finds K,-j and discovers that, since 
s and k both cover the range of & in the first double sum in Eq. (58c), we have 
K ti - Ki, = 0. 
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We may conclude that, at least for a vehicle of transverse deformation configu- 
ration with linearly elastic hinge torques and no external forces or torques, the 
vibration equation 7y + Jy + Ky = 0 found in Eq. (42) offers the following sym- 
metry properties: 7 and K are symmetric and / is skew- symmetric. Such a case was 
illustrated previously in Fig. 8. Moreover, if all vectors D ,r are parallel in the 
undeformed configuration for q,re£B, then (from Eq. 55b) the matrix ] is zero. 
We may call such configurations rectilinear transverse deformation cases. Figures 9 
and 10 illustrate examples of this kind, identified respectively as polar and equa- 
torial rectilinear transverse deformation cases. 

In the quest for special cases with demonstrable symmetry properties, we might 
next consider the configuration for which g* = g for all i e £P; g • ea° = 0, and in the 
nominal state, g • D rk = 0 and u° • D rfc = 0 for all r,keff>. Figure 11 illustrates this 
configuration, which is here referred to as the meridional deformation case. 

In order to examine the structure of the equations of motion in this case, we 
choose b°,b“,b® such that 


to 0 = nb? 

(59a) 

= b“ for ie<P 

(59b) 


D rfc = D rk b' for r, k e £ 8 


(59c) 


assuming moreover that in the nominal state, b( = b° for all r e <£?. With these sub- 
stitutions, the matrices ] and K which first appeared in Eq. (42) take the form of 
Eqs. (60a, b) (see Eqs. 43 and 44, and note that 

= {b k ) T b, 


b* = {b*} r 


0 

1 

0 


for any ke£3): 


and 


hi 



dk k 

d yi 


3 $j_ 

dyj 


Kij = - 6^2 e ik 
US' 




A 

dC‘ k 

dyj 


A k + 


sAn 

d Vi J 


dyj 


(60a) 


(60b) 
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Fig. 10. Example of "equatorial transverse deformation case” in nominal configuration 



i b <*5 i \ *2 

Fig. 11. Example of "meridional deformation case” in nominal configuration 


46 


JPL TECHNICAL REPORT 32-1565 



These equations simplify most substantially when, as a further restriction, it is 
assumed that in the nominal state, the principal axes of the augmented body 
barycentric inertia dyadic parallel b°,b°, and b° for all ke£B. (Such a case is 
illustrated in Fig. 11.) In this case, Eqs. (59), (60a), and (49) can be combined to 
provide 


A 

idT k ~ 3 F e ~ ~ 

Jii = ~2 tikb T A— + D k ,b 3 — E jk$™b 3 - EjfcO (b^bj - b^b^) 

ktff ' ( 9 Yj h, 

+ 0)7 2 nD k ,b 3 Ej, (bibt -f bj> 2 ) b 3 D ek 

sc£B—k 

+(M 2 e jr D rk D kr (Ub'b 3 - btf) b^A - 

rt£B-k ) °Y> 

= ~2 E ik bl (61) 

fce^p ( 0 y i 9 y » ) 0 y> 

This special case is comparable in its simplification of } to the transverse deforma- 
tion case (see Eq. 55). 

The matrix K in the meridional deformation case, with the indicated restrictions 
on the inertia dyadics, has elements given by Eqs. (59), (60b), (50), and (41a) in 
combination. The result is 


A 

r 3 F> 


Kj> — —2 e ikbl 

ktQ 


[T k + 2 DkrbsF ' - n *b t 9>*b t 
\ «ss 


A A A 

, „ r7r , „ , 3 T k „ „ T / 3F* , 3 C k ‘ = \ 

+ 0 2 &Dk r b 3 b l b i b 3 D rk ] + — h 2 F) ks b 3 ( — — 4- — F* 1 


. r 3& 0 


r r 3^°" 1 

L 3y; & * + bl 3y/ J 


A 

— 0C' 10 

— n 2 b^ k — — hi 

3y> 

5 36 


+ 04? 2 n*D w fc, + 

/ 3y> iL oyj 3y> J 

, r RC-t , r SC° J ~n , „ I 

bl L 3 y; &1 _ bl 3 yJ} aD,k [ 3 y/ 

[3* + 2 Dtrbafl + —■ + 2 D ks b 3 (■— + F* ) 
1 0 y> 3y> \ 0 y> 0 y> / 


= -2 

US} 


n 2 b, 


r RC 0 * „ , „ dC ko i 

bir* 1 ~3irJ 




r ~ rC° s 

+ On 2 n*D k ,D sk b-A \-^-b 2 + b 1 ^~bA\ 

s((H-k L 3yy 3yy J) 3yy 

By returning to Eq. (34) and comparing with Eq. (59b), we can determine that 


(62) 


A 

3 C rm 

3yy 


(e/wi E;r) ^ ( e ;m Ey r ) ^2 


(63) 
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so that all partial derivatives of direction cosine matrices in Eq. (62) are propor- 
tional to b 2 . With the substitution of Eq. (63) into Eq. (62), we find 


K if = -2 e ik b T J(e ik ~ en)b 2 [ 1* + 2 D kr b 3 F ' + 

tec? { L re® J 

~ 

~ r dF e ~ = 1 

+ 2 Dk,b 3 — b (e,-, — e jk ) b 2 F * 

~a L hi J 

- A [(«,* - «,.) bjb,#* + (e,o - «ifc)SA]| - -^7 

= ~2 ti*b T J(e ik - eji )b 2 1" T* + 2 D kr bsf r + — 1 
kt3> * L nSB J 


+ 2 

*as 


A 

~ r 3 F® ~ = i 

D k! b 3 1 + (ey« — £ jk )b 2 F‘ I 


- n 2 & 2 [(«,* - eyo) (*? - *“)] | — 


Since e, 0 = 0 for all f, and b T 2 h 2 (any vector) = 0, and b T 2 b 2 = 1, this expression 
reduces to 

K-a = ~ 2 e \kb\b 3 2 F) k i 

ktff 



+ n 2 2 «««/* (*? - *5*) - 

ibfcC? 


(64) 


Finally, we have found, in this meridional deformation case, with restricted inertias 
and all local deformations small, a class of system for which a symmetric K matrix 
is a realistic possibility. In particular, if external forces are zero and the hinges are 
elastically constrained as in Eq. (45), we can write Eq. (42) as 


Iy + Ky = 0 


(65) 


where K and I are symmetric matrices with elements given respectively for 
!,;'€(£? by 


Kij = n 2 2 e ik £ jk (*? - *“) + 8i jki (66) 

kif-P 

and, from Eq. (36d), 

lij =%i = b 2T 2 eireir<f> Tr b 2 
n&> 

2 2 e ir e it D, r Dr.(U~b 3 b T 3 )b 2 

rt& stff-r 

= 2 e ir ejr*r "OH 2 2 e iT £j a D, r D rt (67) 

rap ft# 1 
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In summary, we observe that in every case the linear, constant-coefficient equa- 
tions of motion which are to provide the basis for coordinate transformations have 
the form (see Eq. 42) 

ly + 7y + Ky = 0 


where I is a nonsingular symmetric matrix. 


Moreover, J = 0 and K is symmetric in the following special cases : 

(1) Hinge torques are elastic, and 

(2) External forces and external torques are zero, and 

(3) Either 

(a) The nominal spin rate is zero, or 

(b) Only rectilinear transverse deformations occur, or 

(c) Only meridional deformations occur and all augmented body principal 
axes are suitably aligned. 

Furthermore, J is skew-symmetric and K is symmetric for the general transverse 
deformation case under conditions (I) and (2) above. 


It should perhaps be emphasized that the conditions outlined have been shown 
only to be sufficient for the noted symmetry properties; they are not, in general, 
necessary conditions. 


From the basic relationships of analytical dynamics, it can be shown that for con- 
servative, holonomic systems, the linearized matrix variational equations obtained 
from 


d / dL\ 
dt / 



(i = 1, • • • , n) 


(68) 


(where L is the Lagrangian and <7, a generalized coordinate) always have the 
structure 

Mq + Gq +Nq=0 (69) 


where q = [q 1 ' ' q n ] T , M and N are symmetric, and G is skew-symmetric. It 

may be anticipated that these symmetry properties are also present for Eq. (42) 
for holonomic, conservative systems, despite the preference given in this report to 
an Eulerian rather than Lagrangian formulation of the equations of motion. 


In the following section, it will be demonstrated that the indicated symmetry 
properties of Eq. (69) or Eq. (42) can provide computational advantages in the 
determination of coordinate transformations which will permit these equations to 
be replaced by a set of uncoupled scalar equations. 


D. Transformation to Large- Deformation Modal Coordinates 

As established in the preceding section, the equations to be adopted as a basis 
for coordinate transformation for a given substructure can be written as the matrix 
second-order differential equation (Eq. 42) 


Iy + Jy + Ky = 0 
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where y is of dimension v by 1, I is nonsingular and symmetric, and, in particular 
cases, ] and K have special properties. This equation can also be written in the form 

Pf = QV (70) 


where, in terms of matrix partitions, we define 


-M p4 Eo-. L ^ 


Since 7" 1 exists, it is also possible to write Eq. (42) in the standard form 

r = Br 


(71) 


where 



1 0 

U " 

-I'K 

-i- 1 / E 


with U the v by v unit matrix. 


It is our objective in this section to find a new set of coordinates to replace 
y, y, and y in the partially linearized equations of motion. [See Eq. 30c, where 
y 2 , y\ • • , y v+l (for odd N) each represent column matrices containing flexible 
substructure deformation variables represented generically in this section by y.] 
The rationale for coordinate transformation is to minimize coupling among the 
new coordinates in the equations of motion, in order to make coordinate truncation 
feasible. As a practical compromise with this goal, we might instead seek trans- 
formations which uncouple the second-order differential equations in Eq. (42); but 
even this goal is not attainable unless the matrices 7 and K have special properties. 
We must therefore often be satisfied with a transformation which uncouples the 
first-order scalar differential equations in Eq. (70) or its equivalent, Eq. (71). 
Although it may be more convenient computationally to work with Eq. (70), 
thereby avoiding the matrix inversion of 7 in Eq. (71), the latter provides a simpler 
basis for discussion of the properties of the system. 


If the Jordan canonical form of the matrix B in Eq. (71) were nondiagonal, it 
would be impossible to find a transformation which would uncouple the associated 
scalar equations. Such a situation would, however, imply powers of t in the general 
solution of Eqs. (71) and (70), and this would in itself tell us that the dynamic 
response of the spacecraft is unacceptable. We therefore assume in what follows 
that the Jordan form of B is diagonal. In other words, we assume that there exists 
a similarity transformation matrix <t> such that 

S^BS- = X (72) 

where X is a diagonal matrix with elements X t , ■ • • ,X 2l .. Pre-multiplication of 
Eq. (72) by <£> produces 


B$ = <J>X 
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indicating that the columns <t>\ • • ■ , of <t> satisfy 

B& = = \.j& j = 1, ■ ■ ■ ,2v (73) 

so that A, is an eigenvalue of B and & is the corresponding eigenvector. The exis- 
tence of <t> _1 implied by Eq. (72) indicates that the eigenvectors of B must be inde- 
pendent for B to be diagonalizable. 


Post-multiplication of Eq. (72) by <t> -1 produces 

Q^B = A ^- 1 

so that the rows of <£‘ l are left eigenvectors of B; that is, if 4 ,lT , • • • , St 2vT are the 
rows of then 

^’ T B = \ i tr’ T = j — 1, • • • , 2v 

With this interpretation, it is convenient to replace 4> _1 by the matrix 'fr r , where 
the columns of 'fr are left eigenvectors ’Jr 1 , • • • , 4 r2v , normalized so that = U. 
It may be more convenient computationally to construct by assembling left eigen- 
vectors rather than by inverting <£. 


We can now diagonalize Eq. (71) with the transformation 


r = <t>Y 


(74) 


followed by pre-multiplication by 'Sr T . The result is 

* r <I>Y = ¥ t B & Y 


or equivalently, 


<t>- 1 <t>Y = <t> _1 B4>Y 


or 


Y = AY 


(75) 


Having satisfied ourselves that Eq. (71), and therefore Eq. (70), can be diago- 
nalized if the eigenvectors of B are independent (so that <I>~ 1 exists), we should 
note that for this independence, it is sufficient (but not necessary) that we have 
distinct eigenvalues. For if the smallest set of dependent eigenvectors of B is desig- 
nated <J> 1 , • • • ,<t> r , then they must be related by 


2 = 0 (76) 


in which none of the «, is zero. Pre-multiplication by B then produces (using 
Eq. 73) 


2 « ; B4> , — 2 aj\j& — 0 (77) 

i=x ;=i 
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Multiplying Eq. (76) by A r and subtracting Eq. (77) gives 


2 otj (A y - A r ) = 0 
j=i 

If (A, — A r ) =\= 0, then this equation contradicts Eq. (76), which asserted that the 
smallest set of dependent eigenvectors has dimension r. Therefore, the hypothesis 
of Eq. (76) is invalid unless the eigenvalues A, and A r are identical; that is to say, 
if the eigenvalues are distinct (nonrepeated), then the eigenvectors must be inde- 
pendent. The converse is not true; that is, eigenvectors can be independent even 
when eigenvalues are repeated. Nonetheless, we shall in what follows find it con- 
venient to concentrate on the special (and most common) case in which the eigen- 
values of the system of Eq. (70) or Eq. (71) are distinct. 


Now we can direct our attention to Eq. (70), which we may prefer to the equiva- 
lent Eq. (71) for computational reasons. Equation (70) admits a solution 


r = &e Xit 


(78) 


where <!>' and Ay satisfy 


(Q - A yP) = 0 (79) 

as may be confirmed by substituting Eq. (78) into Eq. (70). Equation (79) requires 

|Q — AyP| =0 (80) 

which produces 2v values of the scalar Ay (/ = 1, • • • , 2v). Because Q and P are 
real, complex roots of Ay appear as complex conjugate pairs. We refer to these 
scalars Ai, • • • , A 2 v as the eigenvalues of the differentia] operator in Eq. (70), and 
< t> 1 , • ■ • , <t ,2v are the corresponding eigenvectors. Since Eqs. (70) and (71) are 
equivalent, they both admit the solution in Eq. (78), and these eigenvalues and 
eigenvectors also belong to the matrix B in Eq. (71). In any case, the eigenvectors 
can be solved (from Eq. 79 or Eq. 73) only to within a multiplicative constant, and 
we might select this constant (normalizing eigenvectors) differently for Eq. (70) 
than for Eq. (71). 


If we introduce the transformation r = <t>Y (as in Eq. 74) into Eq. (70), we find 

P$Y = <?<t>Y (81) 


We now require a matrix pre-multiplier to take Eq. (81) into the uncoupled form 
of Eq. (75). 


We can obtain the necessary matrix formally 


p t v’ = g r r 


(82) 
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which is sometimes called the adjoint of Eq. (70). The eigenvalues A'- and eigen- 
vectors $■'' of Eq. (82) must satisfy 


(Q T - AjP T ) = 0 (83) 

and 

| Q T — X'jP T | = 0 (84) 

in parallel with Eqs. (79) and (80). Since the value of a determinant is unchanged 
by transposition, Eqs. (80) and (84) yield the same roots, and the eigenvalues of 
Eq. (70) and its adjoint, Eq. (82), are identical. Transposition of Eq. (83) produces 
(after changing the arbitrary index / to r and replacing X' r by A r ) 

V T (Q - A r P) = 0 (85) 

so that the eigenvectors 4> r ' of the adjoint equation are the left eigenvectors of the 
original Eq. (70) 

The significance of the left eigenvectors depends upon an orthogonality property 
which we can establish by pre-multiplying Eq. (79) by $> r ' T , to find 


$ r,T Q& = $,r' T p& X; 


( 86 ) 


Post-multiplying Eq. (85) by produces, for comparison, 


<t> r,T Q & = <t> r,T P<V\ r 


(87) 


Subtracting Eq. (86) from Eq. (87) yields 


0 = $ r "P<I>>' (A r - Aj) 


which, for A, =j= Ay, requires the orthogonality relationship 

^' T P^i = 0 r=j=; (88a) 

Equations (86) and (88a) now combine to produce 

<pr ,T Q$> =0 r =f= j (88b) 

If now we construct the matrix <t>' whose columns are • • • , < t >21 '', then 

the orthogonality conditions in Eq. (88) indicate that, if all eigenvalues are distinct, 
and are diagonal matrices. This means, of course, that <t>' T is the pre- 

multiplier we need in order to diagonalize Eq. (81). We choose to normalize the 
left and right eigenvectors in such a way that 


$ ,r P4> = U 


( 89 ) 


Then, the equation (from Eq. 81) 


<t>' T P<t>Y - 4>' r g$Y 
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must (by comparison with Eq. 75) become 

Y = AY (90) 

In summary of the general case, we have, demonstrated formally that Eq. (71) 
can be transformed to the uncoupled form Y = AY provided only that the eigen- 
vectors of B are independent, and we have demonstrated directly that the original 
Eq. (70) can also be transformed into Y = AY if its eigenvalues are distinct. For a 
formal proof that Eq. (70) can be transformed into Y = AY even when eigenvalues 
are repeated as long as eigenvectors are independent, we can rely upon the equiva- 
lence of Eqs. (70) and (71). 


Let us now direct our attention to special cases for which transformation is sim- 
plified. If, in the original second-order equation (Eq. 42), the matrix K (as well 
as I) is symmetric and / is skew-symmetric, then, in Eq. (70), the matrix P is sym- 
metric and Q is skew-symmetric. Equation (83) then becomes 

(-<?- x;P)4M' = 0 


or (since A) = Ay) 


(Q + AyP) = 0 (91) 

In comparing this result with Eq. (75), we first recall that the complex eigen- 
values occur in complex conjugate pairs (since Q and P are real). Next, we observe 
that the equation (Eq. 42) 


ly + Jy + Ky — 0 


possesses the first integral 


y T Iy + y T Ky = constant 


(92) 


as may be confirmed by pre-multiplying Eq. (42) by y 7 ' and observing that the 
scalar y T Jy must be zero, since its transpose is y T ] T y = ~y T Jy and any scalar must 
equal its transpose. The first integral in Eq. (92) is a reflection of the absence of 
nonconservative forces in the dynamical system represented by Eq. (42) in this 
restricted case. 


Equation (78) provides us with a solution of Eq. (42) when the former is written 
in partitioned form as 


A 

r = 


y 1 r ~\ A 

— - — — f- gM — 

L y J L _ 


that is to say, y has the solution 


y = 


(93) 
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Since y is real, and <f>’ and A, are generally complex, the conjugate pairs of eigen- 
values and eigenvectors must be so combined that pairs of solutions such as Eq. (93) 
appear as 


y = arf’ e x>t + a*rt’*e k * s ‘ (94) 

where a, is an arbitrary scalar multiplier and * denotes complex conjugate. To 
appreciate the reality of y in Eq. (94), and the possible complexity of y in Eq. (93), 

A 

let A j = aj + ioj and a, = c, + id jy and write Eq. (94) in terms of its real and 
imaginary parts: 


y = e a ’ 1 {( Cj + idj) (R (cos a jt + isin^t) 

+ (C; — idj) (Re<t>’ — i3m<t>’) (cos ojt — isin <jjt ) } 

= e a ’* {Re<j>’ [cj cos ojt + dj sin a jt] 

+ 3m<t>’ [dj cos ajt — Cj sin <x,t] } (95) 


Here i = ( — 1) % . 


The solution for y in Eq. (95) must satisfy the first integral in Eq. (92). If I and K 
are both positive semi-definite (or both negative semi-definite), then a,- = 0 in 
Eq. (95), and the eigenvalue A ,• is purely imaginary for all j. If I and K are not both 
positive semi-definite (or both negative semi-definite), then (by Liapunov’s theo- 
rem) the null solution of Eq. (42) is unstable. As previously, we reject the latter 
possibility (which would indicate unacceptable spacecraft performance). Now we 
can proceed with the comparison of Eqs. (91) and (75), knowing that for stable 
systems with the assumed symmetry properties of 7, J, and K, we have purely 
imaginary eigenvalues. Then, Eq. (91) becomes 


(Q - A *P) = 0 


indicating that <!>'' = and 




(96) 


Thus, in this special case, the left eigenvectors need not be calculated as an extra 
computational task, but may be recorded by inspection of the (right) eigenvectors. 

The second special case of Eq. (42) which is of practical interest requires 7 = 0 
and 7, K symmetric. In this case, we can work directly with the second-order differ- 
ential equations, constructing the v by v matrix cf> whose columns are <f> 1 , • • • , <j> v 
(see Eq. 93) and introducing the transformation 

y = <t>V ( 97 ) 

to obtain 

Itjyrj + — 0 (98) 
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Again we seek a pre-multiplier which uncouples these equations, and again we 
will turn to the transpose of the matrix of left eigenvectors. In this case, however, 
we can discover upon substitution of Eq. (93) into the restricted version of Eq. (42) 
that 


(XJJ + K) ft’ = 0 (99) 

so that an imaginary eigenvalue A, corresponds to a real eigenvector ft’. Moreover, 
ft’ is a left eigenvector as well as a right eigenvector, since, because of the sym- 
metry of I and K, we have (changing j to r for later convenience) 

ft rT (W + K) = 0 (100) 

If now we pre-multiply Eq. (99) by <ft rT and subtract the result from the post- 
product of Eq. (100) by ft’, we get 


ft rT Kft* - ft rT Kft’ + ft rT Ift* (A? — A?) = 0 


For A?=f= A?, we thus have the new orthogonality conditions 


and 


ft rT Ifti = 0 

r =H 

(101a) 

<ft rT Kft> = 0 


(101b) 


We further observe that pre-multiplying Eq. (99) by <ft iT produces 


ft’ T Kft’ = - a )ft’ T lft* 

which, by virtue of the imaginary character of A ; = iuj, becomes 

ft’ T Kft> = <r^ft’ T Ift’ (102) 

If we choose to normalize the eigenvectors in ft such that ft’ T Ift ’ — 1 for all /, 
then pre-multiplication of Eq. (98) by ft T yields 

■q + <r 2 r) — 0 (103) 

where 

r o-i o 1 

a 

a 2 — » 

_ 0 *o? _ 

We have proven this result formally only for the case of distinct (nonrepeated) 
eigenvalues, but the general argument preceding it applies whenever the eigen- 
vectors are independent. Although there are other special cases for which special 
transformations can replace the general procedure developed here for Eq. (70), 
these have been examined previously in this series (Ref. 21), and will not be 
re-examined here. It should perhaps only be mentioned that it is a common prac- 
tice among structural dynamicists to accommodate structural damping not by gen- 
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erating explicitly a symmetric matrix / for Eq. (42), but instead by modifying 
Eq. (103) to the form 


r j + 2£<rij + <r 2 r) — 0 (104) 

where £ is a diagonal matrix whose nonzero elements (the “modal damping ratios”) 
are established by experiment or by the judgment of an experienced analyst. This 
practice amounts to assuming the matrix 7 in Eq. (42) to be a linear combination of 
I and K, which seems to be a reasonable supposition for a wide class of structures. 


E. Hybrid-Coordinate Equations 

The matrix equation for unrestricted motion (Eq. 1) has been written as Eq. (30) 
for the special case in which certain of the angles of relative rotation and their 
time derivatives remain small enough to justify the neglect of terms above the 
first degree. In Eqs. (30a) and (30b), all linearized angles are collected in the 
single set ff\. In Eq. (30c), these equations are repeated in the form of a single 
matrix equation, and the variables in (A are subdivided into the column matrices 
y 2 , , y A+1 , for which in each case, the hinge labels are sequential. When all 

of the flexible appendages are classified as terminal appendages, the corresponding 
bodies should be given sequential indices, including the highest index in the sys- 
tem, so that in Eq. (30c) all small variables are in y 2 and all unrestricted angles are 
in y 1 . Finally, the partially linearized equations appear in Eq. (30d) as a first-order 
matrix equation, with the same grouping of linearized and unrestricted angles. 

In what follows, we let y k designate a typical submatrix of linearized variables, 
so k = 2, • • • , N + 1. Corresponding to each value of k, there is a linear, constant- 
coefficient equation of the form of Eq. (42), which we now write as 


l k f + /*y* -f- KV = 0 (105) 

We write the equivalent first-order equations (obtained either from Eq. 30 or from 
Eq. 105 directly) in the form of Eqs. (70) and (71), i.e.. 


and 


where 


p*r* = Q } t* 


r k = B k i* 



(106) 

(107) 


(108) 


In the general case, we introduce the transformation 


I* = (109) 

as in Eq. (74). The corresponding transformations for y k and y k individually are 
apparent from Eqs. (108) and (109). 
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In terms of the v by v unit matrix U and the v by v null matrix 0, we define the 
matrix operators 


SSo^tl/.'O]; Sou = [0 1 u ] (HO) 

From Eqs. (108)-(110), we have 

7 k = S^ = 25.^ (Ill) 

and 

y k = S 0 V* = (H2) 

These transformations can be substituted into the right-hand side of Eq. (30d), 
while the corresponding transformation in Eq. (109) is substituted on the left side. 


At the same time, we might substitute 



/ = 2£ 0 r* 

(J k odd) 

(113) 

and 



y k = 2J,i* 

(It odd) 

(114) 


into the right side of Eq. (30d) in order to achieve a more consistent notation. 
Eq. (30d) then becomes 


■Aoo | -Aoi I Ao 2$ 2 I 

Aio An J^A 12 $ 2 | • 

I I t 

A 2 0 I A 2 1 


A 22 <E> 2 


_a« _J. 

A 2 jv ^ A 2 ,y + ,* ff+1 


r l 

Y 2 

f 3 


Ajf+1,0 | Ajy+l, 1 I Ay+i, 2 4> 2 | 


I A s jr I A* ^ + i4> w+1 

J f 


X I A A t + 1 _y +l <j>^ +1 


ytf+i 


2 t(C°* + £ oi ) A* + C«&] - S a ok y k 

Jbefli ktC 

y. T r 1 




r y r r s "i 

^■‘OU ± 

R 3 


R N+1 J 


(115) 
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In anticipation of coordinate truncation, we would like to accomplish all pos- 
sible uncoupling of these equations of motion. We have no hope of uncoupling the 
entire set of equations, but we would like to come as close as possible to the point 
of uncoupling from each other the coordinates which define the oscillatory defor- 
mations of an individual flexible substructure. We know from Eq. (89) that 

Q’krpkQk _ u ( 116 ) 


and we can readily identify P k as the nominal value of A kk appearing in Eq. (115) 
for k even. Thus, we can obtain the equations of motion most nearly suitable for 
coordinate truncation by replacing Eq. (115) by an equivalent equation in which 
the row partition corresponding to Y k is multiplied by $>' kT > the transposed matrix 
of adjoint system eigenvectors. 


The result is 


_ _Aco_ _J Aox ( A 0i £ 2 _ j •_ _|_ ■Ao, 1 | 

A 10 ! _A„_ ' 'I 

4>' 2T A 2l ~~[ $' 2T A 22 $ 2 ! 

r 

; | ; ! ; ; 

$'‘ V4l, A» rtl o | <f>' JV+lT A w+1> j | *" r+ ‘ r A ir+llt **l* ■ •l«t>' yt,r A jV+1 , jm <}>^ 


_A»4P_ 

<f /2T A 2 „J 4> ,2T A 2l | <£ ,,2T A 22 $ 2 Jj • I $ ,2T A 2 , y+1 4> y+1 - 


r 1 

Y 2 


y^v+i 


" 2 [(C« 

_*c£B 

+ C ok ) A k + C° k h k ] - 2 flom " 
ksS 

f SLr 1 1 


1 

R 1 J 

— 


|~ 2,V*> 2 Y 2 J 

— 


L R * +1 J 


(117) 


The matrices are diagonal (/ even) when the corresponding matrix A n 

assumes its nominal value, P’. If its deviation from this value is small, and if it can 
be established that indirect coupling among the coordinate complex conjugate 
pairs in Y’ is small, then, for some purposes of dynamic analysis, it may be appro- 
priate to truncate the matrix Y’ to the smaller matrix Y’, preserving only those 
conjugate pairs of coordinates with frequencies in the domain of interest and with 
significant influence on those aspects of the dynamic response which are of interest 
in a particular case. The truncation of Y ; is accompanied by truncation of the cor- 
responding 2v by 2v eigenvector matrices <f> ; and to the rectangular matrices 
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<t>» and The result is a reduced-dimension set of equations for which the solu- 
tion preserves the salient dynamic characteristics of Eq. (117). These equations 
are recorded as Eq. (118). 


— 1 

o 

o 

A 0 i 

1 A 02 $ 2 

-i 1 - 

. 1 

— 

1 ~ 

Aio 

An 

A_12^ _ 

--U- 

& ,2T A 2 0 

$' 2r A 21 

!' 27 A 22 5 2 


1 

' 1 


H 

1 : 


_ V™ T Ay +1 ,o | 

*' N+lT Ay +1 


5 2 l • • 


• J_ _Ao^i^ +1 _ 


r 1 

f2 


yiv+i 






\+i,»+i 




r 2 [(C ok + C ok )A * + C ok A k ] - 2 a ok y k 1 

kiss keff 

L R 1 J 


rja^j 


SJ^Y- 1 1 

J 


A special case of Eq. (118) which is of particular interest is the case for which 
all angles and their derivatives may be assumed small. Recall that, in this case, 
it has been suggested that a chain of three imaginary massless bodies be introduced 
into the mathematical model, with the last of these (labeled S„) representing the 
“mean motion” reference frame of the system. In Eq. (117), the only remaining 
variables in this case are <u° and Y 2 , and the first of these represents the inertial 
angular velocity of the mean motion frame. If u° is retained as a system variable, 
then it is important that the “rigid-body mode” (zero-frequency mode) coordinates 
that arise in the transformed variable matrix Y 2 be discarded in the truncation to 
Y 2 because «>° and the rigid-body modes are redundant. 

A second special case of practical interest arises when in Eq. (105), the matrix J 
is zero and K is symmetric. As established in Section HID, in this case, the trans- 
formation in Eq. (109) can be replaced by the simpler transformation 

y* = (119) 

in parallel with Eq. (97). Differentiated terms in Eq. (30c) then become 

y* = <t> k v k y k = (120) 

and the transformed equations remain second-order differential equations. If such 
transformations are applied to all of the small variables in Eq. (30c), and appro- 


60 


JPL TECHNICAL REPORT 32-1565 



priate row partitions are multiplied by <f> 2T , ■ • • , <j> N+lT , the transformed equations 
have the structure shown in Eq. (121). 


a 00 + 5 00 


a 01 + & 01 


a 02 ^ ' 


a 0N + a 0 


( 

±_a™ 5^ ±J^L _ 2T_ [_ a^_+& ix j_ 


5 °.n+ 1 j 1 w + 1 


$ 2T (a 20 + a 20 ) $ 2T (a 21 + & 21 ) ? 2 T a 22 ? 2 |- <f> 2T (a™ + a 2 *) I <f, 2 T a 2 ’ y+1 J 

• i • i : ! • i : ! 






V 


_J^ _ J»Hh ^ 1 ' _J 5' VW + 2™ 1 &.K+tpM 


^ +ir (a w+1 -° + a y+1 ’°)| 4> N + lT (a N+1 ' 1 + S A,+, ' 1 )| ^»+» T gw+i,*^*l 


|^ + ir( 5 w +1 ,jr + 8»+i,W)j^r+i*3W+i,ff+i^W+i 


_f _ 


2 [(C ok + C° k ) A * + C ok fa] 

leffi 

fl l («°,r l , ‘ • ' » y w . y\ • * • , y", • 

— 2 Ookyk 

ktS 

• • ,’rj** 1 , fj x , ■ • 

• t) 

^ r R s (<Ay\ • • • , y w , y 1 , • 

• • y N r} 1 ‘ 
>7 > 7 > 

• • ' • 

■ t) 

R*(»W. ' • ' , y w , y\ • 

• • * v rj 1 ‘ 

• • ■ • 

■ 

_^ T R"+'(o,°,y\ ■ ■ ■ ,y»,y\ • 

• • n l • 

> 7 > 7 > 


• _ 


This result can be compared to the more general alternative displayed in Eq. (118). 
For computational convenience, Eq. (121) may be rewritten in first-order form, 
providing once again a set of equations having the structure of Eq. (118). In this 
special case, however, we can define the new variables as 

Y k = ( k even) and r* = °dd) (122) 

so that the transformation matrices become 



In this special case, for which Eq. (105) becomes Iy + Ky — 0 with K symmetric, 
there is no accommodation of energy dissipation in the mathematical model. This 
represents a departure from physical reality that can have serious consequences. 
Yet, it is a very difficult task to determine the proper damping characteristics to be 
assigned to any structural joint or connection, and it is virtually impossible to assign 
the proper characteristics of an analytical joint, introduced by the analyst to repre- 
sent discretely a structural deformation which is actually distributed over an elastic 
body. The difficulty of incorporating the damping characteristics of a nonspinning 
discretized flexible structure can be met with the introduction of modal damping, 
as that concept has long^ been used by structural dynamicists. When the vehicle is 
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spinning, the application of this concept is not entirely straightforward, although 
with some reinterpretation the basic idea can still be applied. 

For a nonspinning elastic body whose modal-coordinate free oscillations are 
described by ?j + = 0 as in Eq. (103), the introduction of modal damping takes 

the form of the incorporation of a viscous damping term in each of the uncoupled 
modal vibration equations, to obtain if + 2 £ 0 - 7 ? + — 0, as in Eq. (104). The 

analyst s task is to assign values to the diagonal elements • • • , £ v of the matrix f, 
on the basis of experiments or related experience. These values generally range from 
0.005 to 0.05 for the class of elastic bodies used in space vehicles. In terms of the 
equations of motion for the nonspinning multiple-rigid-body systems considered 
in this report, the incorporation of modal damping would involve the introduction 
into Eq. (121) of the matrix —2 Qo’rj’ as a part of ^R 1 , for / = 2, • • • , N + 1. 
This step is equivalent to the assumption that P in Eq. (105) is a polynomial in l k 
and K k (generally, P is treated as a linear combination of I k and IP). 


For a spinning system, Eq. (121) is replaced by Eq. (118), and the incorporation 
of modal damping becomes a more ambiguous procedure. The basic objective is to 
replace the imaginary eigenvalues of the undamped structure by eigenvalues with 
appropriate negative real parts. As shown in Ref. 21 (pp. 45 and 46), this amounts 
to replacing the quantities typified by A, — la, and A* = — ter,- in the diagonal 
matrix A in Eq. (90) by — £/ <r, + t<r, and — £,<r> ~~ X- The matrix A in Eq. (90) came 
from <t and in Eq. (118), neither A nor <t> ,T Q<t> appears explicitly. There will 
be found, however, on the right side of Eq. (118) a product which has the nominal 
value $' kr Q k >l> k = A\ so it is appropriate to augment Eq. (118) by the matrix product 


"fX o * 





Y k 



1 

b 

o 

) 




as a means of incorporating modal damping. It should be noted that this step no 
longer implies the relationship between I k , K k , and the damping terms in P which 
were implied in the nonspinning case. Here V represents the number of modal 
coordinates retained after truncation. 


F. Sample Problem Formulation 

I. Partial linearization. In this section, we shall extend the sample problem 
considered in Section HE (see Fig. 4) as an illustration of the unrestricted discrete- 
coordinate equations. Now we shall examine a set of partially linearized equations 
for this system, and explore the possibilities of coordinate transformation and 
truncation. 

In what follows, we shall assume that in Fig. 4, the angles y 3 , y 6 , y?, y 8 , and y 9 
remain small, while yj is a prescribed function of time and w° has a nominal value 
in the direction of g 1 . All other circumstances are as described in Section IIE. This 
is a rather complex example, possessing most of the difficulties of the general case. 
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In order to illustrate special cases, we shall occasionally consider restricted versions 
of this example. 


Equations (25) provide the first-order equations of kinematics and dynamics of 
this system without restriction on the size of any variables. These equations are 
equivalent to the kinematical equations in Eq. (24) and the dynamical equations 
in Eq. (1). 


For the partially linearized case, we replace Eq. (1) by Eq. (30c) as our starting 
point, retaining Eqs. (24) without change. Since y x is prescribed, y, does not appear 
on the left side of the version of Eq. (30c) to be written here. The variables in the 
differentiated column matrix on the left side consist of o>° and 

y 1 = [y 2 y3y 4 ] r ; y 2 = tys ye y? y 8 ye]; y 3 = [yio] (124) 


and the equations of motion are 


” a 00 + a 00 

a 01 + a 01 ! a 0 - 1 a 03 + a 03 ' 


* 0 

0) 

— ~ 

J 1 


— 

a 10 + a 10 

a 11 4- & 11 \ a 12 1 a 13 + a 13 



a 20 + a 20 

a 21 + fc^a^Ja 23 + a 23 


_Yl_ 

a 30 + 5 30 

a 3 ' + c \ 31 1 a 32 1 a 33 + a 33 


_ y 3 _ 


" 2 l(C ok + c° k ) A k + C°*&] - (floi + fl 0I ) V, 

—kfi& 

_R‘V,/,y 

R^./.^y^yhy^t) 

R 3 (<» # > y 1 . y 2 > y 3 > y 1 . y 2 > y 3 > *) 


(125) 


as a special case of Eq. (30c). Explicit values of the submatrices in the partitions on 
the left side of Eq. (125) must be constructed from the comparison of Eqs. (30a, b, c) 
and the substitution of expressions found in Eqs. (36) and (37). The matrices a 00 
and a 00 are available, respectively, from Eqs. (36a) and (37a), while all other matrix 
partitions must be expanded before the expressions in Eqs. (36) and (37) can be 
used. For example, the matrix a 01 + a 01 can be expanded as 

a 01 4- a 01 = [a„ 2 + a„ 2 j a 03 4- a „ 3 ! a„ 4 4- a M ] (126) 


and a 10 + S 10 is the transpose of this expression. Explicit values for the 3 by 1 
matrices in the partitions of Eq. (126) are available from Eqs. (36b) and (37b) 
As a second example, the matrix a 12 can be expanded as 




^26 

027 

O28 

029 

a 12 = 


036 

O37 

O38 

^39 


. 

0-4 6 

O47 

O48 

^49 


(127) 
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and the scalars in this matrix are available from Eq. (36d). A third example of 
particular interest is the 5 by 5 matrix 



O55 

o 56 

O57 

068 

O59 


^65 

Gee 

067 

o 6S 

069 

a 22 = 

O75 

076 

O77 

0?s 

O79 


Os 5 

086 

Os 7 

088 

Og 9 


_ O95 

Ose 

O97 

O98 

099 _ 


(128) 


the elements of which can also be obtained from Eq. (36d). 


In order to obtain explicit expressions for the matrix partitions on the right side 
of Eq. (125), one must recognize the definitions in Eq. (124) and then examine 
Eqs. (30a, b). The top 3 by 1 partition on the right side of Eq. (127) is identical to 
the right side of Eq. (30a) except for the transfer of the term (a ln + a m ) yi from 
the left side of Eq. (30a)— this because y! ( t ) is prescribed. The 3 by 1 matrix R 2 is 
from Eqs. (124) and (30b), given by 


g 5T 2 ^ [&*& + (C 6 * + &*) A*] + + 0 5 - (a 5 , + L) y\ 
g ,;T 2 [C*A* + (C 6k + &*) A*] + r 6 + & - (a 61 + o 61 ) yx 


(129) 


L g 97- 2 [C° k A k + (C 9 * + ch A*] + r 9 + ^ - (a 91 + L) y, J 

ktff 

Expressions for R 1 and R 3 are exactly parallel, differing only in dimension and 
numerical indices, as required by the definitions in Eq. (124). 


As in the sample problem formulation in Section IIE, it is convenient for subse- 
quent analysis and computations to rewrite Eq. (125) as a first-order matrix equa- 
tion, and to combine it with a set of kinematical equations, such as Eq. (24). This 
can be accomplished exactly as in Eq. (25), that is, by 


Vx = W 


(130a) 


where 


X- [# t 0 2 O 3 to° (a® a)® y 2 y 2 • • • y 10 y 10 ] r 

= [e r ja ) ° T !y lT !y 2r !y 3, '] :r (130b) 

l I I 

Because we intend in this section to restrict the results in Section IIE by partial 
linearization, we now express V and W in terms of partitioned matrices which 
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reflect the linearization. These results (to be compared to the unrestricted counter- 
parts in Eqs. 26 and 27) are given by Eqs. (130c). 


Vx = 


U I _0 [ 0_j_ _0_ j_0 |_0 j_0 j 0 _ 

0^5 00 + &>° 1 0 ijz 01 4- I 0 I a 02 ' 0_j_a° 3 + a 03 

_o _o_ r uf jO |~0 UTo | 0 _ 

Ola 10 -)- S 10 |o l a 11 + S 11 ' 0 I 0 1 a 13 + 


"oT 




0 0 , _0_ | U\ 0 | oj 


0 


0 a 20 + a 20 1 Oi a 21 + a 21 10 | a 22 ! Oj a 23 + & 23 

_0j 0 | 0 | 0 | 0 | 0j_ £j 0 


0 I a 30 + a 30 


0 I a 31 + a 31 


oT a 3 


0 I a 33 + a 3 


*•» 

y 3 


Po>° 


2 [(C ofc + C ° k ) A* + C*&] - (fl 01 + a 01 ) V. 

-W r, 

y 



R 2 


R 3 


(130c) 


If the partially linearized Eq. (130c) is to have any immediate computational 
advantage over its unrestricted counterpart in Eqs. (25)-(27), then this advantage 
is most apt to be found in the simplification of V, since it is the necessity of 
repeatedly inverting this matrix (actually or effectively) which makes the numerical 
integration of Vi = W such a costly operation. In order to minimize this task, we 
can in either case partition the matrix into constant and variable submatrices in 
order to take advantage of the algorithm for matrix inversion in terms of matrix 
partitions (see Ref. 27, p. 640, or illustration in Ref. 21, p. 49). 

As we examine the matrix V in Eq. (130), and compare with Eq. (26), we are 
reminded that the caret over a matrix is an indication that the matrix may depend 
upon the linearized variables, while the bar over a matrix indicates the absence of 
small-deformation variables; in both classes of matrices, however, we must expect 
variation with other kinematical variables. In order to detect the possibility of 
constancy of the coefficient matrices in Eq. (130c), we must examine the explicit 
expressions for their elements, as found in Eqs. (36) and (37). These elements must 
therefore be examined individually to establish conditions for their constancy. 

From Eq. (36a), the matrix a 00 (which is the same as a n0 ) is seen to have a defi- 
nitely constant part given by <t> 00 (the inertia matrix of the augmented reference 
body) and a generally variable part involving the direction cosines C°* for ked 1 . 
Only if all relative motions are assumed small (so C ok = U for all k) do we have a 
constant matrix a 00 . The matrix a 00 (or & 0 o) from Eq. (37a) is even more disagree- 
able, since it experiences time variations even for small deformations. If a>° is also 
small (even if <u° is large), terms such as a 00 a>° can be ignored; only then does the 
matrix a 00 + a 00 in V become constant. It is, in fact, clear that if all variables in x 
are assumed arbitrarily small (except perhaps 9, which is coupled only to a>°, which 
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may be large if nearly constant), and if there are no kinematically prescribed vari- 
ables of large amplitude, then the matrix V must be a constant, to be inverted 
only once in advance of numerical integration. This is, however, not the case for 
the example under consideration. 

Matrices in the second row-partition and second column-partition of V are avail- 
able in terms of Eqs. (36b, c) and (37b, c). These matrices also will not be effec- 
tively constant except in the extreme case noted, in which all kinematical variables 
are small (except possibly 6, 0, and a>°, which may be large providing that 8 and <a° 
are small). 

All other submatrices in V are composed of elements defined by Eqs. (36d) and 
(37d). The scalar a ik in Eq. (36d) can be constant only if the path elements such 
as ei r adopt zero values when multiplied by time-varying direction-cosine matrices. 
For example, the submatrix a 12 shown in Eq. (127) contains the scalar a :i9 , which, 
from Eq. (36d), is 

^30 - g 3T 2 e*re»rC 3 '<i>”C'°g° 
nff 

- QWg 3T 2 2 e tr *n (C 3r D’ rT C‘ r D rB - C 3s D sr D r ‘ T ) C' 9 g 9 (131a) 

re& stS’-r 


Recall from Def . 27 that the path elements are defined by 


A 


e S k — 


! 1 if /t g lies between 6 0 and 
0 otherwise 


Inspection of the configuration for this sample problem, shown in Fig. 4, now 
reveals that e 3r = 0 unless r = 3 or 4 and e 93 = 0 unless s = 9. (This information is 
also available from Table 2.) Thus, the first series of terms in Eq. (131a) vanishes, 
and we retain only 

a 3n = -Q7lg 3T [(C 33 D 93T C 9S D 39 - C 39 D 93 D 38r ) C 39 


+ (C 34 D 9iT C a *D* 9 - C 39 D 94 D 49T ) C 49 ] g 9 


(131b) 


Further simplification is afforded by the interpretation of Def. 35, which provides 
D 9 3 = D 94 = D 90 ; D 39 = D 30 ; D* 9 = D 40 


Moreover, Eqs. (32), (33) may be interpreted to provide 

C 33 = U ; C 93 = C 03 ; C 94 = C 04 ; C 39 = C 30 ; C 49 = C 40 

Thus, Eq. (131b) becomes 


a 3 o = -tyYlg 3T [(D 90T C 03 D 30 - C 30 D 90 D 3 ° T ) C 30 


+ (C 34 D 90T C 04 D 40 - C 30 D 90 D 40T ) C 49 ] g 9 (131c) 
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Evidently, in this case, the matrix in square brackets varies with time, so that only 
for very special configurations could a 39 be constant. (For example, if the matrix in 
square brackets has a principal axis remaining parallel to g 9 while g 3 remains per- 
pendicular to g 9 , then a 30 = 0; for the sample problem, we don’t have such good 
fortune.) 

Generalizing from this experience with a 3g , we can see that, except for very 
special configurations, matrices a'’ (i odd, j even), as typified by a 12 , are not con- 
stant, because in every case there is at least one unrestricted angle lying between 
each body in the set (P — (ft and any body in the set (ft. The same result emerges 
for a’ k (/, k odd), and the careted matrices all vary with time unless they are zero. 

Thus, the only matrix partitions in V in Eq. (130c) which retain the possibility 
of constancy are a 22 and (of course) the unit matrix U. As shown explicitly in 
Eq. (128), the elements of a 22 are the scalars found in Eq. (36d) for i, k e (ft. 
Because every term in a lk is multiplied by c ir , and i is confined to the set <71, we 
can discard all terms with superscripts r unless either re (ft or 6, is separated from 

by a joint whose index is in (ft. In application to the sample problem in Fig. 4, 
this rules out all r for which y r e y 1 , leaving only re (ft and r = 10. Moreover, the 
path element e*, is also zero unless either se< ft or 6, is separated from 6„ by a 
joint whose index is in (ft. Since s ={= r, we have in this sample problem confined $ 
to (ft. Thus, the expression for a,; c in this case becomes 

a* = g iT 2 e lr ek£ iT * ,r C rk g* + g i '«i.ioe» 1 ,oC i ' ,0 * 1# ' 10 C‘ # *‘g* 

rt(A 

- Wl g iT 2 2 (C ir D srT C ,r D r ‘ - C' s D sr D raT ) C rk g k 

rtCA SifA-r 

-Wlg iT 2 £ ir «)fe,io(C i ’'D 10 ' r7 'C 10 ’’'D’-' 10 - C i ' 19 D 10 ' r D’-' 19,, )C’--' c g fc 

- <Yng iT 2 ei,ioe». (C i ' l0 D s ’ 1 ° 7 'C , ' 10 D ll> ’® - C’®D*’ 10 D 10 '* T ) C 10 ’ h g k (132) 

All terms in Eq. (132) are constant except those involving C i,1 ° and C l0,k , and 
these are given for all i,ke(ft by C 10 * = C 10 6 and C’ 10 = C G1 ° = (C 10,6 ) 1 ’. From 
Eq. (4), we have 

C 10 6 = U COS y 10 - g 10 sin yio + g I9 g 10T (l - COS yin) (133) 

which obviously varies with time. Although many of these time-varying direction- 
cosine matrices are filtered out of Eq. (132) by the path elements, such terms 
survive when either ie {5,6} or ke {5,6}. Thus, the lower right 3 by 3 matrix in 
Eq. (128c) is constant. 

In generalizing from this example, we can conclude that the matrix a” is con- 
stant for even whenever the angles in y’ belong to a terminal appendage, and not 
otherwise except for very special configurations. For a flexible substructure which 
is not a terminal appendage, the corresponding matrix a’’ (;' even) can be parti- 
tioned into a constant part and a time-varying part. 

Partial linearization with the immediate objective of accelerating the numerical 
integration process would appear to be justified only if substantial and identifiable 
portions of the V matrix in Vi = W are replaced by constants in the linearization 
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process. For the sample problem, the V matrix is 24 by 24, and even without linear- 
ization, half of the terms are partitioned into constant (often zero) matrices. The 
suggested linearization in y 5 , • ■ ■ , y 9 renders constant only a 3 by 3 partition of 
a 22 , changing only nine of the 288 potentially variable elements to constants. This 
situation would improve significantly if y 10 and its derivatives were also assumed 
small, since then y 2 would embrace y 3 , the substructure would become a terminal 
appendage, and a new 6 by 6 constant matrix a 22 would emerge, changing 36 of 
the 288 potentially variable elements to constants. 

Partial linearization begins to offer significant advantages in computation when 
the number of bodies in a flexible substructure is large, and particularly when the 
substructure is a terminal appendage. In this same combination of circumstances, 
more dramatic computational advantage can be realized by transforming to dis- 
tributed coordinates, if one can justify truncation of high-frequency modes. 
Although these advantages are not very significant for the small sample problem 
chosen here, this will continue to be the vehicle for illustrating the general theory. 


2. Coordinate transformation. From Eq. (130c), we can extract 


a 22 f + S 20 <1° + & 21 f + & 23 'f = R 2 (a. 0 , y\ y 2 , y 3 , y\ y 2 , y 3 , t) ~ a 20 &° - S 21 y 3 - 5 23 y 3 

(134) 


All terms on the left side of Eq. (134) depend on y 2 or its time derivatives, as does 
the term R 2 on the right. Since coordinate transformations are to be based on a 
homogeneous equation in y 2 , which is obtained from Eq. (134), we must examine 
Eq. (129), in which we have explicitly displayed R 2 , and move those (careted) 
terms depending on y 2 or its derivatives to the left side of Eq. (134). If we then 
record only the homogeneous counterpart to Eq. (134), we have 


a 22 f + a 20 o>° + a 21 f + a 23 f 


g 5T 2 *5* [C 5 *& + - & 51 y\ 

ktff 


= 0 


(135a) 


A. A. — 


g 9T 2 *0* [C 9 *A* + C*A k ] +£ 9 - 8 01 y, 


htff 


In comparison of this particular example with the general theory in Section IIIC, 
we can identify Eq. (135a) as the specific counterpart of the general set of equa- 
tions established by Eq. (39). As demonstrated in the general arguments following 
Eq. (39), it is next appropriate to replace the linear, variable-coefficient equations 
in Eq. (135a) by a corresponding set of linear, constant-coefficient equations, which 
will have the structure of Eq. (42), that is, 


ly + /y + Ky — 0 


(135b) 


In this application, 


= [y. 


y« y7 ys 


y *] 3 
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nominal 


and the elements of J and K are as defined in Eqs. (43) and (44). 

More explicitly, the elements of the matrix I are established by the nominal 
(constant) values of the elements displayed in Eq. (128). Although any “nominal” 
might be selected, it may be most convenient to choose the values when all 
direction-cosine matrices are the unit matrix U, so that all angles for / e <g? — 
are zero. The element in Eq. (128) is then replaced by a ih available from the 
nominal value in Eq. (36d) as 


r A - 

j O'ij 


2 


'g 1 


Q4?g iT 2 2 e is (UD srT D'° - D°'D r ° T ) gi (i, j e (ft) (136a) 

rt£? stS’-T 


Similarly, the elements of ] in Eq. (43) become 

„ dk k 8rj 

Jij — —g‘ 2 ~ -s- 

keff °y > °y> 


(i, i £ (ft) 


(136b) 


for the indicated selection of the nominal case. With the substitution of Eqs. (49) 
and (46), and the recollection of items (7) and (8) from the input list at the begin- 
ning of Section HE, Eq. (136b) becomes 


hi = g iT 2 «i* [4>%W + £j k g i $ kk m <> + 7?<P kk r. ik g> 

ktff 

-cm 2 £js [D ks (g^° + 5°g0 D sk + ( UD* kT D k * - D sk D k,T )7Z°g’] ft d,8 t/ 

(136c) 


where 8 j ,• is the Kronecker delta, and <5° is the nominal value of w°. 

The elements of K in Eq. (135b) are obtained from Eq. (44) in similar fashion, 
incorporating Eqs. (46), (47), and (50) and replacing all direction-cosine matrices 
by U. The result is 

= -g iT 2 {(ejk - *»)W + Cm 2 D k ’ZWD’ k ) 

Uf? tfUS-k 

ft Ejk [(g'S° - +®°#**gV] 

ft cm 2 D k * [(*,. - e jk ) g^V* - £ ,-„ (g’S° - ^g')^ 0 

it® -It 

+ S» fttfej, g>)] D *} + Mi, 

or 

Ki j = -g iT 2 cikUiig’ (3°4> Wc 5° - 042 2 D k WD sk ) 
keff it®— fc 

ft £ jk [S°(^ Wc g ; - g’<t>“)«° + 0?2 2 (g’D ks - D k ’g’)^D ak ] 

it® — k 

ft cm 2 *, S D**S°S° g’D 5 *} + Mi, (136d) 

it® — lc 
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The next step is to find a coordinate transformation of Eq. (135b) which will 
permit this equation to be replaced by a system of uncoupled equations; for guid- 
ance in this procedure, we consider Section HID. In order to determine which of 
the coordinate transformations in Section HID apply to this sample problem, we 
must establish the structure of matrices J and K in Eq. (135b); we already know 
that I is real, symmetric, and nonsingular, and that / and K are real. 

In order to illustrate the two most useful coordinate transformations, this sample 
problem will be continued along two branches: (1) We shall assume the nominal 
value of o)° appearing in Eqs. (136c, d) to be zero, and the structural connections 
to be purely elastic, so that these equations become / a — 0 and Kij = fc.-S,,; and 
(2) we shall examine the most general case represented by Eqs. (136c, d). 

In the first case, Eq. (135b) becomes 


ly + Ky = 0 (137) 

as in Eq. (65). The transformation y = </>?;, recorded as Eq. (97), applies directly, 
and Eq. (137) can be rewritten as 


rj + <t 2 t) = 0 (138) 

as indicated by Eq. (103). Because a 2 is a diagonal matrix, the five scalar equations 
implied by Eq. (138) can be solved independently. If we were to consider inhomo- 
geneous counterparts to Eqs. (137) and (138), introducing forcing functions on the 
right side of both equations, then it would still be possible to obtain the five scalar 
solutions in Eq. (138) independently, and it might be feasible for certain purposes 
to approximate the total response in all five of the y; of Eq. (137) by 

y = <l>v (139) 

where rj and <f> are truncated versions of rj and $. In the extreme case, Jj might be 
a single scalar, perhaps 7) U and <f> the corresponding column of <f>, which, for 171, 
would be the first column. 

Coordinate truncation can be an enormously valuable device for improving the 
efficiency of dynamic analysis, but it can also be a dangerous oversimplification. 
In particular, for this sample problem, one must guard against eliminating a dis- 
tributed coordinate which might be excited to large and potentially destructive 
values by base motions of corresponding frequency, and one must avoid truncation 
of a coordinate whose modest response might degrade the pointing accuracy of the 
spacecraft or a spacecraft component. This sample problem is further complicated 
by the body 6 W (see Fig. 4), which can perform large rotations relative to the con- 
tiguous body 6 e - If 6 U) is an instrument with pointing requirements, it might neces- 
sitate the retention of appendage modal coordinates which affect 6 0 only negligibly; 
if 6 W is an actively controlled device, it might destructively force the appendage 
at the natural frequency of an otherwise insignificant mode. If 6 10 is massive and 
subject to gross departures from the nominal value assumed for the modal analy- 
sis, this will seriously jeopardize the entire truncation process. In any case, the 
truncation operation must be undertaken cautiously, with systematic evaluation 
of consequences. 
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Once some level of truncation has been chosen (perhaps tentatively), we return 
to the partially linearized equations (Eq. 130c), and substitute the truncated trans- 
formation y 2 = $rj (see Eq. 139 or Eq. 121), together with y 2 = $rj and y 2 = $rj. 
The fifth and sixth row-partition equations must then be multiplied by 4> r in order 
to obtain the reduced set of equations. We may observe that the fifth row-partition 
equation now says ~i T 4>Tj — 4> T 4>rj, and we can gain computational advantage by 
replacing this with Utj — Uij. We might wish also to incorporate modal damping 
in our system with the addition of —2 ~£orj to the term 4> T R 2 . 

In the more general case, for which <d° =j= 0 nominally, we have no simplification 
of Eq. (135b), so that coordinate transformation must be preceded by its repre- 
sentation in the first-order form 


Pr = Qr (140) 

as in Eq. (70). The appropriate transformation is 

r = 4>Y (141) 

as in Eq. (74) or Eq. (109). With the same truncation rationale as previously, we 
can instead substitute the approximation 

r = iY (142) 

into Eq. (130c), where r appears in the several forms 

|~-]=r; y 2 = [U;0]r; y 2 = [0 j U] T 

and rewrite Eq. (130c) in the form established by Eq. (118). 

Although the additional algebraic complexity of the more general case with a 
nominally spinning base is easily accommodated, there remain two considerations 
which make this case much more difficult to deal with than the nominally non- 
spinning case; the problem areas involve computation and modeling. 

The computational difficulty arises because the elements of the eigenvectors 
which comprise the columns of 5 in Eq. (142) are in general complex numbers, 
whereas the eigenvectors $ in Eq. (139) are all real. Computation of <t> is thus a 
significantly more time-consuming task than is the computation of although 
recent experience (Refs. 28, 29) suggests that even the general problem can be 
handled efficiently with the proper eigenvalue-eigenvector program. 

More troublesome problems must be faced in the initial construction of an 
appropriate multiple-rigid-body mathematical model of a spinning, flexible space- 
craft. The first major task is the discretization of the spacecraft. This involves not 
only subdividing the vehicle into rigid-body segments but also assigning the proper 
description to the connection characteristics. In the given sample problem, we must 
assign di and k, (i = 5, • • • , 9) as well as system geometry and mass properties. 
This is necessary whether the vehicle is spinning or not, but the assignment may be 
complicated by spin. 
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Unless there is actually a physical damping device at a physical hinge, it is prob- 
ably best to ignore damping torques at hinge connections until transformation to 
modal coordinates has been accomplished; incorporation of a modal damping 
matrix can then be accomplished both for spinning and for nonspinning structures, 
as described at the end of Section IIIE. 

The stiffness characteristics of each connection are defined by a single scalar, so 
that k 5 , ■ ■ • , k a must be specified for this sample problem. The determination 
of these spring constants is a routine chore when the structure is not spinning, 
although the calculation can be greatly complicated if variations in spacecraft 
thermal distortions result in variations in the undeformed nominal state of the 
vehicle. One must wonder, however, whether the structural properties of a non- 
spinning structure are changed by the steady-state forces borne by structural mem- 
bers when the vehicle is spinning. 

As noted in Ref. 25, it is necessary when dealing with a spinning elastic structure, 
described conventionally in terms of finite elements, to augment the elastic stiffness 
matrix with a new matrix called the geometric stiffness matrix or the preload stiff- 
ness matrix. The purpose of the new matrix is to reflect the influence on structural 
response of steady-state loads in the members— loads induced, in this case, by spin. 
Without such correction terms, there would, for example, be no stiffness in the 
direction of the spin axis for a finite element model of a spinning cable, and we 
know from experience that a cable, even without bending moment capacity when 
at rest, gains effective stiffness from spin. If the spinning cable (or a spinning 
beam) is modeled as a continuum, the partial differential equations include terms 
which reflect the foreshortening of the structure in the direction normal to the 
spin axis, and in this way the steady-state load in the structure changes the effective 
stiffness. We should examine the possibility of the necessity of such correction 
terms in the stiffness characteristics of a multiple-rigid-body model of a spinning 
flexible structure. 

Let us examine this question once again in the context of the cable disposed per- 
pendicularly to the spin axis. Now we idealize the cable as a chain of rigid-body 
segments, connected by hinge joints, as depicted in Fig. 11. We imagine that, when 
the vehicle is not spinning, the cable has no bending stiffness, so that in the absence 
of spin there are no springs at the joints. Should we artificially introduce “effective 
spring constants” when the vehicle spins? In this case, the answer is clearly nega- 
tive, since the foreshortening is automatically accommodated by the selection of 
variables of rotation rather than translation. 

Limited investigation suggests that the lesson of the segmented cable can be 
generalized. It is not necessary to augment the elastic stiffnesses of the hinge con- 
nections of a multiple-rigid-body model of an elastic structure to accommodate 
steady-state loads due to spin. Since the determination of the geometric stiffness 
matrix for a finite element model is a task of major proportions, this would seem to 
offer an extra advantage to the multiple-rigid-body model. 


IV. Summary and Conclusions 

A. Summary 

This report is confined in scope to the dynamic analysis of mathematical ideal- 
izations consisting of hinge-connected rigid bodies in a topological tree. The equa- 
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tions of unrestricted motion of the system are recorded as Eq. (1), in a matrix form 
that follows naturally from the previous developments of Hooker and Margulies 
(Refs. 9, 14). Several partially linearized forms of these equations appear as 
Eqs. (30a) through (30d), with the objective of providing approximations of Eq. (1), 
for which numerical integration is facilitated. Equations (1) and (30) are formu- 
lated in terms designed to minimize the demands on the user of a digital computer 
program for the numerical integration of these equations. This report includes the 
conceptual outlines of such a program and illustrative examples for its use; how- 
ever, we reserve for later documentation the detailed listing of the computer pro- 
gram and instructions for the user. 

The bulk of this report is devoted to the exposition of what seems to be an 
entirely new idea, which is based upon the application of a coordinate transforma- 
tion to a subset of the deformation variables in the partially linearized Eqs. (30), 
and subsequent coordinate truncation. The resulting hybrid-coordinate equations 
appear as Eq. (118). Because the distributed (modal) coordinates in this formula- 
tion are linear combinations of “small” variables which represent relative rotations 
of contiguous pairs of rigid bodies, their use does not imply that the deformations 
of the elastic substructure modeled by a chain or branch of elastically connected 
rigid bodies are “small” in any global sense, but only that the strains or local 
deformations are “small.” Thus, we have accomplished a discretized version of 
the equations of elasticity for large deformations with small strains. The distributed 
coordinates introduced here are called large-deformation modal coordinates. We 
expect these results to have applicability to a class of spacecraft for which deforma- 
tions are so large as to preclude the use of modal coordinates in the usual sense of 
structural dynamics (as in Ref. 21), and yet we can gain the advantages of computa- 
tional efficiency normally associated with hybrid-coordinate analysis (as opposed 
to discrete-coordinate analysis). 

While Eq. (118) provides an alternative to Eqs. (30) which has advantages in 
computational efficiency, the former equations demand a much more complex 
specification of input information, including mode shapes and frequencies gen- 
erated in a separate computer program for eigenvalue-eigenvector analysis. For 
this reason, the decision has been made at JPL to develop a distinct computer pro- 
gram for the numerical integration of Eq. (118), as opposed to its successively 
more comprehensive counterparts in Eqs. (30) and (1). This program will be docu- 
mented separately when it becomes available. 

B. Projection 

As noted in the Introduction, the determination of the influence of spacecraft 
nonrigidity on mission performance is not a single problem but a family of prob- 
lems, and a family of solution procedures is required for the efficient resolution 
of these problems. In this report, we have tried to deal with those vehicles for which 
a multiple-rigid-body tree model is appropriate. In a report now in preparation, 
we treat those vehicles for which it seems preferable to adopt an idealization con- 
sisting again of hinge-connected rigid bodies in a topological tree but with a small- 
deformation elastic appendage attached to each of the rigid bodies.* These equa- 
tions are still quite tractable, and are amenable to digital computer numerical 
integration. More general formulations will be developed if it becomes clear that 
they meet the criteria of need and reasonable efficiency of simulation. 


*The equations of motion appear in vector-dyadic form in Ref. 30. 
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Appendix A 

Derivation of Discrete-Coordinate Dynamical Equations" 

In addition to the symbol definitions adopted in Section IIA, consider the 
following: 

Def. Al. Let p k be the position vector of c k with respect to CM. 

Def. A2. Let the system of forces applied to by 6j for ke£B and / e £8* be equiva- 
lent to a resultant force f ki applied at the labeled point (/ij or /i k ) com- 
mon to 6k and 6j, plus a torque T ki . 

Def. A3. Let t ki be the kinematical constraint torque applied to 6 k by 6 h for kel 3 
and j e ffS k ; that is, let the total torque T l£i applied to 6 k by its neighbor 
6j be given by 

T*> = t ki + 8 iyk T k g k - 

so that if k > /, then g k • T k > = Tk , and if k < /, then g’ • T kj — — t,-. 

Def. A4. Let oo'- be the inertial angular velocity of 6k- 

All dynamical information regarding the motion of the n + 1 rigid bodies whose 
indices comprise the set £3 must be contained in the translational equations 

F'+Sf"-»i f (X + p') = 0 re£3 (A-l) 

and the rotational equations 

T* + 2 T*' + 2 P w X f k > - l*-cb* - to* X l fc -oo* = 0 Jfce£3 (A-2) 

;e®k fe£ B* 

The objective is to recast these 2(n + 1) vector equations into n + 6 scalar inde- 
pendent equations which do not involve unknown scalar components of the kine- 
matical constraint torques t ki or the constraint forces f ki . To this end, sum all equa- 
tions from Eq. (A-l) to find 

2 [F r + 2 f rs ] = 2 [m r X + m r p r ] 

rt£B st£Br rt£B 

The action-reaction principle (Newton’s third law) provides 

2 2 f r8 = 0 (A-3) 

«S5 W®r 

and the mass center definition furnishes 

2 m r p r - 0 (A-4) 

reffl 

so that the sum of Eqs. (A-l) becomes simply 

2 F r = 2 m r X (A-5) 

re£B re£B 

*The derivation following is the natural progeny of the derivations by Hooker and Margulies 
(Refs. 9 and 14), and represents a special case of that presented in Ref. 30. 
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or, with Defs. 25 and 30 from Section IIA, 


F = OftX (A-6) 

This result is the vector counterpart of the matrix dynamical equation recorded 
in Section IIB as Eq. (2), but in addition to verifying that equation, it permits 
Eq. (A-l) to be rewritten as 

2f rs =-F' + m r r2 5-+P v ] re£8 (A-7) 

*as, L tugs J 

Summing selected equations from the set implied by Eq. (A-7) yields 

f’ k = 2 2 f rs = “2 |"(F r - -^7 2 F") - m r p'l (A- 8 ) 

by virtue of the action-reaction principle. (Note that the indicated summation pro- 
duces the translatory equation of motion for the nested set of bodies consisting of 
all bodies on the branch attached to and commencing with 


The action-reaction principle combines with Eq. (A-8) to give 


pi = — ft* = 




kem (A-9) 


and substitution of this result into Eq. (A-2) eliminates all interbody forces from 
the rotational equations. 


The term of concern in Eq. (A-2) can, with Eq. (A-9), be written (using Defs. 32, 
34, and 35 from Section IIA) as 


2 P ki X t ki = 

jeffii 


2p w x 2 

je£8i refliki 




= 2 


L kr X F' - L fcr X 


mr_ 

Q 71 Z 


F' 1 - L*' X m 


ue£B 


rp r 


= 2 L kr X F r + D“ X 2 F“ - 2 L* r X m r p r 

re® re® 

= 2 D* r x F r - 2 L fcr X m r p r (A-10) 

re£B reffi 

It should be noted that, by Def. 32, L** = 0; thus, certain of the terms in the pre- 
ceding summations are zero. 


Combining Eqs. (A-10) and (A-2) furnishes the vector equations 

T*+ 2T‘ f + 2 D fcr X F r — 2L tr Xm r p v -l t -tf-rfxi ,[ *a t = 0 fce£3 

;e£Sk re® re® 

(A-ll) 

Next, sum individual equations from Eq. (A-ll) over those values of k corre- 
sponding to the indices of a nested set of bodies consisting of all bodies on a branch 
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attached to and commencing with the resulting equation, which by virtue 
of the action-reaction principle involves no interbody torques except T 8,Vi , is 
given by 


T «v. + 2 [T k + 2 (D*= r X F r - L* r X m f p r ) — \ k • — to* X I* • w*] = 0 se^ 

(A- 12) 

In what follows, 2 will be replaced by 2 which, by Def. 27, is equivalent. 
Definition A3 permits the substitution 

T 8Ar * = t 8S < + r s g* (A- 13) 

so that dot-multiplication of Eq. (A-12) by g s yields the scalar equations 
r« + g 8 • 2 [T & + 2 (D* r X F r — h kr X m r p r ) — \ k ‘ & k — oa k Xl k ‘ a> k ] =0 se£P 

ktff rim 

(A- 14) 

Equation (A-14) provides n scalar equations of the desired character, being free of 
any involvement with constraint forces or torques. The required three additional 
scalar equations are readily obtained by summing the individual vector equations 
represented by Eq. (A- 11) over all values of k in £8. Since the action-reaction prin- 
ciple in this case eliminates all interbody torques, the result is 


2 [T k + 2 (I> cr x F r - L kr X m r p r ) - l k ’j> k - 0 )* X !*•«*] = 0 (A-15) 

keffl re® 

In order to prove that Eqs. (A-14) and (A-15) establish the validity of Eq. (1) of 
Section IIB, one must first use available kinematical relationships to express p r , (*>*, 
and oX (for r, k e£3) in terms of the kinematical variables appearing in Eq. (1). 


To eliminate p r , let Q r > be the set of indices of bodies lying on the path between 
& r and 6 U and note, from Defs. A1 and 32, that 

p r - p> = L’ r + 2 (L er - L 8; ) - L r > (A-16) 

se(9rj 

See Fig. A-l for an illustration of this relationship. Having recorded Eq. (A-16), 
one can recognize more easily that the right side is unnecessarily complex; since, 
from Def. 32, L’’ = L rr = 0, and for any index s in the set £8 — (9 r j the sum L 8r — L 8J 
is zero (as in Fig. A-l for s = x,y, or z ), Eq. (A-16) can be replaced by 

pr - pi = 2 (L sr - L 8 >) r,/e£8 (A-17) 

st £ B 

with the obvious consequence, 



rrij 




L 8 ’) 


r,j : t£8 


(A-18) 
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For any given value of r, there are n + 1 equations implied by Eq. (A-18); 
summing these provides 


2 


nij 

<vn 


p’ 


2 9 ’ = 2 2 (L ,r — L*^) 

jdB V 11 K® y " 1 


(A-19) 


Equation (A-4) removes the second term on the left, and the definition of 0% 
(Def. 30) simplifies the first, yielding 


P r = 2 2^f(L--L*0 

m «£B ni 

Reversing the summation sequence and expanding produces 


(A-20) 


= sU"- 2&L"1 (A-21) 

aSL jc£B K '' 1 J 

Definitions 33 and 35 now permit the representation 

p r = 2 [L sr + D 88 ] = 2 D sr (A-22) 


The term in the dynamical equations (A-14) and (A-15) which requires the sub- 
stitution of p r from Eq. (A-22) is 


- 2 L kT X m r p r = -2 L* r X m r 2 D ,r = -2 L' cr X m r (D kr + 2 D 8r ) 

re® re® se® re® «®-Jt 

(A-23) 
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The first term in parentheses in Eq. (A-23) contributes (with the use of Def. 35) 
-2 L* r X m r D kr = 2 (D* 4 - D*') X m r D kr 

rtfR rt£3 

= D“X2 m&r - ~ ( 2 D k ' X m r D kr ) (A-24) 

re® dt ret 3 

The interpretation of b k as the mass center of the augmented body 6 ( permitted 
the recognition in Section IIA of the relationship 

2 m r D*' = 0 (A-25) 

which eliminates the first term on the right of Eq. (A-24). The second term involves 
the inertial time derivatives of vectors T> kr , which are fixed in so that, with 
Def. A4, one can write 

T)kr = o>*= X D* r (A-26) 

Equation (A-24) thus becomes 

d 

-2L*'X m r D Sr = --r[2D‘ , X (nirto^ X D* r )] 

re® dt 

= - { [ 2 rn r (D kr • D* r U - D* r D* r )] ' “*} 

at re® 

= --^(K k -a k ) (A-27) 

where U is the unity dyadic and the dyadic K fc is defined by 

K*=2w, (D kr • D kr U - D kr D kr ) (A-28) 

re® 

Thus, Eq. (A-23) becomes 


-2L*'X m r p r = - ~ (K k • «*) - 2 2 m r L kr X D sr (A-29) 

re® dt re£8 se®- k 

The second term on the right in Eq. (A-29) can be expanded to find (using Eq. A-25 
and L kk = 0) 

-2 2 m r L kT X D sr - 2 m r L kr X D rr 

—k $e£B—k—r fe£ZJ-k 

= -2 2 m r L fcr X D 8r + 2 2 m.L lr xi)'* 

reSi-k re®-k te®-r 

= -2 2 m,L*'XDn 2 2 mJJ r XT> rt (A-30) 

re®— k te®—k—r re®-k S(85-r 

where in the final step, the indices of the summed quantities in the first term have 
been exchanged, taking advantage of the fact that the double summation ranges 
over a set of indices r, s which is symmetric in the sense that a, b belongs to the 
set if and only if b, a belongs to it as well. This observation (first noted in Ref. 9) 
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permits the factorization of the terms in Eq. (A-30), and the recognition that 
li kk — 0 allows the further simplification 

2 [ 2 m s L kr X D rs - 2 m g L ki X D ri ] =22 m, ( L kr - L* s ) X D™ 

«£S-fc «f8-r xtB-k -t rdB-k st(B-r 

(A-31) 

The quantity L kr — L * 8 is zero for any index s corresponding to a body which 
lies anywhere on the branch that begins with 6^ kr and includes for any other 
index s, the quantity D r8 is also D r!r . Thus, D rs can be changed to D rt . At the same 
time, s e £8 — r is replaced by sefi, since there is no contribution when s — r. 
With Defs. (33) and (35), Eq. (A-31) then becomes 

2 2 m e (L kr - L* 8 ) XD r ’= 2 [2 m 8 (L/ Cr - L k ‘)] X B' k 

re ®— k st£B—r re®-k 

= 2 [042L' cr - 2 m s L ka ] X b rfc 

re® -Jr se® 

= 2 [QHL ir + (vnD kk ] x i) rk = (yyt 2 D kr x b rt 

r«£3-k re£B— fc 

(A-32) 


Finally, Eqs. (A-29)-(A-32) combine to provide 

-2 L* r X m r p r = - ~ (K k • o) fc ) + Q 11 2 D tr X D 1 * (A-33) 

re® re®— J: 

thereby permitting Eqs. (A-14) and (A-15) to be written in a more useful vector- 
dyadic form as 

T„ + g 8 • 2 e S k [T fc + 2D s, XF r - <t> kk • 0 ) k - to* X 4> u • w k 

kt{? re® 

+ 072 2 D tr X D 1 *] = 0 se<£> (A-34) 

re®— * 

and 


2 [T* + 2 D fcr X F r - • co fc - X • w" + 012 2 D* r X b rk ] = 0 

k«£8 refB rt®-k 

(A-35) 

where the relationship between K* in Eq. (A-28) and <t> kk in Def. 36 has been noted 
to provide K fc + l fc = <t> kk . 

Observing that D rk is fixed in 6 r , we can substitute 

D rfc = dr X D rfc + w X (co r X D rfc ) (A-36) 

The final sums in Eqs. (A-34) and (A-35) thus become 
047 2 D kr xi) rk = (Vn 2 { - D* r X (D rt X U') + D*' X [u r X (u' X D'*)]) 

re£B — k 

= 042 2 {— D rt D* r • d> r + D rfc *D fcr (i) r 

rcfB-Jt 

+ D ir X[o'X(w'XD rt )]) (A-37) 


JPL TECHNICAL REPORT 32-1565 


81 



The final form of the vector-dyadic equations of motion is available upon sub- 
stitution into Eqs. (A-34) and (A-35) of Eq. (A-37) and the kinematic expansions 


and 


co* = 0 ° + 2 e rfc y r g r 

riff 


(A-38a) 


cb* — cb° + 2 (g r yV + “’’ X g r y r ) (A-38b) 

Tiff 

The result of these substitutions into Eq. (A-35) may be written in the form (with 
judicious switching of index labels to permit later comparisons) 


% 2 2 (D>* * D^U - D'*D* ; ')> • <b° 

Iceffl iiS3 -k 

+ 2 [ < t >rr - 07l 2 (D ,r • D r, U - D jr D ri )] • 2 e* r g*y* = 

nff} ii(R-r luff 

2 {T* + 2 D k} x F> -<!>**• 2 e rk yr (o # + 2 e ar y*g a ) X g' 

IttCB >'e£B rei? 

- (co° + 2 £rtg r yr) X I*** • (co° + 2 £rtg r yr) 

rtf? 

+ 042 2 D*' X [(co» + 2 Er;g r yr) x { (co° + 2 trjg'yr) X D>*}] 

jlffi-k riff riff 

+ <m 2 (D» • D*’U - D>'*D*>) • 2 enyr (co° + 2 * 8 ry 8 g s ) X g r ) 

(A-39a) 

If this equation is written in the symbolic form 

a 00 • <a» + 2 a°*y* = 2 A* (A-39b) 

loS 1 fceCB 

then, by comparison with Defs. 38, 39, and 41, it becomes apparent that Eq. (A-39) 
is equivalent to the first three scalar equations recorded in Section IIB as Eq. (1). 
This recognition requires the identification of the matrix a 00 in Def. (38) as the 
representation of the dyadic a 00 in vector basis {b 0 }, i.e., a 00 = {b°} T a 00 {b°). 
Similarly, it must be recognized that a 0 * = {b °} T a ok and A* = {b *} f A*. 


The remaining scalar equations in Eq. (1) must be confirmed by Eq. (A-34) after 
substitution of Eqs. (A-37) and (A-38). With some switching of indices to facilitate 
comparison with Eq. (1), Eq. (34) may be written in the form 


g s * 2 £ S r [4> rr - 04? 2 (D ;r • D r, U - D ,r D r> )] • co° 

riff ii 33 -r 

+ g s * [ 2 e.r [®’ r - 071 2 (D ir • D ri U - D”D ri )] • 2 eyrg’y, = 

r iff iiffS-r jiff 


g 8 • 2 e sk (T* + 2 D* J X F> - O** • 2 e rk y r {(o° + 2 ^ysg s ) X g r 

luff il& Tiff siff 

- [(«* + 2 e rk g r yr) X O** • (u° + 2 ertg'yr) 

reZP red 1 

+ OK 2 X [(tt° + 2 erjg'yr) X + 2 trjgjr) X D**}] 

jlffl-k riff riff 

+ 071 2 (D J * • D* ; U - D>'*D*’) • 2 e Ti y r («° + 2 e 8 rg s y 8 ) + g r } 

ii&-k riff siff 

+ r s (se^P) 


(A-40a) 
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with the symbolic equivalent (for se£P) 


a 80 * “° + 2 o.yy’i = g s * 2 e.kA* + r 8 (A-40b) 

ji(P kiff 

Once the vectors in this scalar equation are recast as matrices and comparison is 
made with Eq. (1) and Defs. 39 and 40, it becomes clear that with Eq. A-40 we 
have established the validity of the last n scalar equations in Eq. (1). Thus, Eqs. (39) 
and (40) together constitute a proof of Eq. (1). 
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Appendix B 

Exact Scalar Equations for a Three-Body Example 


In order to reveal in explicit detail the structure of the scalar equations of motion 
derived in generic matrix terms in the body of this report, these matrix equations 
are expanded here for application to the special three-body system illustrated in 
Fig. B-l, in the absence of external forces and torques. 

Body A 0 has a planar slit through its center c in the b“, b® plane, and a dumbbell 
body 6\ has its center attached to A a at c, so that A 1 can rotate relative to A 0 only 
about a hinge axis parallel to b° (which is also labeled bj). 

Body An is a solid cylinder embedded in a cavity in A 0 , with the cylinder axis (and 
the single degree of rotational freedom) established by b® (which is also labeled 
b;)- Point c is the mass center of the combination A 0 + An, and principal axes of 
inertia of A 0 + A 2 are defined by b°,b“,b®. The point c is also the mass center of 
Ai, and hence of the system A 0 + A l + An. 

The equations of motion are given by Eq. (1) as 


#00 #01 #02 


P «o H 
1 1 ) 


- A® + C 0X A X + C" 2 A 2 “ 

#10 #11 #12 


Vi 

= 

U lT A x + r r 

#20 #21 #22 


_ 72 _ 


U lT A 2 + r. 


A 

where the symbol u x = [1 0 0] T represents g 1 and g 2 , since g 1 - g 2 = b® = b) = b 2 . 
From Def. 41 in Section IIA, we have simply 


A® = -S 0 $®®c® (B-2) 

because both A , and A 2 are supported at their mass centers, which are also their 
barycenters, so that D 20 = D 10 = 0. 



Fig. B-l. Three-body example with mass-center 
connection points 
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The 3 by 1 matrices A 1 and A 2 are, from Def. 41, 


A 1 = -<t> 11 y 1 C lo S°C 01 U 1 - (C 10 S°C 01 + yjU 1 ) 4> n (C 10 (D° + yjti 1 ) (B-3) 

and 

A 2 = -$ 22 y 2 C 20 35°C 02 n 1 - (C 20 S°C 02 + yjti 1 ) 4> 22 (C 20 0 )° + y.U 1 ) (B-4) 

where advantage has been taken of the identities tP-tt 1 = 0 and D 21 = D 12 = 0. 

Further simplification of A 1 and A 2 is afforded by the substitution from Eqs. (4) 


and (5) of 

C 01 = U cos yj + u 1 sin yj + u 1 u lT (1 — cos yi) (B-5) 

C° 2 = [/ cos y 2 + n 1 sin y 2 + tFu 1 * (1 — cos y 2 ) (B-6) 

C 10 = C otT C 20 = C 02T (B-7) 

which provide 

C^U 1 = U 1 COS yj + 0 + U 1 (1 — COS y,) = U l (B-8) 

C 02 !* 1 = U 1 COS y 2 + 0 + U 1 (1 — COS y 2 ) = M 1 (B-9) 

With Eqs. (B-8) and (B-9), A 1 and A 2 become 

A 1 - — y l 4 ,ll C l °co 0 U l - (C 10 S°C 01 + yiM 1 ) (C 10 W ° + y lM l ) (B-10) 

A 2 = -y 2 4. 22 C 2<, <o°u 1 - (C 20 w°C 02 + y 2 tT) <E> 22 (C 2 V> + y.u 1 ) (B-ll) 


Before expanding Eq. (B-l) in scalar terms, we must examine the definitions of 
the elements of the coefficient matrix on the left side. From Defs. 38-40, we have 

a 0 o = <E> 00 + C 01 4> U C 10 + C 02 4> 22 C 20 (B-12) 


#01 — 

C 01 <t> n u 1 

a 10 — u ir 4> ll C 10 

(B-13) 

#02 

C 02 $ 22 w l 

a 20 = M ir $ 22 C 20 

(B-14) 

> 11 w 1 ; 

#12 = 

a 2i = 0; a 22 — u lT <t> 2 -*u l 

(B-15) 


If is a dumbbell of length 2 L, with tip particle masses equal to m and a 
massless connecting rod, then 




2mL 2 

0 

0 


0 0 
0 0 
0 2 mL 2 


(B-16) 


since (from Def. 36) 4? 11 = Z u plus terms which in this case are zero. 
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Similarly, for the cylinder A,, we can write 


$22 = 


! 

0 

0 


0 0 
Y 0 
0 /' 


(B-17) 


where / and /' are moments of inertia of the cylinder for the symmetry axis and a 
transverse axis, respectively. 

Use of Eqs. (B-6), (B-7), and (B-17) provides, upon expansion, 

C 02 4> 22 C 20 = <t> 22 (B-18) 


so that Eq. (B-12) becomes 


a 00 = <£ 00 + <t> 22 + C 01 <t> n C 10 (B-19) 

Because c is the mass center of A 0 + A>, c is also the barycenter b 0 . Since $ 00 is 
the inertia matrix of the augmented body A 4> 00 + $ 22 is the inertia matrix of the 
system A» + A 2 about the system mass center at c. Since by hypothesis, b°,b°,b° are 
principal axes of this system, we can write 


$00 _ f _ $22 _ 


h 0 0 
0 h 0 
0 0/ 3 


(B-20) 


where I u I 2 , 1 3 are principal axis moments of inertia. 

Combination of Eqs. (B-16), (B-5), and (B-7) produces, upon expansion, 


C 01 $ n C 10 = 2 mU 


1 

0 


0 


0 


sin 2 -/, — sin yj cos yi 

COS 2 y! 


0 — sin yj COS yi 

With this information, we can expand Eqs. (B-12)-(B-15) to obtain 


(B-21) 


floo 


h + 2mL 2 0 0 

0 l 2 + 2mL- sin 2 yj — 2mL 2 sin cos yi 

0 —2mL- sin y, cos yi I 3 + 2 mL 2 cos 2 y x 


(B-22) 


®01 


2 mL- 
0 
0 


= 2mL 2 u 1 


(B-23) 


Cjo = [2mL 2 0 0] = 2mL 2 u lT 


(B-24) 


^02 


/ 

0 

0 


— Ju 1 


(B-25) 
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a 2 o = [J 0 0] = }u lT 


(B-26) 


fl u = 2 mL 2 ; — a 21 = 0; a 22 = / (B-27) 

Having expressed the left side of Eq. (B-l) wholly in scalar terms, we can return 
to Eqs. (B-2) and (B-lO)-(B-ll) and focus again on the right side. Now we can 
utilize Eqs. (B-16)-(B-21) to simplify the upper partition to 

A 0 + C^A 1 + C 02 A 2 = -WV - y, (C 0I $ u C 10 )7S°u 1 - 7S° (C’W'C 10 ) o>° 

- <o°C ol $ 11 y 1 a 1 - C 01 y 1 U 1 <t >,, C l0 o>° - y lC 01 U 1 <f> it U 1 

- y, (C 02 3> 22 C 20 ) S “U 1 - S° (C 02 4> 22 C 20 ) o>° - 73°C #z $ a2 y 2 « 1 

- y 2 C <,2 U 1 $ 22 C 2 V' - ylC 20 U 1 4> 22 M 1 

= -Z° (<t> 00 + <t> 22 + C 01 $ n C 10 ) m 0 

- y, (C 01 <£ 11 C 10 5* > U 1 + WC 01 ^ 1 !* 1 + C 01 U 1( t ,11 C 10 <j) 0 
+ y,C ,1 U 1 <f> n U 1 ) - y 2 ($ 22 S°« 1 + w°C 02 <t> 2 V 

+ C 02 U 1 <t> 22 C 20 <fl° + y 2 C 20 U 1 cE> 22 U 1 ) 

= ~-«°a 00tl > 0 + Vi (C 0, * 11 C 10 8 1 + 2mL 2 u 1 - C 01 u'$ n C 10 ) <o° 

+ y 2 (<t> 22 5' + /u 1 - C 02 u 1 <5 22 C 20 ) 

But 

C 02 u‘$ 22 C 20 = /'C 02 u 1 C 20 

= /' (u 1 cos y 2 + u'u 1 sin y 2 ) ([/ cos y 2 — u 1 sin y 2 ) 

= /' (u 1 cos 2 y 2 — u'u'u 1 sin 2 y 2 ) = J'u 1 (cos 2 y 2 + sin 2 y 2 ) 

= J'u 1 

so that 

$ 22 m 1 - C 02 M 1 3- 22 C 20 = /'w 1 - J’u 1 = 0 
Moreover, the identity 

- (C 01 u 1 $ ll C 10 ) 5 ' = Co'Q^iPC 10 = C 01 $ ,1 C 10 u 1 
permits further simplification. The final matrix expression is 

A 0 + C 01 A' + C 02 A 2 = -SOflooco 0 + (2C 0l 4> 11 C 10 + 2mL 2 t/)y 1 uV + /y 2 uV 
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Substitution of the scalar expansions of the coefficient matrices (see Eqs. B-21 and 
B-22) yields 


A 0 + C^A 1 A C 02 A 2 = 


0 W 3 J 2 A 2ml/ s (( 03 S A o> 2 c) 02^3 2 t tiU 2 (o> 2 c 2 A < 03 sc) 

— o) 3 (Zi A 2 mL 2 ) — (ou A 2yj) 2 mUsc o>J 3 — Jy 2 A 2mL 2 (o>iC 2 — 2^s 2 ) 

0)3 (/) + 2mL 2 ) + Jfc + 2mL 2 (2yjC 2 — o^s 2 ) 2mL 2 $c(« 1 + 2yi) 


0)1 

0) 2 

0) 3 


(v 

a, - 


~~ I3 ) <t>2<*>3 A 2mL 2 [(co! — o>l) sc A 0>20>3 (s 2 - c 2 )] 
- Z x ) 0 ) 30)3 — Jy 2 0) 3 2rnL 2 ( 0)1 A 2yi) s (o> 3 S A 

* Io) fdiWo A ~1~ 2jyiZj 2 ( 0)1 H” 2v^ r -4- 


a>2 


>.c) 


(B-28) 


where s denotes sin y x and c denotes cos y,. 

Returning to the original equation of motion, Eq. (B-l), we see that all that 
remains of the task of obtaining scalar equations is the determination, from 
Eqs. (B-10) and (B-ll), of the scalars u lT A x and u lT A 2 . This operation is simplified 
by the relationships 

u xT u l - 0; u' T| i> n = 2mL 2 u lT ; u lT $ 22 = Ju lT 

U' r C 10 = U lT COS y! A U lT (1 — COS yj) = U lT ; U lT C 20 = U lT 

u lT 7S°C 01 u l = u'Vw 1 = 0; &'u' = 2 mL 2 u} 

u lT 'S 0 C 02 u 1 = u lT S°w 1 = 0; <l' 22 u 1 = Ju 1 

With these substitutions, Eq. (B-10) provides 

« lT A l = -y, (l/ , ’V 1 C 10 a>°U 1 ) - U ir (C 10 S°C 01 A yjii 1 ) (C 10 0 )° A y 1 U 1 ) 

= -« lT 8 B (C°‘* 11 C 1# ) o. 0 (B-29) 

and Eq. (B-ll) provides (noting Eq. B-18) 

« lT A 2 = -yziu'^&otfu 1 ) - « 1T (C 2 VC 02 A y 2 n 1 ) <I? 22 (C 20 o>° A frit 1 ) 

- -u ,T «5°$ 2 V (B-30) 


Substitution into Eqs. (B-29) and (B-30) of the scalar expansions given by 
Eqs. (B-17) and (B-21) permits the representation 

« lT A> = [0 «>3 ~o> 2 ] (C 01 * 1 ^ 10 ) o° 

= 2 inU [0 s (<u 3 s A « ) 2 c) — c (u> 3 s A w 2 c)] o>° 

= 2 mil [(<o| — o>l) sc A ci) 2 a >3 (s 2 — c 2 )] (B-31) 

and 

u lT A 2 = [0 03 -« 2 ]$ 22 o > 0 =J’[Q «»3 — a, 2 ] to 0 — 0 (B-32) 
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Finally, we can combine Eqs. (B-22)-(B-28), (B-31), and (B-32) in Eq. (B-l), to 
obtain the required set of five scalar equations of motion in the form 

(Z 2 + 2mL 2 )<l 1 + 2mL 2 y 1 + Jy 2 = (Z 2 — I 3 )ta 2 <a 3 + 2 mU [m,m 3 (sin 2 y 2 — cos 2 y 2 ) 

+ (<a 2 — 0 >|) sin yi COS y 2 ] (B-33) 

(I, + 2 mL 2 sin 2 y 2 ) <J 2 — 2mL 2 <1 3 sin y 2 cos yi = (Z 3 — Z,) u> 3 a>i — 7y..m 3 

— 2mL 2 + 2yx) sin y 2 (<b 3 sin y, + cos y 2 ) (B-34) 

(Z 3 + 2mL- cos 2 y 3 ) w 3 + 2mL 2 i> 2 sin y 2 cos y, = (Z, — Z 2 ) u^w-, + Jy 2 <*, 

+ 2mL 2 (on + 2y 2 ) cos y 2 (» 3 sin y, + u>-. cos y,) (B-35) 

2 mL 2 (»! + y 2 ) = 2 mL 2 [(a> 2 — M |) sin y 2 cos y 2 + to-, o >3 (sin 2 y 2 — cos 2 y i)] + t, 

(B-36) 

] (<^1 + y 2 ) = t 2 (B-37) 

Needless to say, it is not recommended that Eq. (1) be used (as in this appendix) 
for the purpose of deriving equations of motion for simple three-body systems 
such as that illustrated in Fig. B-l; the final scalar equations of motion obtained 
here can be derived from first principles by an experienced analyst in much less 
than half the time it will take him to execute the steps outlined in this appendix, 
obtaining the special set of equations from the generic matrix equations. Even for 
a much more complicated system, scalar equations of motion can probably be 
obtained in literal scalar form more efficiently by starting from first principles, 
although as the complexity of the system grows, the question becomes moot. 

The purpose of the generic matrix equations of motion in this report is not to 
give the analyst a better starting point for the derivation of equations of motion 
but to give the digital computer a set of equations in a format it can readily process. 
At the present, the computer is programmed to obtain the scalar equations from 
the given matrix equations numerically, performing the necessary multiplications 
and additions after each integration interval. In the future, as digital computer 
capacity for the manipulation of literal algebraic symbols is improved, it may 
become desirable for these manipulations to be performed symbolically only once, 
in advance of all numerical integrations. At no time is it intended that the analyst 
perform the manipulations, unless (as in this appendix) it is appropriate for demon- 
stration and checking. Equations (B-33)-(B-37) were derived also from direct time 
differentiations of angular momentum vectors, in order to establish some measure 
of confidence in the general equations, and to develop some subjective apprecia- 
tion of the analytical labor involved in each of the alternative derivations. 
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Appendix C 

Linearized Scalar Equations for a Three-Body Example 


In Appendix B, we considered a relatively simple example of a three-body sys- 
tem, and obtained from Eq. (1) a complete set of five scalar equations of motion 
for the system. Now we consider another three-body system with five degrees of 
freedom, and in this case (see Fig. C-l), we no longer make the support points of 
and on coincident with the corresponding mass centers, as they were for 
the preceding appendix (see Fig. B-l). 

The objective of this appendix is to display the structure of Eq. (1) in this par- 
ticular application, noting the complexity introduced by the off-center mounting 
of and 62 on 6 n , and then to extract from these general equations of motion a 
new set linearized in the variables y, and y 2 and their derivatives. These results 
are then to be compared to those obtained from the generic linearized equations 
developed as Eq. (30) in this report. 


In this example (see Fig. C-l), is a cylinder of radius r, and and are 

identical uniform, thin rods of mass m and length 2 L, attached to by hinges 

through points p , and p 2 . Unit vectors b), b", b° are fixed in / 0 , with b“ along the 
symmetry axis. Position vectors from the mass center c 0 of to /u and /i, are 
given, respectively, by 


P°1 _ r b« p02 = — 


(C-l) 


Thus, the barycenter b n coincides with c„, and we have 


D 00 - 0; D 01 = L 01 = p n, ; D 02 = L 02 = p 02 


(C-2) 


In and 6 2 , the barycenters are displaced from the mass centers, such that 


D 1 = (L - R)bl; D'° = - Rbl ; D 22 — — (L — R)b|; D 20 = Rb| 

(C-3) 

where R = mL/Qli. 



Fig. C-l. Three-body example 
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The hinge axes are defined by 


g 1 = g 2 = b“ 


and the corresponding elastic hinge torques by 


(C-4) 


Tjbj = — fcy,bj; Tjb'f — — ky 2 bl 


(C-5) 


The vectors in Eqs. (C-l)-(C-5) appear in the equations of motion in this report 
in the form of matrices or scalars. In this example, it is convenient to express many 
of these quantities in terms of the symbols 


A 

1 " 

A 

■ 0 

A 

' 0 

M 1 = 

0 

0 

; U 2 = 

1 

0 

; U 3 = 

0 

1 


(C-6) 


Then, some of the quantities of interest become 

D°° = 0 D 01 = ru~ D 02 = -rt/ 2 


D 11 


U 2 ((yyi _ 


R 


m 


D'° = —Ru- = D 12 

D 

D 22 = —u 2 (Qlt — m) ■ — • 

m 

D 20 = Rtr = D 21 


S' = g 2 = «' 

r i = — fcy, t 2 = ~ ky 2 (C-7) 

The equations of motion are given by Eq. (1) as 


000 

0OI 

002 


G>° " 


■ A° + C 01 A’ + C 02 A 2 ■ 

010 

0J1 

012 


Vi 

= 

U lT A l — fcy a 

_ 020 

021 

022 


_ 72 _ 


M ir A 2 - fcy 2 


corresponding exactly to Eq. (B-l) for the simpler system in Appendix B. 
From Def. 41 and Eq. (B-7), 

A® = ~sr°<t>°V + OW [ ru 2 C 01 (C 10 s°C 01 + ^t/ 1 ) 2 (-Rt/ 2 ) 
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- rU 2 C 02 (C 20 ra°C 02 + y 2 t?) 2 (Ru 2 ) 

+ (—Uu 2T C 10 U 2 + C 01 U 2 U 2T ) rRyj (<5°C 01 + ^C 01 ?/ 1 )!/ 1 
+ (- Uu- T C 20 u 2 + C 02 u 2 u 2T ) rRy 2 (S°C 02 + y^U 1 ) u 1 


(C-9) 
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This equation simplifies immediately upon noting that 


uHA = 0; tPu 2 = u 3 ; u 2 u 2T — [0 j u 2 ! 0] = U 2 
Further simplification is afforded by the substitution from Eqs. (4) and (5) of 

C 01 = U cos y, + 2 1 sin yi + i Au lT (1 — cos yx) (C-10) 

C 02 = U COS y 2 + 5 1 sin y 2 + u 1 u lT (1 — COS y 2 ) (C-ll) 

C 10 = C 0lT ; C 20 = C 02T (C-12) 

where u 1 u lT = [u 1 j 0 j 0] = U 1 . With the substitutions 

u 2T C 20 u z — u 2T u 2 cos y 2 — u 2T u 1 u 2 sin y 2 + u 2T u l u lT u 2 (1 — COS y 2 ) = COS y 2 (C-13) 

u 2l C 10 u 2 — u 2T u 2 cos yi — util'll 2 sin yi + u 2T u 1 u 1 ' r u 2 (1 — cos y,) = cos yi (C-14) 


C 0 i u 2 u 2T — C 0i U 2 = cos yj U 2 + sin yi u x U 2 (C-15) 

C 02 u 2 u 2T = C 02 U 2 = cos y 2 L7 2 + sin y 2 u'U 2 (C-16) 

C^U 1 = u 1 cos yx + u 1 (1 — COS yx) = u 1 (C-17) 

C 02 !! 1 = u 1 COS y 2 + U 1 (1 — COS y 2 ) = U 1 (C-18) 

C 02 U 2 = U 2 COS y 2 + U 3 sin y 2 (C-19) 

C 0 l U 2 = U 2 COS yi + U 3 sin y, (C-20) 

« 2 C 01 = U 2 COS y! + U 2 U l sin y, (C-21) 

5 2 C 02 = u 2 cos y> + u 2 BP sin y 2 (C-22) 

and such operations as 

tPC’WlPu 2 = ( U 2 COS y 2 + TPu 1 sin yi) tPll 3 

= — u 2 « 1 M 2 siny 1 = — « 1 siny 1 (C-23) 

we can rewrite A° in Eq. (C-9) as 


A 0 = — 3°$°V’ + Q4?rR [u 1 sin y 2 yf — w 2 © 0 ^ ( u 2 cos y t + u 3 sin yi) 

— M 2 3°C 01 U 3 y 1 — U 2 C 01 « 1 C 10 w° (« 2 cos yx + U 3 sin yi) yi 
— U 2 S°£° (u 2 COS y 2 + u 3 sin y 2 ) — (u 2 cos y 2 + iPu 1 sin y 2 ) U l U 3 y 2 
- M 2 S°C 02 U 3 y 2 - U 2 C 02 t?C 20 S° (u 2 COS y 2 + U 3 sin y 2 ) y 2 
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+ ( — (7 COS yi + C/ 2 COS yx + tP[/ 2 sin yx) yx?3 0 « 1 
+ ( — U cos y 2 + U 2 cos y 2 + iPt/ 2 sin y 2 ) y^U 1 ] 
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With the further substitutions 


C 02 u 3 = u 3 cos y 2 — u 2 sin y 2 (C-24) 

C 01 u 3 = u 3 cos yi — u 2 sin y, (C-25) 

u 2 C 01 u 1 C 10 = (ft 2 cos yj + iPu 1 sin yi) (ft 1 cos y, — tPft 1 sin yi) 

= ft 2 # 1 cos 2 yj — u 2 u 1 tfiP sin 2 y t (C-26) 

u 2 C 02 u 1 C 20 — (ft 2 cos y 2 + u 2 u' sin y 2 ) (ft 1 cos y 2 — ft’ft 1 sin y 2 ) 

= tPC’ cos 2 y 2 — u^tPu 1 sin 2 y 2 (C-27) 


A 0 becomes 

A 0 = -WV + CmrR { u 1 (yf sin y, + y| sin y 2 ) 

— tPoftu 0 [u 2 (cos y t + COS y 2 ) + U 3 (sin yi + sin y 2 )] 

— ft 2 ™ 0 [u 3 (y i cos yt + y 2 cos y 2 ) — u 2 (yi sin y a + y 2 sin y 2 )] 

— (u 2 u l cos 2 y! — Wu'iPTi 1 sin 2 yi) (« 2 yi cos yj + t/ 3 yj sin y t ) 

— (tPu 1 COS 2 y 2 — u 2 u 1 u'u 1 sin 2 y 2 ) <3° (u 2 y 2 COS y 2 + U 3 y 2 sin y 2 ) 

— (U — U 2 ) W (y t cos y t + y 2 cos y 2 ) + ft 1 17 W (y , sin yi + y 2 sin y 2 ) } 

(C-28) 


Expansion and multiplication establish the identities 


tru 1 


0 

0 

0 


1 

0 

0 


0 

0 

0 


(C-29) 


Wu'u'u 1 — 


0 

0 

0 


1 

0 

0 


0 

0 

0 


= —u 2 u l 


(C-30) 


U - U 2 = 


1 

0 

0 


0 

0 

0 


o' 

0 

1 


= U 1 + U 3 


(C-31) 


frU 2 = 


0 

0 

0 


0 

0 

1 


0 

0 

0 


(C-32) 
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Although these substitutions would all simplify the computation of A 0 , the second 
and third alone provide 

A 0 = — 3°$ 00 a)° + 0 7irR { u 1 (yrsin y x + y| siny 2 ) 

— [« 2 (cos yj + cos y 2 ) + u 3 (sin y, + sin y 2 )] 

— U 2 a>° [u 3 (y t COS yi + y 2 COS y 2 ) — U 2 (y t sin yi + y 2 sin y 2 )] 

— tPtP 75° [ U 2 (y! COS y t + y 2 COS y 2 ) + M 3 (yj Sin yj + y 2 Sin y 2 )] 

-(C/ 1 + U 3 ) Z°U l (y! COS yj + y 2 COS y 2 ) + u'U^Tfu 1 (y t sin y! + y 2 sin y 2 ) } 

(C-33) 


Further reduction of A 0 requires expansion of <o° as o>° = [« t <u 2 w 3 ] T . This 
explicit representation provides 


U to to 


o>ia>3 a> 2 a>3 — (o>i + <o|) 

0 0 0 
(a>2 “I” <t)§) ““ 0)i0>2 


(C-34) 


U"r - 


ti>2 a>i 

0 0 

0 


0 
0 

(1)3 0) 2 


(C-35) 


w 2 u V = 


0>3 0 (»[ 

0 0 0 

0 0 0 


(C-36) 


( U 1 + U^^u 1 = 

' 1 0 0' 

0 0 0 

■ 0 

<03 



1 

0 0 

1 


0 0 1 

too __ 


L _a, 2 J 



"0 

0 

o' 


0 


" 0 

u'U'Wu 1 = 

0 

0 

0 


<03 

= J 

0 


0 

1 

0 


_ to 2 _ 


0J3 


to.iU 


so that A 0 becomes 

A 0 = — w°<t> 00 ft) 0 + QilrR | u 1 (yi sin y, + yi sin y 2 ) 


(C-37) 


(C-38) 


0) 2 W 3 

0 

(cos yj + COS y 2 ) + 

«i + l»2 

0 

tt) 1 to 2 


Wjft) 3 


(sin y] + sin y 2 ) 


+ 


) 

O O 

1 

(yi COS y! + y, COS y 2 ) — 

COi 

0 

L “2 J 


0)3 


(yi sin Yl + y 2 sin y 2 ) 
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0 

0 

0 


(yi sin y 2 + y 2 sin y 2 ) 



<*>i 

(y, COS yi + y 2 COS y 2 ) + 

0 

0 _ 



0 “ 


" 0 

+ 

Ct>3 

(y i COS y, + y 2 COS y 2 ) + 

0 


ti>2 


<*>3 


(yi sin yt + y 2 sin y 2 ) 


A 0 = — o> 0 4>°V > + CyilrR {u l (yi sin y t + yfsiny 2 ) 

+ (w,a 2 U J — la-MzU 1 ) (cos y t + COS y 2 ) + [(<o 2 + o>I) U 1 + u)iO) 3 « 3 ] (sin y t + sin y 2 ) 
+ 2o) 2 u' i (yj COS y! + y 2 COS y 2 )} (C-39) 


From Def. 41, we also find 

A 1 = -4> u y, (C'VC 01 + yjW 1 ) U 1 - [C 10 S°C 01 + yjU 1 ] 4> 1J [C 10 o>° + y^ 1 ] 

4- Oft [-Ru 2 C 12 (C 20 3°C 02 + y.u 1 )" Ru 2 - Ru 2 C 10 (S°S°) u 2 r 
+ (~ Uu 2T C 21 u 2 R 2 + C 12 u 2 u 2T R 2 ) y 2 (C’VC 01 + fcS 1 ) a 1 ] (C-40) 

and 


A 2 = (C 20 m°C 02 + y M l ) U l - (C 20 w°C 02 + y 2 U’) $ 22 (C 20 a>° + y^ 1 ) 

+ Qn [ — tt 2 C 21 (C l VC 01 + yiU 1 ) 2 R 2 u 2 - Ru 2 C 20 (W) u 2 r 


+ ( - Uu 2 ‘C 12 u-R 2 + C 2 'uhr T R 2 ) y, (C 20 5°C 02 + y,u l ) u 1 ] (C-41) 

These expressions are certainly no simpler than that found previously for A 0 , and 
we must expect that their expansion and reduction to expressions of minimum 
complexity would again be a substantial chore (although one that a computer can 
execute numerically in very little time, since only multiplications and additions 
are involved). 

It is becoming clear that the task of obtaining exact equations of motion for the 
case illustrated in Fig. C-l is much more arduous than for the simpler system 
analyzed in Appendix B (Fig. B-l), although both systems involve three bodies and 
five degrees of freedom. Having made this point, we now attack our second objec- 
tive, and seek an approximation of Eq. (C-8) in which y u y 2 , y,, y 2 , y,, and y 2 are 
small enough to warrant linearization in these variables. 

The matrix A 0 in Eq. (C-39) then simplifies to 


A 0 ^ — w°«l> 00 ( 1 ) 0 + QWrR {2 (<»i<o 2 u 3 — i^wsti 1 ) 

+ (yi + y 2 ) [(<01 + <o5) u 1 + <o 1 w 3 m 3 ] + 2 (-f] + y 2 ) o> 2 u 3 ) 
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For A 1 and A 2 , we can substitute the direction cosine approximations (see 
Eqs. C-10-C-12) 

C 01 a u + yi & C 10 S u - YiM 1 

C 02 S U + y 2 U 1 C 20 SiU — yjM 1 

C 12 = C 10 C 02 a C7 - nu 1 + ysU 1 
C 21 = U + yiU 1 — yjW 1 

into Eqs. (C-40) and (C-41). The desired approximations are 

A 1 — 4 >11 yioi°U 1 - S 0 ^ 11 ^ 0 + y 1 U v Z°^ 11 to 0 

- y^tTS 11 *) 0 ~ y itt 1 * 11 ® 0 +tf>* 11 yi2W’ 

+ 042 [ — U 2 (U — ytU 1 + y 2 u 1 )'5°<i> 0 R 2 ti 2 

— U 2 (<o°y 2 ti’ + y 2 U 1 <i) 0 — y 2 U 1 <5°u) 0 + 'S 0 a! 0 y 2 M 1 ) R 2 U 2 

- flu 2 (U - yjU 1 ) W« 2 r — R 2 (U — u 2 u 2T ) y 2 S°« l ] 

= -S 0 * 11 ,,. 0 - 04?R 2 u 2 3»°S°u 2 - Q«MW« ! 

- 7i [ < t> 11 (0 0 U 1 + u 1 $ 11 (0° + to 0 ^ 11 !! 1 ] 

- y 2 [Q7lR 2 u 2 (5°u l + «*5°) u 2 + 042fl 2 (t/ - u 2 u 2r ) W] 

+ yi [uW'»° - a>°U 1 4> ll m° + WW 

+ WiFiiWii 2 + 044flru 2 iF5 o 3 o u 2 ] 

+ yjI-WWWu 2 + Q7lR 2 u 2 u’Z%°u 2 - Q42R 2 uW , u i u 2 ] 

= -75°$ 11 <u 0 + 04Z (R 2 + Rr) SWu 2 ® 0 

+ y t [^“U 1 - U 1 ^ 11 + ($ 11 U 1 )~] c 0 

+ 72 [04?R 2 w 2 tF + 042R 2 u 2 u 1 u 2 + ORR 2 (1/ — u 2 u 2r ) u 1 ] o>° 

+ 7l [u 1 ®°<t> 11 - - 042 R (R + fJHW]/ 

+ y 2 [ORR 2 ^ 0 ^ 0 ] (C-44) 

and 

A 2 = -<t> 22 y 2 »°U 1 - <S 0 $ 2 V + y 2 U 1 <3°<& 22 (i, 0 

- y 2 a>°U 1 4? 22 <D 0 — 7 2 U 1 $ 22 m° + S°<t> 22 y 2 U 1 <t)° 

- •S°$ 22 7 2 u 1 + 042 [-u 2 (U + 7 xU 1 - y 2 U 1 )S°w°R 2 U 2 

- U 2 (-y^RoV + y^VTi 1 + yiUV + c5°y 1 S 1 ) R 2 U 2 
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- RrW ( U - t r V! )sW - R 2 (U - u 2 u 2T ) y^u 1 ] 

= ~7t°‘b 22 *> 0 - Q??R 2 uW > u 2 - 0 »fimWu ! 

- yi [QTlR 2 ^ 2 (s 0 !? + 2 ‘s 0 ) u 2 + Q7?R 2 (C/ — u 2 u 2T ) t^M 1 ] 

- y 2 [$ 2 W + M 1 $ 22 <d° + s 0 * 22 ^] 

+ yi [ -Q7?R 2 u 2 (uVs°u 2 - 5Wu ! + s^tTu 2 )] 

+ y 2 [tPs 0 * 22 ** 0 - (S°a 1 $ 22 o)° + 3°<l> 22 tiV + 077R 2 w 2 8 1 3 0 S°u 2 

+ Q'}lRm J U 1 m 0 7S°U 2 ] 

= -Z°<S> 22 0>° +Qyi(R 2 + Rr)Wu 2 m ° 

+ yi [Q7lR 2 tPu 3 + (ynRWtPu 2 + Qt?R 2 (C7 - tru 2r ) tP] o. 0 
+ y s [$ 22 tt 1 - «><i > 22 - (wn <■>“ + yi [QftR^VSV] 

+ y 2 [fiWfc 22 - 33°tP$ 22 + <3°$ 22 ir 

- QWR (R + r) tPtP® 0 !* 2 ] < 0 ° (C-45) 

The augmented body inertia matrices are available from Def. 36 as 
4 . 0 o = p + m(D 0lT D°'U - D ol D olT ) + m(D 02T D 02 U - D 02 D 02T ) 

= I 11 + m(D nr D u U - D u D llT ) + (071 - m) ( D lor D 10 U - D 10 D'° r ) 

= 722 4 . m (D 22r D 22 U ~ D 22 D 22T ) + (07? - m) (D 2 ° T D 20 U - D 20 D 2 ° r ) 


With Eqs. (C-7), these matrices become 

$00 = /« + mr 2 (U — U 2 ) + mr 2 (U - U 2 ) = 7° + 2mr 2 (f/ 1 + U 3 ) (C-46) 

$u _ jn + _ m )2 (7/ _ [72) + (Qw - m) R 2 (1/ - (/ 2 ) 

= + m - rn) ((7 - l/ 2 ) =2 Z“ + m - m) (^~j ( U 1 + f/ 2 ) 

and, since 


<f> 2 “ = <|) 11 

defining 


722 = JU = 



(U 1 + U 3 ) 



(4071 - 3m) ( U 1 + U 3 ) = (/ - OftR 2 ) (CR + U 3 ) 



(C-47) 


(C-48) 
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as the moment of inertia of a rod of mass m and length 2 L about an axis normal 
to the rod and passing through one end of the rod. 


The combination of Eqs. (C-44)-(C-48) is simplified by the following identities: 

[Ru 1 = 0 ; — tPL / 1 = (UW) T = 0 

Ihu 1 = u 1 ; Wu 1 = 0 

u~u 3 = U 3 iP ; u 2 ii' = — CRm 3 

U - U 2 U- T = U - U 2 = U 1 + U-- 
SPtPu 2 = 0 

u 1 (U - U 3 ) = u 1 ( U 1 + U 2 ) = tPf / 2 = U 3 u 1 

Now Eqs. (C-44) and (C-45) can be written as 
A 1 a* — (J - Q7lR 2 )m° ( U 1 + C/ 3 ) + CfaR ( R + r) «*3°SV 

4- (7 - Q f nR~) y, (t/ 3 u 1 - tT(7 3 + u 1 ) „» + 2QRR 2 y 2 t/ 3 tPo,° 

+ y, {(/ ~ OWR 2 ) + l/ 3 ) -tfu'U 3 + Z°U 3 u l ] 

- Qt/R (R + r) uru'zPTi 2 } o , 0 + Y MWtPtfu 3 * 0 

= - (7 ~ (YnR-)~ 0 (U l + C/ 3 ) <o° + Q 4 ?R (R + r) «%°u 2 «) 0 
+ 27y 1 [/ 3 M 1 u >° - 2Q4?R 2 (y, - y 2 ) [7 3 wV 

-I- y, {/ [tT«° ([/' + [/■'■) + < 5 ° (CRw 1 - U 1 U 3 )] CO 0 

- QtzR 2 [W (l / 1 + C/ 3 ) + 3 " ( U 'u l - tp( 7 3 ) - o>° 

T QtlRrU 'u' 7 ?'u 2 a a } + y.QRR^^tPo ) 0 

or, after expansion, 

A 1 ^ — (/ — Q 7 ?R 2 ) (oiiW.-jM 1 — C0i<i) 2 u 3 ) — OflR (R + r) (102UI3U 1 — a )](1> 2 U 3 ) 

+ 27y,o, 2 U 3 - 204? R 2 (y, - y 2 ) u) 2 U 3 

+ yi {/ [(w? — 0)1) U 1 + o>io> 3 W 3 ] — OftR 2 [(o)? — 0)3 + o)| + wj) U 1 + o>io> 3 « 3 ] 

~ QRRr (o>i + o> 2 ) u 1 } + y 2 Q4?R 2 [(»>'f + ulju 1 + u>iu> 3 u 3 ] 

= it’ { — (J + CR'lRr) < 1 ) 20)3 + [/ ( 0)1 — o) 2 ) — 04?R 2 (o; + o> 2 ) 

— Q 4 ?Rr (o)T + 0)1)] y! + QRR 2 (a>l + o>l) y 2 ) 


(C-49) 
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+ ft 3 {(7 + Q?lRr) o)! 0 ) 2 + 27o) 2 yi — 2Q47R 2 0) 2 (y a — y 2 ) 
+ (7 — Q74R 2 ) 0 )! 0 ) 3 yi + 09?R 2 o, 1 o, 3 y 2 } 
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and, by parallel construction, 


A 2 ^ u' { — (/ + QnRr) u) 2 k )3 + [7 (u >2 — <■>§) — 04ZR 2 (wi + <»I) — Q/JlRr (i»i + to 2 )] y 2 
+ QTzR 2 (mi + to 2 ) yx) + ti 3 {(/ + tyyiRr) a)ito 2 + 27o> 2 y 2 + 2Q7lR 2 o>2 (yi — f 2 ) 

+ (7 ~ Q4?R 2 ) (0 1 a> 3 y 2 + 0WR 2 it) l u) 3 y,} (C-50) 

Equations (C-43), (C-49), and (C-50) are in final form, ready to be substituted 
into Eq. (C-8) to obtain linearized equations of motion. These equations also 
require the appropriate approximations for the coefficients on the left side of 
Eq. (C-8), as obtained from Defs. 38-40. 

From Def. 38, the linear approximation of a 0 „ is given by 

a 00 =* $0° + $ 11 + yd?* 11 - y^ 11 !! 1 + $ 22 + y 2 tP$ 22 - y 2 $ 22 U l 

- Q?Z {(D 10T D°'U - D l0 D 0lT ) + (D 2 ° T D 02 U - D 20 D 02T ) 

- y, (D 10T u'D 01 U + u 1 D 10 D 0lT ) - r 2 (D 20T u 1 D 02 U + u l D 20 D 02T ) 

+ (D°' r D 10 U - D 01 D 10T ) (C7 - yi tP) + yi u'D 01 T D 10 + y 1 D 0,r i?D’°l7 
4- (D 2 ' T D‘-U - D 21 D 12T ) (17 - yitl 1 ) + y^D 211 © 12 + ( Yl - y 2 ) D 2lT u’D 12 U 

- y m'D 21 D' 2T + (D 02T D 20 (7 - D 02 D 20T ) (C - y 2 u l ) 

+ y 2 u'D 02T D 20 + y 2 D 02 T u 1 D 20 C + (D 12T D 21 l) - D 12 D 2lT ) (U - y,u>) 

+ y 2 S'D l2T D 21 - (y, - y 2 ) D 12T Tl 1 D 21 U - y 1 U 1 D 12 D 2lT ) 

Substituting Eqs. (C-47) and (C-7) simplifies this lengthy expression somewhat, 
providing 

floo = * 00 + 4> u + <I> 22 + (7 - Q%R 2 ) (ti'U 3 - 17 3 C 1 ) (y a + y 2 ) 

-CmR{-2r(U - U 2 ) + r( Vl + y 2 )u 1 t/ 2 
-r(U - U 2 ) (17 - ydT) - ry^ 1 

- R (77 - C7 2 ) (17 - yjU 1 ) - Ryjtt 1 + Ry 2 tP[7 2 

- r{U — U 2 ) (U - y 2 u 1 ) - r y2 tP 

— R(U — U 2 ) (U - y 2 u l ) - Ry 2 U l + RyjtFf/ 2 ) 

= <t> 00 + + <t> 22 + /(m 1 ^ 3 - U 3 U 1 )(y 1 + y 2 ) 

+ CftlR 2 {(U 3 iP - iPU 3 ) ( yi + y 2 ) + 2 (C7 1 + C7 3 ) 

— C 3 ?? ( Yl + y 2 ) — (yi + y 2 ) fTl7 2 + (y a + y 2 ) U 1 } 

+ Q7?Rr {2 (U 1 + l/ 3 ) - ( 7l + y 2 ) tPt7 2 + 2 (17 1 + U 3 ) 
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— XJW (yi + y 2 ) + (yi + y 2 ) U 1 } 

= $00 + $11 + $22 + 7(wi[/ 3 - um 1 ) ( yi + y 2 ) 

+ c ynR 2 [2 U 1 + 2U 3 + ( yi + y 2 ) (u 1 - u'U 2 - WU 3 )] 

+ tynRr [4C7 1 + 4 U 3 + ( yi + y 2 ) (E 1 - E 1 !/ 2 - UW)] (C-51) 


By expanding the constituents of this expression, one can verify that 

$0° + $ 1J + $ 22 4- (2 (ynR 2 + 4(YhrR) ( U 1 + U 3 ) 

= 1° + 2 (J + 2QfiRr + mr 2 ) (U 1 + U 3 ) 

~ 1° + ^2mr- + 4 mrL + 

where the new symbol I is the inertia matrix of the total system with respect to the 
system mass center when yx = y 2 = 0. Thus, we have 

a 00 = f + (yx + y 2 ) [] (u'U 3 - U 3 ?P) + CMR 2 (S 1 - E 1 !/ 2 - E 1 !/ 3 ) 

+ QflRr (a 1 - iPU 2 - UW)] 


■| mL 2 ^ (U 1 + U») = I (C-52) 


But 


u 1 - S 1 !/ 2 = tPU 3 


so we have 

a 00 ^ I + (yx + y 2 ) (/ + OERr) (SW 3 - UW) (C-53) 

Expanding Z, /, and R from Eqs. (C-52), (C-48), and (C-3) provides 


a 00 Z° + 2m 


^r 2 + 2rL + L 2 ^ ( U l + U 3 ) + m^-| L 2 + rL ^ (E 1 !/ 3 -IPE 1 ) ( 7l + y 2 ) 


Z? + 2m^r 2 +2rL + -|-L 2 ^ 
0 
0 



(C-54) 


It may be noted that a 00 is the inertia matrix in basis {b°} of the total system 
relative to the system mass center. Detailed construction of this matrix from physi- 
cal definitions confirms that Eq. (C-54) provides the correct linearized approxima- 
tion of a 00 . 
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Construction of the equations of motion from Eq. (C-8) requires expansion of 
the linearized approximations of a 10 and a 20 from Def. 39. The expansion provides 


a 10 s u 1 ^ 11 (17 - y 1 u 1 ) - Qnu lT [D 10 T D 01 U - D’°D 0lT 

- y 1 D 1 ° r u 1 D 01 l7 - y 1 D l ° T D 01 TP + D 12 T D 21 U - D 12 D 2lT 

- ( Vl - y 2 ) D' 2 T WD^U - yi D 1JT D 21 « 1 + y 2 D 2 lT D 1 W] 


Equation (C-47) yields 

u 1 "* 11 (17 - 7l u l ) = (/ - OhR 2 ) u lT (t/ 1 + L7 3 ) (17 - y^ 1 ) = (/ - OftR 2 ) u ir 
Equation (C-7) yields 

u lT D 10 D 01 = u lT D 12 D 21 = u lT u 1 = 0 D 10 VD 01 = D 12T u 1 D 21 = 0 

Thus, a 10 simplifies* to 

a 10 =* (/ - OftR 2 ) « ir - (ftl (D loT D 01 + D 12T D 21 ) u lT 


= (]- VhR 2 + Q4?R 2 + Q nRr) u lT = (/ + Q??Rr) u lT 
= m^|-L 2 + rL^u 17 (C-55) 

Similar calculations produce 

L 2 + rL^u lT (C-56) 

For the remainder of the terms in Eq. (C-8), only constant-term approximations 
are required. Defs. 39, 40 provide 



002 = 0oi = u'm L 2 + 
a 12 = a 2 1 » -(Ynu lT (D 2lT D 12 U - D 21 D 12T )u 1 
= ^R 2 = ^mL 2 

a 11 ^u lT & 11 u 1 = (J — QflR 2 ) u lT (U 1 + U 3 ) u 1 
/ 1 rn \ 

~ J — (YftR 2 = mL 2 

a T , s* u 1 T $ 22 u 1 = / - Q??R 2 = mL 2 Q- - 


(C-57) 

(C-58) 

(C-59) 

(C-60) 


*Note that in this example, the quantities S l() and a 2u (defined by Eq. 37b and available as the 
variable parts of the linear approximations shown here for a 10 and a 20 ) are zero. This convenient 
result stems from the diagonality of <t‘“ and the parallelism of D 10 , D 01 , D 12 , and D 21 ; it is not a 
general result. 
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Equations (C-54)-(C-60) provide all of the coefficients on the left side of the 
equations of motion (Eq. C-8), and all of the terms on the right are available from 
Eqs. (C-49), (C-50), and (C-43), utilizing also the expressions for C 01 and C 02 . The 
resulting linearized scalar equations are obtained, after simplification of terms, 
as follows: 

h<Oi + (7 + (mRr) (y, + y 2 ) = (h ~ I3 ) 0 * 2^3 + (/ + 0 VlRr) (0 )\ — <o|) (y, + y 2 ) 

(C-61) 

7 2 a > 2 — (7 + CytlRr) (y, + y 2 ) = (I 3 — /,) <o 3 (Oi — (7 + (YHRr) o) 1 o ) 2 (yi + y») 

(C-62) 

h&a — (/ + Q VlRr) (y, + y 2 ) a > 2 = (Ii ~ h) "io >2 

+ 2 (J + C yyiRr ) ai 2 (yi + y 2 ) 4- (Q HRr + }) <o,<i > 3 (y, + y 2 ) (C-63) 

(7 + QfftRr) d>i + (/ — QfiR 2 ) y 1 + Q^R 2 y 2 = — (/ + Q YlRr) to 2 a > 3 

+ [7 (<a! — «* 3 ) ~ QyiR 2 (<af + «|) — QtylRr (caf + ml)] y, + OliR 2 («>! + «i) y 2 ~~ fcyi 

(C-64) 

(7 + 0 TlRr) ax + QWR 2 y‘i + (7 ~ 0 TiR 2 ) y 2 = — (7 + 0 flRr) u> 2 o > 3 

+ [7 (o>2 — »f) ~ QllR 2 (o>i + a> 2 ) — (yplRr (o)l + (a 2 )] y 2 + QT?R z («, + u> 2 ) yi ~~ ky 2 

(C-65) 

where /i, L>, and Z 3 are the diagonal elements of I as defined by Eq. (C-52), so that 
these are the principal axis moments of inertia of the total system when yi = y 2 = 0. 

Note that in Eqs. (C-64) and (C-65) the coefficients of y, and y 2 are interchanged, 
so that their coefficient matrix is symmetric. This result was anticipated in Sec- 
tion IIIC, since the system in Fig. C-l is an example of a “meridional deformation 
case” (see Fig. 11). 

Equations (C-61)-(C-65) have been derived by the laborious process of lineariz- 
ing Eq. (1) from Section II A of this report. In work not shown here, the same 
equations have also been derived from first principles. Equations (C-61)-(C-63) 
may be obtained simply by taking the inertial frame time derivative of the system 
angular momentum for the system mass center. Equation (C-64) is available from 
the dot-product of b, and the equation 

M^ 1 = H^ 1 + me, X P, 

where is the external moment about /«, applied to <f,, H ^ is the inertial angular 

momentum of referred to /«,, c, = Lb), and P, is the vector from the system 
mass center to /i,. Equation (C-65) is available from a similar calculation. The 
derivation of Eqs. (C-61)-(C-65) from first principles is a much more rapid and 

jpl Technical report 32-1565 


102 



more error-free process for an experienced analyst than is the application of Eq. (1) 
for this purpose, although at present a digital computer can accomplish only the 
latter task. 

Although this appendix has served to provide a second explicit check on the 
linear terms in Eq. (1), its primary objective is to check these results against the 
generic linearized equations recorded as Eqs. (30) of Section IIIB. 


For a three-body system with only small relative motions, Eq. (30a) becomes 
(fl 0 o + L) + floiy'i 4- a 02 y\ = A 0 + A 0 -f ( C 01 + C 01 ) A 1 + C 01 A ] 

+ (C 02 + 6 02 ) A 2 + C 02 A 2 (C-66) 

and Eq. (30b) provides 

(a 10 + a 10 ) <i° + + a 12 y 2 — g lT (A 1 + A 1 ) + f t + (C-67) 

and 

(a 2 o + a 20 ) <i)° + a 21 v\ + a 22 y2 — g 2 ^ (A 2 + A 2 ) + r 2 + £ 2 (C-68) 

Having struggled to derive Eqs. (C-61)-(C-65) from Eq. (1), we can now recog- 
nize that we have already executed the steps implied by Eqs. (C-66)-(C-68). 
Equations (C-53) and (C-54) are, in fact, expressions for floo A flon, and Eqs. (C-57) 
and (C-56) provide a 01 , ~a 0 2 , and a 2 o, with the further indication that a 10 = a 20 = 0. 
(This simplification is computationally significant, and its validity in every new 
case should be investigated.) The right-hand sides of Eqs. (C-66)-(C-68) apparently 
also coincide with results already obtained, with A 0 + A 0 given by Eq. (C-43), 
A 1 + A 1 given by Eq. (C-49), and A 2 + A 2 given by Eq. (C-50), and with C 10 = 
C zo = U, C'° = — yiU 1 , and C 20 = —y 2 tP. 

A These results also match those available from Eqs. (33a) and (33b) for C ri and 
C ri (/ = 0 and r = 1, 2), as well^as Eqs. (36) and (37) for the quantities cfon, ^oo, Ooj> 
Am, (fn-j, a„ 2 , a) 2 = a n , a n , a 22 , A 0 , 4°, A 1 , A 1 , A 2 , and A 2 . Thus, the generic partially 
linearized equations in Section IIIB have also been confirmed. 

The final step in the development of hybrid-coordinate equations is the trans- 
formation of the small variables (yi and y 2 ) into a new set of variables 17, and rj 2 , 
and (when justified) the truncation of the new set of variables to obtain for the 
total system a set of equations of reduced dimension. As noted in Sections IIIC 
and D, the quest for an appropriate transformation is, in general, a complex com- 
putational task involving eigenvector calculations. For the simple system under 
consideration in this appendix, however, the necessary transformation can be deter- 
mined by inspection of Eqs. (C-61)-(C-65), or by physical argument. 

We require the transformation 


1 

1 

1 

" 1 

i—H 


1 

£ 

L 

1 

[ 

~ 2 

1 

1 

t— • 

1 


L ^ J 


(C-69a) 
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with the inverse transformation 


Vi 1 1 yi 

= (C-69b) 

_ V2 j [_ 1 — 1 J |_ ya _ 

We can replace yi + y 2 in Eqs. (C-61)-(C-63) by t) U and replace Eqs. (C-64) and 
(C-65) by their sums and differences, to obtain 

Zi<«i + (7 + QiiRr) rj\ = ( I 2 — h) o . 2 <» 3 + (7 + QflRr) (o>| — <al) jp (C-70) 

7 2 *2 — (7 + (mRr) rj i<o 3 = (I a ~ 7,) (o 3 (ui — (7 + Q 7lRr) <ui<o 2 7 ?! (C-71) 

Z 3 *3 — (7 + QilRr) r^iib 2 — {1 1 — I 2 ) t»i<i ) 2 + 2 (7 + QTIRt) io 2 r)i + (7 + Q flRr) < 010 . 31 ?! 

(C-72) 

Jrji + 2 (7 + Q'HRr) <oi = —2 (7 + Q TlRr) <o 2 a> 3 + [7 (o> 2 — o> 3 ) — (YnRr(<al + o>§) — k ] 77 x 

(C-73) 

(7 - 2QMR 2 ) ij 2 = [7 (<02 - «.!) - 2QMR 2 (<of + <ol) - CynRr («.? + «,§) - fc] i ? 2 (C-74) 

Because only the final equation involves rj 2 , and this equation is satisfied by 
i? 2 ==0, it is quite satisfactory for the purposes of dynamic analysis of rotational 
motions of the central body (defined by <»,, <o 2 , o> 3 ) to truncate this system of equa- 
tions, abandoning Eq. (C-74) entirely. 
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